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EXECUTIVE  SUMMARY 


Over  the  last  three  decades,  there  have  been  continuous  efforts  in  developing  failure  criteria  for 
unidirectional  fiber  composites  and  their  laminates.  Currently,  there  exist  a  large  number  of 
lamina  failure  criteria  and  laminate  failure  analysis  methods.  In  this  project,  a  comprehensive 
and  objective  study  of  lamina  and  laminate  failure  criteria  was  performed.  Comparisons  among 
the  commonly  used  failure  criteria  were  made  for  failure  in  unidirectional  composites  under 
various  loading  cases.  From  these  comparisons,  the  characteristics  of  these  criteria  were 
identified  and  discussed.  Further,  with  the  aid  of  some  limited  experimental  lamina  and  laminate 
strength  data  available  in  the  literature  and  new  data  generated  by  the  authors,  an  attempt  was 
made  to  select  the  failure  criteria  and  laminate  analysis  methods  that  are  mechanistically  sound 
and  are  capable  of  accurately  predicting  lamina  and  laminate  strengths  for  states  of  combined 
stresses.  It  was  found  that  those  lamina  failure  criteria  which  separate  fiber  and  matrix  failure 
modes  most  accurately  predict  lamina  and  laminate  strength. 
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1.  INTRODUCTION. 


Over  the  last  three  decades,  there  have  been  continuous  efforts  in  developing  failure  criteria  for 
unidirectional  fiber  composites  and  their  laminates.  Currently,  there  exist  a  large  number  of 
lamina  failure  criteria  and  laminate  failure  analysis  methods.  A  comprehensive  evaluation  of  the 
accuracy  of  these  failure  criteria  in  the  light  of  available  experimental  data  seems  to  be  overdue. 

There  are  two  major  elements  in  the  analysis  of  composite  laminates,  i.e.,  lamina  failure  criteria 
and  laminate  stress  analysis  with  lamina  stiffness  reduction.  Between  the  two,  the  accuracy  of 
the  failure  criterion  is  the  most  crucial  issue. 

Evaluating  these  lamina  failure  criteria  is  a  two  part  process.  The  first  step  is  to  characterize  the 
criteria  in  their  ability  to  predict  failure  in  a  unidirectional  composite  or  a  lamina.  These  are  the 
precise  conditions  for  which  the  criteria  were  designed.  Those  criteria  which  correlate  with 
experimental  data  and  those  criteria  which  are  mechanistically  sound  can  be  identified. 

Secondly,  the  lamina  failure  criteria  must  be  evaluated  in  their  ability  to  predict  the  failure 
strength  of  a  laminate  comprised  of  laminae  with  varying  fiber  orientations.  Endorsing  a  lamina 
failure  criteria  based  on  its  success  with  unidirectional  failure  predictions  is  premature.  In  a 
laminate,  failure  mechanisms  are  more  complicated  (i.e.,  in  situ  laminae  can  exhibit  considerably 
higher  matrix  strength  than  experimentally  determined  through  unidirectional  lamina  tests).  A 
lamina  failure  criterion  must  be  flexible  and  accommodate  the  more  complicated  nature  of 
laminate  analysis. 

In  this  study,  six  failure  criteria  which  appear  representative  of  those  that  have  been  proposed 
over  the  years  are  investigated.  These  failure  criteria  which  appear  representative  of  most  of 
those  which  have  been  proposed  over  the  years  are  maximum  stress,  maximum  strain,  Hill-Tsai, 
Tsai-Wu,  Hashin-Rotem,  and  Hashin  criteria.  Maximum  stress  and  maximum  strain  criteria 
assume  no  stress  interaction.  Hill-Tsai  and  Tsai-Wu  criteria  include  full  stress  interaction. 
Hashin-Rotem  and  Hashin  criteria  involve  partial  stress  interaction.  Existing  lamina  and 
laminate  strength  data  are  used  to  evaluate  these  failure  criteria.  For  some  laminates  under 
certain  loading  conditions,  all  six  criteria  may  predict  similar  results,  and  their  performance 
cannot  be  ranked.  Therefore,  a  number  of  laminates  are  identified  for  which  the  strength 
predictions  according  to  these  six  criteria  are  substantially  different.  Off-axis  coupon  specimens 
were  cut  from  those  laminates  and  tested  in  uniaxial  tension.  Adhesive  film  was  placed  along  all 
the  interfaces  of  the  laminae  to  suppress  failure  due  to  free  edge  stresses.  To  avoid  complications 
arising  from  extension-shear  coupling  in  some  of  the  off-axis  specimens,  special  oblique  end  tabs 
were  used.  These  additional  strength  data  are  used  to  help  rank  the  six  strength  criteria. 

2.  LAMINA  FAILURE  ANALYSIS. 

The  purpose  of  the  lamina  failure  criterion  is  to  determine  the  strength  and  mode  of  failure  of  a 
unidirectional  composite  or  lamina  in  a  state  of  combined  stress.  All  the  existing  lamina  failure 
criteria  are  basically  phenomenological  in  which  detailed  failure  processes  are  not  described 
(macromechanical).  Further,  they  are  all  based  on  linear  elastic  analysis.  Nahast  [1]  and 
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Labossiere  and  Neal  [2]  have  made  an  extensive  literature  survey  of  existing  lamina  failure 
criteria  for  composites. 

A  list  of  the  lamina  failure  criteria  taken  from  references  1  and  2  is  presented  in  appendix  A. 

The  majority  of  the  lamina  failure  criteria  were  developed  for  two-dimensional  stress  states  in 
orthotropic  materials.  Some  of  the  criteria,  such  as  the  Tsai-Wu  criterion  which  is  a  completely 
general  tensor  polynomial  failure  equation,  have  reduced  forms  in  order  to  utilize  two  strength 
properties  for  two-dimensional  stress  states.  In  this  study,  only  such  2-D  criteria  are  included. 
The  in-plane  principal  strengths  in  a  composite  system  are  denoted  as  follows; 

•  X&X:  tensile  and  compressive  strengths,  respectively,  in  fiber  direction. 

•  7  &  F:  tensile  and  compressive  strengths,  respectively,  in  transverse  direction 
(perpendicular  to  fibers). 

•  S:  shear  strength 

For  a  strain  based  analysis,  the  corresponding  failure  strains  are  X^,  X^',  Fg ,  Y^',  and  Se  . 

The  ability  of  a  lamina  failure  criterion  to  determine  mode  of  failure  is  essential  in  bringing  this 
analysis  tool  to  the  laminate  level  (an  individual  lamina  failure  within  a  laminate  does  not 
necessarily  constitute  ultimate  failure).  Modes  of  failure  are  defined  as 

•  Fiber  Breakage  (mode  1):  longitudinal  stress  (Gn)  or  longitudinal  strain  (ejj)  dominates 
lamina  failure. 

•  Transverse  Matrix  Cracking  (mode  2):  transverse  stress  {Ozi)  or  transverse  strain  (£22) 
dominates  lamina  failure. 

•  Shear  Matrix  Cracking  (mode  3):  shear  stress  (X72)  or  shear  strain  (y/2)  dominates  lamina 
failure. 

It  is  important  to  point  out  that  both  mode  2  and  mode  3  are  matrix  failures.  The  two  modes  are 
separated  because  they  are  caused  by  different  stress  components  according  to  some  criteria.  For 
example,  according  to  the  Maximum  Stress  Criterion,  mode  2  should  be  interpreted  as  matrix 
cracking  due  to  G22,  and  mode  3  should  be  interpreted  as  matrix  cracking  due  to  X72. 

2.1  LAMINA  FAILURE  CRITERIA. 

Lamina  failure  criteria  can  be  categorized  into  three  groups. 

•  Limit  Criteria:  These  criteria  predict  failure  load  and  mode  by  comparing  lamina  stresses 
<^11 ,  CT22,  and  Xi2  (or  strains  Zu ,  £22 ,  and  772 )  with  corresponding  strengths  separately.  Interaction 
among  the  stresses  (or  strains)  is  not  considered. 
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•  Interactive  Criteria:  These  criteria  predict  the  failure  load  by  using  a  single  quadratic  or 
higher  order  polynomial  equation  involving  all  stress  (or  strain)  components.  Failure  is  assumed 
when  the  equation  is  satisfied.  The  mode  of  failure  is  determined  indirectly  by  comparing  the 
stress/strength  ratios. 

•  Separate  Mode  Criteria:  These  criteria  separate  the  matrix  failure  criterion  from  the  fiber 
failure  criterion.  The  equations  can  be  dependent  on  either  one  or  more  stress  components; 
therefore,  stress  interaction  varies  from  criterion  to  criterion  within  this  group.  If  the  failure 
equation  contains  only  one  stress  component,  then  the  failure  mode  corresponds  to  that  particular 
direction;  otherwise,  the  failure  mode  can  be  determined  as  is  done  with  the  interactive  criteria  by 
comparing  stress/strength  ratios  of  the  satisfied  equation. 

In  the  interest  of  keeping  the  project  manageable,  it  was  necessary  to  choose  only  a  representative 
subset  of  the  criteria  listed  in  appendix  A.  In  total,  six  lamina  failure  criteria  were  selected  to  be 
examined  in  further  detail.  From  the  limit  criteria  group.  Maximum  Stress  and  Maximum  Strain 
were  chosen.  From  the  interactive  criteria  group,  Hill-Tsai  (also  called  Tsai-Hill  or  Azzi-Hill) 
and  Tsai-Wu  were  chosen.  In  an  AIAA  Failure  Criteria  Survey  [3],  80%  of  the  respondents  said 
they  utilized  one  of  these  four  lamina  failure  criteria.  Figure  1  shows  the  breakdown  for  each 
criterion.  Maximum  Strain  is  most  commonly  used  at  30%  with  Maximum  Stress  next  at  22%. 
Hill-Tsai  and  Tsai-Wu  usage  came  in  at  17%  and  12%  respectively.  The  popularity  of  these  four 
criteria  and  the  fact  that  they  are  the  most  generalized  and  representative  of  their  respective 
groups  was  the  basis  for  their  inclusion. 

The  final  two  criteria  come  from  the  separate  mode  group.  Both  Hashin  and  Hashin-Rotem 
criteria  provide  for  separate  treatment  of  matrix  and  fiber  failure  modes  while  maintaining  some 
degree  of  stress  interaction  for  the  individual  modes. 

The  six  lamina  failure  criteria  which  are  considered  are: 

•  Limit  Criteria: 

Maximum  Stress: 

=  1  fiber  failure 

Jx 


—  =1  transverse  matrix  cracking 
Y 


^  72 

=  1  shear  matrix  cracking 


(1) 
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Maximum  Strain: 


1  fiber  failure 


=  1  transverse  matrix  cracking 


=  1  shear  matrix  cracking 


Interactive  Criteria: 


Hill-Tsai: 


Tsai-Wu: 

11  ^2^  22  ^iP  +  F22<3  22^  +  ^^12^  22  ^66'^  12^  ~  ^ 
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Separate  Mode  Criteria: 


Hashin-Rotem: 


=  1  fiber  failure 

Jx. 
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— j  +f— — j  =1  matrix  failure 
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Hashin: 
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A  u 


^7  =  1  fiber  failure  (compression) 


+(~^)  =-?  matrix  failure 
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Lamina  Failure  Criterion 

FIGURE  1.  RESULTS  OF  AIAA  FAILURE  CRITERIA  SURVEY 


For  Maximum  Stress,  Maximum  Strain,  Hill-Tsai,  and  Hashin-Rotem,  the  criterion  is 
generalized  for  either  tensile  or  compressive  stresses;  the  corresponding  (tensile  or  compressive) 
strength  value  must  be  chosen  based  on  the  sign  of  the  applied  stress.  The  Tsai-Wu  criterion  is 
designed  for  use  in  all  quadrants  of  the  stress  plane;  thus,  it  may  be  used  directly  without 
modification  for  different  stress  signs.  The  Tsai-Wu  criterion  requires  a  biaxial  test  to 
experimentally  determine  the  interaction  term  Fn-  It  has  been  suggested  to  use  F  12  =  \li2XX'), 
which  reduces  Tsai-Wu  down  to  the  Hoffman  criterion.  Narayanaswami  and  Adelman  [4]  found 
this  term  to  be  insignificant  for  the  most  part,  and  suggested  setting  it  equal  to  zero.  Cui  et  al. 
[5]  also  found  that  F12  =  0  gave  adequate  accuracy  for  engineering  purposes.  Thus,  to  avoid 
ambiguity,  Fn  is  set  equal  to  zero  in  the  present  study. 

The  Hashin  criterion  listed  here  is  a  slight  modification  of  the  2-D  criterion  presented  in  his  1980 
paper  [6].  In  that  paper,  Hashin  suggested  using  a  combination  of  both  axial  and  transverse  shear 
strengths  Sp,  and  5t  for  the  compressive  matrix  equation.  Since  it  is  difficult  to  find  transverse 
shear  strength  values  in  the  literature,  the  tensile  equation  given  is  used  as  the  compressive 
equation  by  simply  replacing  Y  with  Y'. 
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2.2  COMPARISON  AMONG  LAMINA  FAILURE  CRTTERTA. 


The  theoretical  comparisons  of  the  six  lamina  failure  criteria  discussed  in  the  following 
paragraphs  allow  each  criterion  to  be  evaluated  on  the  basis  of  purely  mechanistic  reasoning. 
Separating  the  theoretical  comparisons  from  the  correlation  with  experimental  data  helps  focus 
on  objectively  evaluating  just  the  mechanics  of  the  criteria.  The  ability  to  accurately  predict  data 
is  addressed  in  later  sections. 

A  FORTRAN  program  was  written  in  order  to  efficiently  analyze  all  six  lamina  failure  criteria 
under  all  possible  loading  conditions.  The  code  was  designed  to  be  flexible.  It  allows  for  an 
expandable  database  of  different  material  properties.  Laminate  configurations  and  loading 
conditions  are  easily  manipulated.  The  user  has  a  variety  of  output  files  to  choose  from.  Ready- 
to-plot  output  files  allow  the  user  to  immediately  plot  and  interpret  the  results.  The  code  together 
with  representative  examples  is  included  in  appendix  B. 

2.2. 1  Bidirectional  Stress  Plane. 

A  series  of  failure  envelopes  for  combined  stresses  is  presented  to  graphically  show  the 
characteristics  of  the  six  selected  lamina  failure  criteria.  These  envelopes  are  composed  of 
failure  stresses  normalized  by  the  lamina’s  respective  tensile  strengths  X  and  F  or  shear 
strength  5.  For  these  graphs,  the  material  properties  of  the  AS4/3501-6  graphite/epoxy  system 
tested  by  Sun  and  Zhou  [7]  is  used.  Table  1  lists  elastic  and  strength  constants. 

TABLE  1.  MODULI  AND  STRENGTH  VALUES  FOR  THE  AS4/3501-6  GRAPHITE- 
EPOXY  SYSTEM  [7] 


El 

138.90  GPa 

X 

2206.0  MPa 

E2 

9.86  GPa 

X' 

-2013.0  MPa 

Gi2 

5.24  GPa 

Y 

56.5  MPa 

V/2 

0.30 

Y' 

-206.8  MPa 

Ply  Thickness: 

0.132  mm 

S 

110.3  MPa 

Figure  2  is  a  plot  of  the  selected  criteria  in  a  (Su  -  G22  stress  plane  {ii2  —  0).  The  Maximum 
Stress  envelope  is  a  simple  rectangle  bounded  by  the  failure  loads  ±a/;  and  ±022-  Again,  because 
these  loads  are  normalized  with  X  and  Y,  compressive  Gn  (quadrants  n  and  HI),  and  compressive 
G22  (quadrants  HI  and  IV)  failure  segments  do  not  equal  unity  (i.e.,  X  >  1X^1  and  F  <  IF1  for  this 
case). 

For  the  analysis  using  the  Maximum  Strain  criterion,  failure  strains  were  calculated  from  the 
strength  parameters  using  a  linear  relationship: 
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-  -  -  -  Hill-Tsai 
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FIGURE  2.  COMPARISON  OF  LAMINA  FAILURE  CRITERIA  UNDER  Cji  -  C22 
BIAXIAL  STRESS 

The  Maximum  Strain  envelope  is  close  to  that  of  the  Maximum  Stress  but  is  slightly  skewed  due 
to  the  effect  of  Poisson’s  ratio.  There  is  considerably  more  skewing  in  the  vertical  {G22)  direction 
because  V12  »  \2i  in  unidirectional  fiber  composites. 

Both  the  Hill-Tsai  and  Tsai-Wu  criteria  allow  quadratic  stress  interactions;  therefore,  each  has  a 
curved  failure  envelope.  Both  of  these  criteria  match  up  with  the  two  limit  criteria  for  all  four 
unidirectional  loading  cases  {±Gn  with  G22  =  0  and  ±^22  with  Gu  =  0)  as  expected.  The  Tsai-Wu 
criterion  is  a  continuous  curve  throughout  all  four  quadrants.  The  only  parameters  that  vary  are 
the  stress  terms.  The  Tsai-Wu  criterion  includes  linear  stress  terms.  Hill-Tsai,  on  the  other  hand, 
is  a  purely  quadratic  criterion.  In  order  to  account  for  differences  in  tensile  and  compressive 
strengths  commonly  found  in  fiber  composites,  this  criterion  uses  the  appropriate  strength  values 
in  each  quadrant  (X  or  X'  and  Y  or  Y'  accordingly).  Though  both  are  interactive,  Tsai-Wu  and 
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Hill-Tsai  produce  different  failure  envelopes  in  the  stress  plane.  In  the  compressive  O22 
quadrants,  the  Tsai-Wu  failure  envelope  extends  beyond  the  longitudinal  strengths  X  and  X'. 

Finally,  both  the  Hashin  and  Hashin-Rotem  criteria  reduce  to  the  Maximum  Stress  criterion  in 
the  Gi/  -  G22  plane  since  T;2  =  0. 

A  plot  of  the  selected  criteria  in  a  Gh  -  Xn  stress  plane  {G22  =  0)  is  shown  in  figure  3.  The 
Maximum  Stress  envelope  in  this  stress  plane  is  again  a  rectangle,  bounded  by  the  failure  loads 
±Gii  normalized  by  X  and  ±Xi2  normalized  by  S.  The  Maximum  Strain  criterion  predicts  exactly 
the  same  loads  as  the  Maximum  Stress  criterion.  The  Hashin-Rotem  criterion  also  reduces  to  the 
Maximum  Stress  criterion  in  this  stress  plane.  Again,  Hill-Tsai  and  Tsai-Wu  failure  envelopes 
intersect  the  other  three  criteria  for  the  four  unidirectional  loading  cases  i±Gii  with  X12  =  0  and 
±ii2  with  Gu  =  0).  In  the  biaxial  loading  regions,  the  two  interactive  criteria  are  nearly  identical. 


•  EiU-Tsai 

—  Moimum  Strew  &  Stnun,  Hwhin-Rotem 


—  —  Tiai-Wu 


Hashin 


FIGURE  3.  COMPARISON  OF  LAMINA  FAILURE  CRITERIA  UNDER  Gn  -  Xj2 
BIAXIAL  STRESS 
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The  linear  stress  term  Gy;  in  Tsai-Wu  produces  a  slightly  higher  or  lower  failure  load  than  in 
Hill-  Tsai,  depending  on  the  quadrant.  Interestingly,  the  Hashin  criterion  reduces  to  Maximum 
Stress  with  a  compressive  (5n  and  reduces  to  Hill-Tsai  with  tensile  Gyy. 

Figure  4  contains  a  plot  of  the  six  criteria  in  the  022  -'^12  stress  plane  (G/y  =  0).  Maximum  Stress 
and  Maximum  Strain  are  identical  rectangles  showing  ±<522  normalized  by  Y  and  ±Xy2  normalized 
by  S.  Tsai-Wu  and  Hill-Tsai  produce  curved  envelopes  due  to  022-112  interaction.  Again,  the 
linear  G22  term  in  Tsai-Wu  produces  a  different  shape  than  Hill-Tsai,  pushing  the  failure 
envelope  beyond  the  lamina  shear  strength  S.  Hashin  and  Hashin-Rotem  in  this  stress  plane 
match  Hill-Tsai  exactly  considering  Ojj  =  0.  As  expected,  all  six  criteria  intersect  at  the  four 
unidirectional  loading  cases  (±G22  with  Xj2  =  0  and  ±Xy2  with  G22  =  0). 


Hill-Tsai,  Hashin-Rotem,  Hashin 
Maximum  Stress  &  Strain 


2 

1.5 

1 

0.5 

m 

'"2  0 

-0.5 

-1 

-1.5 

-2 

^-3-2-1012 


HGURE  4.  COMPARISON  OF  LAMINA  FAILURE  CRITERIA  UNDER  G22  -  Xy2 
BIAXIAL  STRESS 
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2.2.2  Off-Axis  Loading. 


The  six  lamina  failure  criteria  may  also  be  characterized  by  their  failure  predictions  of  a  lamina 
subjected  to  off-axis  loading.  Figure  5  shows  a  schematic  defining  off-axis  loading  in  this  case. 
The  angle  between  the  applied  load  g^x  and  the  lamina  fibers  is  defined  as  6,  thus,  the  lamina 
stresses  Gjj,  G22,  and  Xj2  must  be  determined  through  transformation  of  Gxx  to  the  lamina 
coordinate  system.  For  instance,  Gn  =  Gxx  and  G22  =  0  for  0  =  0°.  Likewise,  G22  =  Oxx  and 
O/;  =  0  for  9  =  90°.  Figure  6  shows  the  predictions  of  the  six  criteria  for  a  graphite  epoxy  lamina 
having  the  material  properties  given  in  table  1.  It  is  immediately  evident  all  six  criteria  predict 
very  similar  failure  stress  Gxx  over  the  entire  off-axis  range.  Figure  7  zooms  in  on  the  region 
0  =  0°  -  10°  to  facilitate  the  discussion. 


HGURE  5.  OFF-AXIS  LOADING  OF  A  LAMINA 


HGURE  6.  COMPARISON  OF  LAMINA  FAILURE  CRITERIA  FOR  OFF-AXIS  LOADING 
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HGURE  7.  DETAILED  COMPARISON  OF  LAMINA  FAILURE  CRITERIA  FOR 
OFF-AXIS  LOADING 

All  six  lamina  failure  criteria  predict  a  failure  load  of  X  at  0  =  0°  and  F  at  0  =  90°.  The 
Maximum  Stress  criterion  predicts  three  separate  failure  regions  representing  the  three  possible 
modes  of  failure;  fiber  breakage,  transverse  matrix  cracking,  and  shear  matrix  cracking.  As  the 

off-axis  angle  rotates  from  0°  to  90°,  the  stress  distribution  in  the  lamina  varies.  Over  the  course 
of  fiber  rotation,  failure  is  predicted  by  three  different  equations  thus  producing  the  three  separate 
regions.  Between  0  =  0°  and  2.9°,  fiber  breakage  is  predicted  (mode  1).  At  0  =  0°,  the  failure 
load  is  simply  X.  As  0  increases  in  this  region,  the  predicted  failure  load  actually  increases  only 
slightly,  see  figure  7.  This  is  because  the  Gu/X  ratio  remains  dominant  even  though  the  Gn 
component  in  the  lamina  decreases;  thus,  a  larger  applied  load  is  necessary  to  satisfy  the 
dominant  equation.  At  the  critical  angle  0  =  2.9°,  the  Xjz/S  ratio  becomes  dominant;  therefore, 
the  failure  mode  switches  to  shear  matrix  cracking  (mode  3).  This  region  continues  until  0  = 
27.1°  where  failure  mode  switches  again  to  transverse  matrix  cracking  (mode  2)  with  the  G22IY 
ratio  becoming  dominant.  This  region  continues  through  90°  where  the  failure  load  is  simply  F. 
It  is  important  to  note  that  these  critical  angles  representing  changes  in  failure  mode  are  specific 
to  this  material  system.  Using  a  different  system  would  change  these  critical  angles,  though  for 
most  fiber  reinforced  composites,  similar  transition  angles  should  be  expected. 
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The  Maximum  Strain  criterion  produces  a  failure  curve  similar  to  that  of  Maximum  Stress.  From 
0  =  0°  to  2.9°,  failure  occurs  in  mode  1.  Due  to  Poisson’s  effect,  Maximum  Strain  predicts  a 
slightly  higher  failure  load  than  Maximum  Stress  in  this  region.  In  the  shear  region  (mode  3), 
Maximum  Strain  results  are  identical  to  the  Maximum  Stress  results.  Maximum  Strain’s  shear 
region  extends  to  29°  where  failure  switches  to  mode  2.  Again,  Poisson’s  effect  slightly 
increases  the  Maximum  Strain  prediction  compared  with  Maximum  Stress  until  they  meet  at  0  = 
90°.  Using  measured  ultimate  strains  that  include  nonlinear  effects  would  insignificantly  alter 
the  character  of  this  failure  curve.  Since  the  differences  in  the  failure  predictions  of  Maximum 
Strain  and  Maximum  Stress  are  so  slight,  they  are  plotted  as  one  curve  in  both  figures  6  and  7. 

The  two  interactive  criteria,  Hill-Tsai  and  Tsai-Wu,  give  exactly  the  same  predictions  in  the  three 
failure  regions  as  the  Maximum  Stress  criterion,  with  critical  angles  at  2.9°  and  27.1°  for  mode 
1-3  and  mode  3-2  transitions,  respectively.  Because  these  criteria  are  completely  interactive, 
their  failure  curves  remain  smooth  throughout  the  entire  off-axis  loading  case.  This  is  of  special 
note  in  the  fiber-shear  dominated  region  from  0  =  0°  to  roughly  0=10°  seen  in  figure  7.  In  this 
area,  both  Cn  and  have  significant  contributions.  Because  Z  »  S,  a  criterion  which  couples 
these  stress  components  (e.g.,  Hill-Tsai  and  Tsai-Wu)  will  predict  a  noticeably  lower  value  in 
this  area  than  a  limit  criterion  (Maximum  Stress  and  Strain).  For  the  region  where  both  Xn  and 
022  have  significant  contributions,  the  difference  between  the  limit  and  interactive  criteria 
diminishes  because  Y  and  S  are  of  similar  magnitude.  The  two  interactive  criteria  eventually 
converge  with  the  two  limit  criteria  at  0  =  90°  (only  O22  exists). 

The  two  separate  mode  criteria  exhibit  characteristics  of  both  the  limit  and  interactive  criteria. 
They  yield  the  same  three  failure  regions  as  the  other  criteria.  Due  to  its  ability  to  separate 
modes,  Hashin-Rotem’s  failure  prediction  is  identical  to  Maximum  Stress  in  the  fiber  {On) 
dominated  region  0  =  0°  to  2.9°.  After  the  mode  1-3  transition  at  2.9°,  Hashin-Rotem’s  failure 
prediction  begins  to  move  away  from  Maximum  Stress  and  towards  the  prediction  of  Hill-Tsai 
(See  figure  7).  Note  that  as  0  approaches  90°,  the  Ou  component  becomes  insignificant,  leaving 
only  the  O22  and  X12  components.  Hill-Tsai  and  Hashin-Rotem  are  identical  in  the  O22  -  '^12  plane, 
thus  the  merging.  A  final  failure  load  of  Y  is  predicted  at  0  =  90°  as  expected. 

The  Hashin  predictions  initially  coincide  with  Hill-Tsai’s  predictions  since  they  both  couple  Ou 
and  Xi2-  After  the  mode  1-3  transition,  the  Hashin  criterion  continues  to  predict  the  fiber 
dominant  behavior  which  couples  Ou  and  Xn,  suggesting  shear  failure  instead  of  fiber  failure. 
Since  Hill-Tsai  begins  to  account  for  the  growing  O22  term,  the  Hashin  prediction  becomes 
slightly  higher.  At  about  0  =  10°,  the  predicted  failure  mode  of  the  Hashin  criterion  switches  to 
matrix  failure  and,  with  On  vanishing,  the  prediction  merges  with  Hill-Tsai  at  0  =  90°. 

2.2.3  Pure  Shear. 

The  six  lamina  failure  criteria  are  compared  by  their  failure  predictions  in  a  pure  shear  loading 
situation.  The  angle  between  the  applied  shear  load  and  the  lamina  fibers  is  defined  as  0  as  in  the 
unidirectional  off-axis  loading  previously  discussed.  The  shear  loading  is  given  by  ±Zxy  {Oxx  = 
Oyy  =  0).  A  lamina  composed  of  the  material  in  table  1  is  used  for  this  analysis. 
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Figure  8  shows  a  +Txy  loading  case  for  a  single  lamina  rotated  at  an  angle  6.  All  six  criteria 
predict  a  failure  load  of  5  at  0  =  0°  and  90°  as  expected.  The  criteria  are  all  symmetric  about 
6  =  45°  where  a  failure  load  of  approximately  -Y'  (compressive)  is  predicted.  Due  to  this 
symmetry,  the  discussion  will  go  through  0  =  45°  only.  The  criteria  all  predict  shear  matrix 
cracking  from  0°  to  30.8°  and  transverse  matrix  cracking  from  30.8°  to  45°.  As  the  lamina  fiber 
direction  rotates  from  0°  to  90°,  the  lamina  shear  stress  +X72  goes  down  as  the  other  two 
components  Gji  and  G22  become  larger;  for  positive  shear.  On  becomes  more  tensile  and  O22 
becomes  more  compressive  though  they  have  equal  magnitudes.  For  0  =  0°,  X12  =  Xxy  and 
ct;/  =  C22  =  0.  For  0  =  45°,  T12  =  0  and  On  =  -G22  =  '^xy-  The  stress  components  On  and  G22 
always  remain  of  equal  magnitude,  thus  fiber  failure  never  occurs  (Y' «  X). 


% 


& 


Loading  Angle 


HGURE  8.  COMPARISON  OF  LAMINA  FAILURE  CRITERIA  FOR  POSITIVE 
PURE  SHEAR 

For  the  Maximum  Stress  criterion,  the  mode  3-2  transition  point  dramatically  alters  the  failure 
prediction.  The  G22IY'  equation  becomes  dominant  after  the  transition  at  0  =  30.8°.  Because  G22 
is  smaller  at  30.8°  than  at  45°,  this  equation  actually  requires  a  larger  ultimate  %xy  to  satisfy  the 
dominant  equation.  This  explains  the  drastic  change  in  the  failure  curve.  Maximum  Stress  and 
Maximum  Strain  differ  only  in  this  transverse  failure  region  due  to  Poisson’s  effect  and  are 
plotted  as  one  curve.  They  both  predict  a  maximum  failure  load  at  the  shear-transverse  failure 
transition. 
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The  interactive  criteria  produce  smooth  curves  for  this  loading  case,  even  through  they  do 
differentiate  between  modes  of  failure.  The  linear  terms  in  Tsai-Wu  gives  its  failure  curve  a 
different  character  than  Hill-Tsai,  though  both  clearly  eliminate  the  cusp  formed  by  the  limit 
criteria.  Both  separate  mode  criteria  and  Hashin  and  Hashin-Rotem  match  up  closely  with  Hill- 
Tsai.  The  separate  mode  criteria  deviate  slightly  due  to  the  small  ai;/X  contribution  found  in  the 
Hill-Tsai  criterion  near  0  =45°  where  Cu  becomes  a  maximum.  All  four  interactive  and  separate 
mode  criteria  reach  a  maximum  at  45°  where  <522  is  dominant  and  at  its  maximum  ( =  -Y' ). 

Figure  9  shows  a  -Xxy  loading  case  for  a  single  lamina  rotated  to  an  angle  0.  T  -ike,  the  +Xxy  loading 
case,  all  criteria  are  symmetric  about  0  =  45°  and  switch  from  a  shear  matrix  cracking  failure  to 
transverse  matrix  cracking  at  0  =  13.4°.  At  0  =  0°  and  90°,  a  strength  of  -5  is  predicted  for  all 
criteria  while  a  strength  of  approximately  -Y  (tensile)  is  predicted  at  0  =  45°.  The  applied  load 
Xxy  transforms  into  lamina  stress  components  as  in  the  +Xxy  case,  but  because  the  sign  of  the 
applied  shear  is  opposite,  Ou  is  compressive  and  (J22  is  tensile  after  transformation.  It  is 
interesting  to  note  that  the  mode  3-2  transition  occurs  much  quicker  in  the  -Xxy  case  than  that  of 
+Xxy  because  Y  <  Y'.  This  allows  the  dominate  shear  region  to  switch  over  to  transverse  failure 
sooner. 

The  limit  criteria  predict  an  increasing  failure  load  in  the  shear  failure  region.  This  is  due  to  a 
decrease  in  the  magnitude  of  the  component  Xj2,  even  though  the  X12/S  ratio  is  dominant;  higher 
loads  are  predicted  in  order  to  satisfy  the  dominant  equation.  Maximum  Stress  and  Maximum 
Strain  differ  slightly  in  the  transverse  failure  region  due  to  Poisson’s  effect.  Hill-Tsai,  Tsai-Wu, 
Hashin-Rotem,  and  Hashin,  due  to  their  allowance  of  stress  interaction  between  022  and  X12,  all 
predict  smooth  failure  curves,  eliminating  the  sharp  jumps  seen  in  the  limit  criteria.  Hill-Tsai, 
Hashin-Rotem,  and  Hashin  are  nearly  identical  with  only  a  small  On/X  contribution  in  the  Hill- 
Tsai  criterion  separating  them.  Tsai-Wu’ s  failure  curve  varies  from  the  other  interactive  criteria 
due  to  its  linear  terms,  primarily  the  C22  term.  Again,  because  the  shear  loading  is  negative,  the 
lamina  is  weakest  at  0  =45°,  in  contrast  with  the  positive  shear  case  in  which  it  is  strongest 
at  45°. 

In  contrast  to  off-axis  loading,  in  the  case  of  pure  shear  loading  these  criteria  predict  quite 
different  lamina  strengths. 
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FIGURE  9.  COMPARISON  OF  LAMINA  FAILURE  CRITERIA  FOR  NEGATIVE 
PURE  SHEAR 


2.3  COMPARISON  WITH  EXPERIMENTAL  DATA. 

The  ability  of  the  lamina  failure  criteria  to  correctly  predict  failure  strength  can  be  evaluated  by 
comparing  with  experimental  results.  Among  other  factors,  the  accuracy  of  these  criteria 
depends  on  availability  of  reliable  material  strength  data,  i.e.,  X,  X',  Y,  Y'  and  S,  or  the 
corresponding  ultimate  strains.  Except  for  longitudinal  and  transverse  tensile  strengths  X  and  Y, 
good  measurements  of  the  compressive  and  shear  strengths  are  not  easy  to  obtain,  which  makes 
an  objective  assessment  of  the  lamina  failure  criteria  all  the  more  difficult. 

Assuming  that  reliable  uniaxial  strength  properties  are  available  evaluating  the  failure  criteria, 
failure  loads  of  a  lamina  must  be  determined  for  a  combined  state  of  stress;  i.e.,  at  least  two  of 
the  three  stress  components,  Ojj,  (J22  and  X;2  must  be  present.  The  off-axis  tension  test  offers  the 
simplest  way  to  produce  a  combined  state  of  stress. 

7  3  1  T  amina  Failure  Criteria  Comparison  With  Off-Axis  Tension  Data. 

Many  authors  have  performed  off-axis  unidirectional  lamina  tensile  tests.  This  test  is  performed 
by  loading  a  uniform  coupon  specimen  to  failure.  Because  all  the  lamina  failure  criteria  predict 
very  similar  failure  loads,  correlation  with  experimental  data  can’t  be  used  as  a  means  of  ranking 
these  criteria.  Tests  of  a  boron-epoxy  system  by  Pipes  and  Cole  [8]  illustrate  this  point.  The 
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strength  properties  of  this  material  system  are  provided  in  table  2.  Using  those  values,  the 
theoretical  predictions  of  the  six  lamma  failure  criteria  along  with  the  experimental  data  are 
plotted  in  figure  10.  Only  angles  of  15°,  30°,  45°,  and  60°  were  tested.  For  this  range  of  angles, 
matrix  failure  dominates  the  strength.  It  is  seen  that  the  interactive  and  separate  mode  criteria 
yielded  better  predictions  than  the  limit  criteria. 


FIGURE  10.  COMPARISON  OF  LAMINA  FAILURE  CRITERIA  TO  OFF-A^S  DATA 


TABLE  2.  STRENGTH  VALUES  FOR  THE  BORON-EPOXY  MATERIAL 
SYSTEM  IN  [8]. 


X 

1296.2  MPa 

X' 

-2489.0  MPa 

Y 

62.1  MPa 

Y' 

-310.3  MPa 

S 

68.5  MPa 

Other  sets  of  off-axis  unidirectional  data  yield  similar  conclusions.  Hashin  and  Rotem  [9]  tested 
a  glass-epoxy  system  at  a  number  of  off-axis  angles.  The  Hashin-Rotem  criterion  correlates  with 
the  data  nearly  perfectly,  though  other  criteria  (such  as  Hill-Tsai  and  Hashin)  would  have  been 
just  as  close. 
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2.3.2  Lamina  Failure  Criteria  Comparison  With  Tubular  Specimens. 


In  this  section  the  use  of  uniaxial  (hoop  wound)  tubular  specimens  is  discussed.  Use  of  a  tubular 
specimen  allows  a  biaxial  state  of  stress  to  be  applied  to  a  composite  laminate.  Tubular 
specimens  also  eliminate  the  free  edge  effect  found  in  flat  coupon  specimens,  as  has  been 
verified  by  many  authors  including  Colvin  and  Swanson  [10]. 

Wu  and  Scheublein  [11]  generated  biaxial  lamina  data  using  a  graphite-epoxy  (Morganite  II) 
system  with  material  constants  shown  in  table  3.  Figure  1 1  shows  the  predictions  for  the  material 
system  versus  experimental  data  (Cu  -  C22  plane).  Clearly  all  criteria  match  up  at  the  four  axis 
intercepts  since  those  data  points  were  used  to  generate  the  failure  envelopes.  The  data  in  the 
tensile-compressive  quadrant  are  not  sufficient  to  show  any  distinction  between  the  various 
criteria. 

-  -  -  -  HUl-Tsai 

"  Maximum  Stress,  Hashin-Rotem,  Hashin 

-  -  -  -  Maximum  Strain 

-  —  Tsai-Wu 

□  Glass-Epoxy  Data  from  Wu  and  Scheublein 


a  ^  Ultimate  (MPa) 


HGURE  11.  COMPARISON  OF  LAMINA  FAILURE  CRITERIA  TO  G//  -  G22  DATA 
FROM  WU  AND  SCHEUBLEIN 

TABLE  3.  STRENGTH  VALUES  FOR  THE  GRAPHITE-EPOXY  MATERIAL 
SYSTEM  IN  [11]. 


X 

1027.3  MPa 

X' 

-710.2  MPa 

Y 

43.4  MPa 

Y' 

-125.5  MPa 

S 

72.4  MPa 
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Jiang  and  Tennyson  [12]  used  tubular  specimens  to  characterize  the  failure  of  an  IM7/8551-7 
graphite-epoxy  material  system.  The  elastic  and  strength  constants  for  this  system  are  shown  in 
table  4.  The  tubes  were  tested  in  all  three  biaxial  planes,  <5n  -  G22,  (Sii  - 112,  and  C22  -  '^12- 

TABLE  4.  MODULI  AND  STRENGTH  VALUES  FOR  IM7/855 1-7  GRAPHITE- 
EPOXY  [12]. 


Figures  12,  13,  and  14  show  the  predictions  based  on  the  material  properties  versus  the 
experimental  strengths.  Again,  the  limited  data  inhibits  any  attempt  to  distinguish  the 
performance  among  the  criteria. 

- Hill-Tsai 

. .  Maximum  Stress,  Hashin-Rotem,  Hashin 

-  -  -  -  Maximum  Strain 

—  —  Tsai-Wu 


□  IM7/8551-7  Data  from  Tennyson  and  Jiang 


-2000  -1000  0  1000  2000  3000 
CTjj  Ultimate  (MPa) 


FIGURE  12.  COMPARISON  OF  LAMINA  FAILURE  CRITERIA  TO  On  -  On  DATA 
FROM  JIANG  AND  TENNYSON  [12] 
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Ultimate  (MPa) 


-  -  -  -  Hill-Tsai,  Hashin  (+  ) 

Maximiim  Stress  &  Strain,  Hashin'Rotem,  Hashin  ( *  ) 

—  —  Tsai-Wu 

□  IM7/8551-7  Data  from  Tennyson  and  Jiang 


-1500  -1000  -500  0  500  1000  1500  2000  2500 

Ultimate  (MPa) 


HGURE  13.  COMPARISON  OF  LAMINA  FAILURE  CRITERIA  TO  an  -  On  DATA 
FROM  JIANG  AND  TENNYSON[12] 
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-  -  -  -  Hiil-Tsai,  Ha8hin>Rotem,  Hashin 
— Maximum  Stress  &  Strain 

—  —  Tsai-Wu 

□  1M7/8551>7  Data  from  Tennyson  and  Jiang 
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FIGURE  14.  COMPARISON  OF  LAMINA  FAILURE  CRITERIA  TO  c>22-c>j2  DATA 
FROM  TENNYSON  AND  JIANG 

Swanson  et  al.  [13]  obtained  strength  data  for  an  AS4/55A  unidirectional  composite  in  the 
<^22  -  't/2  plane.  Table  5  lists  the  strength  properties.  Figure  15  contains  theoretical  predictions  of 
the  six  lamina  failure  criteria  compared  with  the  experimental  data.  The  plot  shows  that  Hill- 
Tsai,  Tsai-Wu,  Hashin-Rotem  and  Hashin  (all  including  O22  -  'I12  stress  interaction)  predict  the 
data  very  well  for  tensile  022.  However,  for  the  combination  of  -022  and  X72,  only  Tsai-Wu 
performs  well,  hi  fact,  the  test  data  indicate  that  lamina  shear  strength  increases  as  the  a22 
component  becomes  compressive. 

In  order  to  verify  the  aforementioned  phenomenon,  two  more  independent  sets  of  <322  -  in  biaxial 
data  were  analyzed.  The  first  set  is  a  T800/3900-2  graphite-epoxy  tested  by  Swanson  and  Qian 
[14].  The  second  set  is  from  tests  by  Voloshin  and  Arcan  [15]  using  glass-epoxy  (Scotch-Ply 
Type  1002).  Both  sets  of  material  strength  constants  are  also  given  in  table  5.  Figures  16  and  17 
show  the  predictions  from  the  six  lamina  failure  criteria  versus  experimental  data.  The  trend  of 
an  increasing  shear  strength  as  the  022  term  becomes  more  compressive  is  again  seen.  Further 
discussion  on  this  phenomenon  will  be  given  in  the  next  section. 
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TABLE  5.  STRENGTH  VALUES  FOR  MATERIAL  SYSTEMS  [13-15]  USED 
IN  C22  - 112  BIAXIAL  FAILURE  COMPARISONS 


AS4/55A 

1 

C 

k. 

c 

5cotch-Ply 
rype  1002) 

1 

T800/3900-2 

X 

not  provided 

X 

1108.0  MPa 

X 

T 

X' 

-617.8  MPa 

X' 

Y 

Y 

Y 

Y' 

-94.7  MPa 

Y' 

-137.30  MPa 

Y' 

-200.0  MPa 

S 

51.8  MPa 

S 

36.92  MPa 

S 

100.0  MPa 

-  -  -  -  HUI'Tsai,  Hashin-Rotem,  Haahin 
— —  Maxifflum  Stress  &  Strain 
- Tsai-Wu 


□  AS4/55A  Data  from  Swanson  et  al 


HGURE  15.  COMPARISON  OF  LAMINA  FAILURE  CRITERIA  TO  <522  -  Xn  AS4/55A 
DATA  FROM  SWANSON,  MESSICK,  AND  TIAN 
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Ultimate  (MPa) 


-  -  -  -  Hill’Tsai,  Hashin-Rotem,  Haahin 
—  Maximum  Stress  &  Strain 

- Tsai-Wu 

□  Voloshin  &  Arcan  Data 


HGURE  17.  COMPARISON  OF  LAMINA  FAILURE  CRITERIA  TO  (J22  -  X12 
GLASS-EPOXY  DATA  FROM  VOLOSHIN  AND  ARCAN 


2.4  ANALYSIS  OF  THE  PHYSICAL  BASIS  FOR  LAMINA  FAILURE  CRITERIA. 

All  existing  lamina  failure  criteria  for  composites  are  phenomenological  or  macromechanical  in 
approach.  In  other  words,  they  are  more  or  less  curve-fitting  techniques.  In  order  to  judge  then- 
adequacy  in  failure  prediction,  we  need  to  understand  the  failure  mechanisms  in  the  fiber/matrix 
system. 

It  is  reasonable  to  say  that  there  are  three  failure  modes  in  composites,  namely,  fiber  failure, 
matrix  failure,  and  fiber/matrix  interfacial  failure.  Since,  in  general,  the  stresses  in  the  fiber  and 
matrix  are  different,  their  respective  failures  are  determined  by  different  strengths.  Thus,  a 
conceptually  sound  lamina  failure  criterion  must  start  by  separating  the  stress  states  in  the  fiber 
and  matrix. 
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2.4. 1  Fiber  Failure. 


Assume  that  the  stresses  in  the  fiber  are  C//,  <52!  and  xj.  For  convenience  of  discussion,  we 
take  the  stress  quadratic  form  as  the  failure  criterion: 


J  2 


,0 


22 


J  2 
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11 


CT  f  X  ,J  2 

'^22  \  t  12 


(8) 


where  Xf,Yf,  and  Sf  are  the  fiber  strengths.  Transverse  isotropy  in  the  strength  of  the  fiber  is 
assumed. 

For  unidirectional  fiber  composites  with  fiber  volume  fraction  Cf,  the  fiber  stresses  are 
approximately  related  to  the  composite  stresses  as 


(T  /  —  ^  ^  /y  /  _  /y 

^11  “  ^  '  ^22  “  ^  22»  ^12  “  ^ 


12 
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(9) 


The  longitudinal  composite  strength  X  is  related  to  the  longitudinal  composite  strength  Xf  as 

Z 


(10) 


Using  equations  9  and  10,  equation  8  can  be  expressed  in  the  form 


(11) 


The  values  of  G22  and  Xn  are  limited  by  Y  and  S,  respectively.  Since  Yf»  Y,  Z  »  F,  and 
Sf»  S,  the  fiber  failure  criterion  of  equation  1 1  can  be  approximated  by 


-^11  ^ 


(—)  =1 
X 


(12) 


This  justifies  the  fiber  failure  criterion  (equation  8)  used  in  Maximum  Stress,  Maximum  Strain, 
and  Hashin-Rotem. 


2.4.2  Matrix  Failure. 

Matrix  failure  is  recognized  as  matrix  cracking  along  the  fiber  direction.  If  cracking  occurs  in 
the  matrix,  then  all  three  matrix  stress  components  ,  G22"' ,  and  Xn’"  are  to  be  included  in  the 
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matrix  quadratic  failure  criterion: 
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(13) 


It  is  generally  agreed  that  the  fiber/matrix  interface  is  the  weaker  surface,  and  so-called  matrix 
cracking  may  actually  occur  along  the  fiber/matrix  interface.  Hence,  failure  is  governed  by  the 
interfacial  stresses  C22"  and  'Z12"  .  If  Cu”'  is  neglected  and  matrix  stresses  and  strengths  are 
assumed  as 


'  22 


m  _ 
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equation  13  reduces  to 

(^/+(^/=-(  (15) 

I  O 

From  the  above  discussion,  it  is  obvious  that  Y  and  S  in  equation  15  should  not  be  interpreted  as 
only  the  tensile  and  shear  strengths,  respectively,  of  the  neat  resin. 

Thus,  from  the  consideration  of  failure  mechanisms  in  the  composite  system,  the  criteria  for  fiber 
failure  and  matrix  failure  should  be  separated.  Equations  12  and  15  turn  out  to  be  exactly  those 
proposed  by  Hashin  and  Rotem  in  equation  5. 

2.4.3  Generation  of  Failure  Envelopes  in  the  Stress  Planes. 

In  composites,  tensile  strengths  X  and  Y  are  very  different  from  the  compressive  strengths  X' 
and  Y'.  For  example,  in  using  the  Hill-Tsai  criterion,  proper  strength  values  X  or  X'  and  Y  or  Y' 
must  be  selected  based  on  the  stress  quadrant  to  be  analyzed.  This  has  been  considered 
inadequate,  and  a  single  equation  was  desired  with  the  result  of  the  Tsai-Wu  failure  criterion.  In 
2-D  plane  stress,  Tsai-Wu  is  essentially  Hill-Tsai  with  additional  linear  terms  in  Gn  and  022- 
These  linear  terms  allow  Tsai-Wu  to  account  for  compressive  and  tensile  stresses.  Such  a 
polynomial  cannot  be  related  to  the  concept  of  deformation  energy. 

In  view  of  the  different  tensile  and  compressive  failure  mechanisms,  there  is  no  reason  for  the 
lamina  failure  envelope  to  be  described  by  a  single  equation  as  suggested  by  Tsai-Wu.  It  is 
difficult  to  argue  that,  for  example,  failure  of  a  composite  under  biaxial  tension  should  depend  on 
its  compressive  strength  properties  and  vice  versa.  Although  mathematically  more  convenient, 
such  a  practice  as  adopted  by  Tsai-Wu  may  cause  unreasonable  failure  predictions.  As  shown  in 
figure  2,  Tsai-Wu  suggests  that  a  compressive  stress  G22  would  increase  the  longitudinal  strength 
of  the  composite.  There  are  no  known  mechanistic  reasons  to  support  this.  The  fact  171  >  171 
causes  the  translation  of  the  failure  ellipse  to  the  said  position. 
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2.4.4  On  the  Maximum  Strain  Criterion. 


On  the  <5n  -  Xn  and  G22  -  t/2  planes  (see  figures  3  and  4),  the  Maximum  Stress  and  Maximum 
Strain  criteria  predict  identical  results.  However,  for  biaxial  loading  in  the  <3ii  -  G22  plane  (see 
figure  2),  these  two  criteria  differ  significantly.  The  Maximum  Strain  criterion  predicts  that  for  a 
tensile  longitudinal  stress  Gn,  the  tensile  transverse  stress  G22  would  be  greater  than  Y  in  order  to 
fail  the  composite.  Specifically,  for  Cu  near  X,  the  a22  required  to  cause  failure  is  approaching 
27.  If  the  transverse  strength  of  the  composite  is  controlled  by  the  fiber/matrix  interfacial 
strength,  then  this  is  not  possible. 

It  is  concluded  that  the  Maximum  Strain  criterion  is  not  adequate  for  predicting  the  transverse 
matrix  cracking  failure  mode  where  Cn  is  present. 

2.4.5  Dependence  of  Shear  Strength  on  Compressive  Normal  Stress 

In  all  the  existing  lamina  failure  criteria,  the  lamina  strengths  X,  Y,  and  S  are  assumed  to  be 
constants.  However,  from  the  three  g22  -  X12  biaxial  plots  in  section  2.3.2. 1  (figures  15-17),  there 
is  strong  evidence  that  when  the  composite  is  subjected  to  a  combined  022  -  Xn  loading,  it 
becomes  stronger  when  g22  is  compressive.  More  specifically,  for  given  g22  =  ±go,  the  shear 
stress  X;2  at  failure  corresponding  to  a22  =  -Go  is  appreciably  greater  than  the  shear  stress  Xj2 
corresponding  to  g22  =  +go. 


This  behavior  indicates  that  a  compressive  fiber/matrix  interfacial  normal  stress  (which  is 
proportional  to  G22)  would  create  a  greater  fiber/matrix  interfacial  shear  strength.  To  reflect  this 
behavior,  the  matrix  failure  criterion  of  equation  15  may  be  modified  to 


where 


(—)\( — — 
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=  1 


G22<0 
0  G  22  >0 


(16) 


The  term  p,  plays  a  role  similar  to  friction  coefficients.  Equation  16,  denoted  the  modified  matrix 
criterion,  still  yields  the  expected  values  of  g22  =  Y  atX]2  =  0  and  'T;2  =S  at  G22  =  0. 

In  the  absence  of  Gju  the  linear  stress  terms  in  the  Tsai-Wu  criterion  (see  equation  4 )  produce  a 
failure  envelope  in  the  g22  -  X/2  plane  that  exhibits  a  characteristic  strengthening  effect  similar  to 
that  of  equation  16.  However,  this  effect  would  be  reversed  in  the  Tsai-Wu  criterion  in  cases  in 
which  the  transverse  compressive  strength  171  is  less  than  the  transverse  tensile  strength  171. 
Thus  the  Tsai-Wu  criterion  can  be  made  to  fit  the  data  reasonably  well,  although  with  no  physical 
basis.  In  fact,  simply  increasing  the  transverse  tensile  strength  of  the  composite  system  would 
translate  the  half  ellipse  in  figures  15-17  up  and  away  from  the  data.  The  use  by  the  Tsai-Wu  of 
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the  same  compressive  and  tensile  strengths  criterion  in  all  four  stress  quadrants  is  again 
considered  suspect. 

Figure  18  shows  AS4/55A  data  of  [13]  (from  figure  15)  plotted  with  the  Tsai-Wu  criterion  and 
the  modified  matrix  criterion,  equation  16.  It  is  clear  from  the  comparison  that  the  modified 
matrix  criterion  fits  the  data  as  well  as  Tsai-Wu.  In  this  particular  case,  (i  =  0.6  fits  the  data  well. 
Figures  19  and  20  are  simply  replots  of  figures  16  and  17,  respectively,  with  the  modified  matrix 
criterion  added.  Only  one  p,  value  (which  most  closely  fits  the  data)  is  chosen  for  plotting. 
These  plots  also  show  the  accuracy  of  this  modified  matrix  approach  in  the  C22  -  '^12  plane.  More 
work  in  this  area  must  be  done  to  correlate  |X  with  S  and  Y'. 


—  —  Tsai-Wu 

—  -  -  -  Modified  Matrix  Criterion  w/ p=0.5 

—  -  —  -  ModiHed  Matrix  Criterion  w/ji;=0-6 
"■■■  ModiHed  Matrix  Criterion  w/  p=0.7 

□  AS4/55A  Data  from  Swanson  et  al 


-1^2  Ultimate  (MPa) 


FIGURE  1 8.  COMPARISON  OF  LAMINA  FAILURE  CRITERIA  AND  THE  MODIFIED 
MATRIX  CRITERION  TO  022  -  X72  AS4/55A  DATA  FROM  SWANSON, 
MESSICK,ANDTIAN 
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—  -  »  -  Hill-Tsai,  Haahin-Hotam,  Hashin 

■  Maximum  Stress  &  Strain 

—  —  Tsai-Wu 

—  •  —  *  Modified  Matrix  Criterion  w/  ^^.8 
□  Voloshin  &  Arcan  Data 


Ultimate  (MPa) 

HGURE  20.  COMPARISON  OF  LAMINA  FAILURE  CRITERIA  AND  THE  MODIFIED 
MATRIX  CRITERION  TO  022  -  in  GLASS-EPOXY  DATA  FROM 
VOLOSHIN  AND  ARCAN 

2.4.6  Concluding  Observations  on  Lamina  Failure  Criteria. 

From  the  comparison  with  experimental  data  and  mechanistic  reasoning,  we  conclude  that  a 
Separate  Mode  Failure  criterion  in  the  form 

— ^  =  1  for  fiber  failure 
X 

(—^)  +( - — — )  =1  for  matrix  failure  (17) 

Y  S-iia  22 

is  most  suitable  for  lamina  failure  prediction. 
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A  modified  form  of  Tsai-Wu  failure  criterion. 


— —  =  1  for  fiber  failure 

F2O 22  +  ^22^  22  +  1^5  22  +  ^66^  n  =  ^  for  matrix  failure  (18) 

is  also  adequate  for  most  advanced  composites  (for  which  lY’l  >  lYI)  based  purely  on  its  relative 
success  in  matching  test  data  rather  than  mechanistic  reasoning. 

3.  LAMINATE  FAILURE  ANALYSIS. 

Classical  laminate  strength  analysis  is  based  on  the  assumption  of  a  two-dimensional  stress  field 
in  the  laminate.  Laminate  failure  is  the  eventual  result  of  progressive  failure  processes  taking 
place  in  the  constituent  laminae  under  loading.  Conceptually,  a  ply-by-ply  failure  analysis  should 
yield  the  desired  failure  load  for  the  laminate.  In  reality,  however,  the  failure  mechanisms  in 
laminates  are  a  great  deal  more  complicated  than  those  in  a  unidirectional  composite  under  plane 
stress.  In  addition  to  the  three  intralaminar  failure  modes  (fiber  failure,  matrix  tension,  and 
matrix  shear  failure)  occurring  at  the  lamina  level,  three-dimensional  failure  mechanisms  are 
present  in  the  laminate,  the  most  notable  of  which  include  delamination  and  failure  induced  by 
free  edge  singular  stresses.  Classical  laminate  strength  analysis  is  restricted  to  those  laminates 
whose  failure  is  not  dominated  by  3-D  failure  modes. 

The  effects  of  free  edge  stresses  are  usually  treated  separately  from  classical  laminate  failure 
analysis.  It  is  thus  generally  assumed  that  the  laminate  is  either  free  from  free  edge  stresses  or 
laminate  failure  does  not  initiate  from  the  free  edge.  Some  authors  have  utilized  tubular 
specimens  to  avoid  the  effect  of  free  edge  stresses.  The  use  of  layers  of  film  adhesive  at  the 
interlayers  can  also  toughen  the  interface,  forcing  failure  to  occur  in  in-plane  modes.  This  latter 
approach  is  taken  in  this  study  to  enable  the  use  of  laminate  coupon  specimens  for  testing 
laminate  strength. 

As  lamina  failure  is  progressive  in  nature,  the  progressive  loss  of  lamina  stiffnesses  must  also  be 
accounted  for  in  the  laminate  analysis.  However,  the  local  stress  concentration  effect  due  to 
matrix  cracks  is  usually  neglected  except  for  laminates  with  thick  laminae  such  as  a  [0/90g/0] 

laminate.  This  local  stress  concentration  effect  on  laminate  strength  was  discussed  by  Sun  and 
Jen  [16]. 

3.1  STIFFNESS  REDUCTION. 

Some  of  the  laminate  failure  analysis  methods  consider  a  laminate  capable  of  load  bearing  after 
an  individual  ply  within  the  laminate  has  failed.  These  methods  require  a  procedure  for 
“discounting”  the  failed  ply  and  reducing  the  laminate  stiffness.  Two  methods  for  achieving  this 
were  developed  for  the  present  study,  the  Parallel  Spring  Model  and  the  Incremental  Stiffness 
Reduction  Model. 
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3.1.1  Parallel  Spring  Model. 


Each  lamina  is  modeled  with  a  pair  of  springs  representing  the  fiber  (longitudinal)  and  matrix 
(shear  and  transverse)  deformation  modes.  The  entire  laminate  is  modeled  by  grouping  together 
a  number  of  parallel  lamina  spring  sets  as  shown  in  figure  21.  When  fiber  breakage  occurs,  the 
longitudinal  modulus  is  reduced.  When  matrix  cracking  occurs,  the  shear  and  transverse  moduli 
are  reduced.  The  value  to  which  the  moduli  are  reduced  was  arbitrary  although  it  was  commonly 
set  equal  to  zero. 


i 
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HGURE  21.  SCHEMATIC  OF  THE  PARALLEL  STIFFNESS  MODEL 

This  model  is  also  capable  of  differentiating  between  types  of  matrix  failure  if  desired;  i.e.,  the 
transverse  and  shear  moduli  can  be  reduced  separately  depending  on  the  specific  type  of  matrix 
failure  mode.  The  model  which  reduces  Ei  for  fiber  failure  and  E2  and  Gn  for  either  transverse 
or  shear  matrix  failure  is  denoted  the  PSM.  The  model  which  reduces  Ei  for  fiber  failure,  E2  for 
transverse  matrix  failure,  and  E2  and  G72  for  shear  matrix  failure  is  denoted  the  PSMg.  The  idea 
behind  the  PSMs  is  that  a  transverse  matrix  failure  doesn’t  necessarily  inhibit  the  ability  of  the 
lamina  to  carry  significant  shear  loads.  Creating  these  two  different  reduction  models  has  little 
micromechanical  basis  and  is  done  mainly  for  curve  fitting  purposes. 
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3.1.2  Incremental  Stiffness  Reduction  Model. 


To  avoid  the  sudden  jump  in  strain  at  ply  failure  seen  in  the  Parallel  Spring  Model,  a  model 
resembling  the  bilinear  hardening  rule  in  classical  plasticity  can  be  formulated.  T  aminatP 
stiffness  reduction  is  achieved  similar  to  the  Parallel  Spring  Model.  However,  it  is  assumed  that 
the  reduced  laminate  stiffness  governs  only  the  incremental  load-deformation  relations  beyond 
immediate  ply  failure. 

Both  of  these  stiffness  reduction  models  have  flexibility.  Instead  of  reducing  the  appropriate 
moduli  suddenly  after  a  ply  failure,  a  nonlinear  function  such  as  exponential  function  may  be 
used  to  gradually  reduce  these  values.  This  progressive  softening  approach  may  model  certain 
laminates  better  than  others,  i.e.,  those  laminates  whose  failure  is  dominated  by  matrix  cracking. 

For  most  fiber-dominated  composites,  setting  the  stiffness  constants  directly  to  zero  after  the 
corresponding  mode  of  failure  is  simple  and  unambiguous.  The  use  of  such  reduction  can  be 
justified  by  regarding  the  laminate  analysis  to  be  at  the  in-plane  (x,  y)  location  where  all  ply 
failures  would  occur.  Consider  a  90°  lamina  (within  a  laminate)  containing  a  number  of 
transverse  matrix  cracks,  as  shown  in  figure  22.  The  90°  ply  still  retains  some  stiffness  in  the 
loading  direction  {E2  direction  locally).  However,  the  assumption  is  made  that  ensuing  0°  fiber 
failure  will  occur  at  the  weakest  point.  This  point  is  where  matrix  cracking  has  occurred  in  the 
90°  plies  or  where  locally  E2  =  0.  Thus,  it  is  acceptable  to  reduce  E2  directly  to  zero  after 
transverse  matrix  cracking  for  the  ultimate  strength  analysis. 


nGURE22.  SCHEMATIC  OF  LAMINATE  WITH  MATRIX  CRACKS 

Since  matrix  cracks  are  discrete,  the  portion  of  a  failed  lamina  between  two  cracks  would  still 
contribute  substantially  to  the  laminate  stiffnesses.  It  is  obvious  that  such  drastic  lamina  stiffness 
reduction,  if  assumed  to  be  true  over  the  whole  laminate,  would  overestimate  the  ultimate  strains 
of  the  laminate.  In  fiber-dominated  laminates,  the  effect  of  matrix  cracks  on  the  overall  laminate 
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stiffness  is  usually  very  small.  It  is  reasonable  to  estimate  the  laminate  ultimate  strains  by  using 
the  virgin  laminate  stress-strain  relations  and  the  laminate  failure  stresses  obtained  from  the 
laminate  failure  analysis. 

3.2  LAMINATE  FAILURE  ANALYSIS  METHODS. 

As  with  lamina  failure  analysis,  a  variety  of  laminate  failure  analysis  methods  have  been 
proposed.  Following  is  a  description  of  each  methodology. 

3.2. 1  Plv-Bv-Plv  Discount  Method. 

This  is  a  very  common  method  for  laminate  failure  analysis.  The  laminate  is  treated  as  a 
homogeneous  material  and  is  analyzed  with  a  lamina  failure  criterion.  Laminated  plate  theory  is 
used  to  initially  calculate  stresses  and  strains  in  each  ply.  A  lamina  failure  criterion  is  then  used 
to  determine  the  particular  ply  which  will  fail  first  and  the  mode  of  that  failure.  A  stiffness 
reduction  model  is  used  to  reduce  the  stiffness  of  the  laminate  due  to  that  individual  ply  failure. 
The  laminate  with  reduced  stiffnesses  is  again  analyzed  for  stresses  and  strains.  The  lamina 
failure  criterion  predicts  the  next  ply  failure  and  laminate  stiffness  is  accordingly  reduced  again. 
This  cycle  continues  until  ultimate  laminate  failure  is  reached. 

A  number  of  definitions  have  been  proposed  on  how  to  determine  ultimate  laminate  failure.  One 
common  way  is  to  assume  ultimate  laminate  failure  when  fiber  breakage  occurs  in  any  ply. 
Another  way  is  to  check  if  excessive  strains  occur  (i.e.,  yielding  of  the  laminate  stiffness  matrix). 
Matrix-dominated  laminates  such  as  [±45]s  may  fail  without  fiber  breakage.  Others  have 
suggested  a  “last  ply”  definition  in  which  the  laminate  is  considered  failed  if  every  ply  has  been 
damaged.  For  this  project,  the  laminate  failure  is  defined  as  occurring  when  either  fiber  breakage 
occurs  in  any  ply  or  the  reduced  stiffness  matrix  becomes  singular. 

3.2.2  Sudden  Failure  Method. 

In  highly  fiber  dominated  composite  laminates,  effect  of  the  laminate  stiffness  reduction  due  to 
progressive  matrix  failures  on  the  laminate  ultimate  strength  is  insignificant.  This  suggests  that 
in  such  laminates  the  progressive  stiffness  reduction  seen  in  the  previous  method  may  be 
unnecessary,  and  laminate  failure  may  be  taken  to  coincide  with  the  fiber  failure  of  the  load¬ 
carrying  ply  (the  ply  with  fibers  oriented  closest  to  the  loading  direction).  To  perform  this 
analysis,  a  lamina  failure  criterion  is  chosen  and  the  failure  load  is  determined  by  calculating  the 
load  required  for  fiber  failure  in  the  dominant  lamina.  No  stiffness  reductions  are  included  in  the 
process.  The  laminate  strength  predicted  by  the  sudden  failure  method  is  usually  higher  than  the 
laminate  strength  predicted  by  the  ply-by-ply  discount  method. 

The  strength  analysis  program  listed  in  appendix  B  includes  both  the  ply-by-ply  discount  method 
and  the  sudden  failure  method. 
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3.2.3  Hart-Smith  Criterion: [17-1 91  The  Truncated  Maximum  Strain  F.nvp.Inpp. 

A  laminate  failure  criterion  was  proposed  by  Hart-Smith  [17-19]  based  on  the  assumption  that 
for  a  fiber  dominated  laminate,  failure  can  be  attributed  to  shear  failures  of  the  fibers,  and  that 
laminate  failure  can  be  treated  as  a  projection  of  a  multiaxial  fiber  failure  criterion  onto  laminate 
stress  space.  The  Hart-Smith  criterion  does  not  require  a  ply-by-ply  discount  procedure.  The 
failure  envelope  is  given  in  the  strain  space  that  corresponds  to  an  “extended”  Tresca  (maximum 
shear  stress)  yield  criterion. 

Initially,  the  Hart-Smith  failure  criterion  was  based  on  experimental  results  (see  reference  17  and 
figure  14  of  reference  19)  on  in-plane  shear  failures  in  ±45°  laminates  which  gave  strengths 
about  half  of  what  is  expected  when  the  shear  stress  is  resolved  into  pure  tension  and 
compression  in  the  fiber  directions.  Hart-Smith  attributed  the  low  laminate  strength 
measurements  to  the  presence  of  biaxial  stresses  which  were  believed  to  induce  shear  failures  in 
the  fibers.  Hart-Smith  also  makes  frequent  reference  to  observations  such  as  those  described  in 
reference  17  on  conically  shaped  fracture  surfaces  of  tension  loaded  carbon  fibers  that  were 
interpreted  as  shear  failures.  From  these  observations  Hart-Smith  concluded  that  in  many  cases 
laminate  failure  can  be  reduced  to  shear  failures  in  the  fibers.  Denoting  L  and  T  as  the  in-plane 
longitudinal  and  transverse  directions  with  respect  to  the  fibers  and  N  the  out-of-plane  normal 
direction,  it  was  suggested  that  these  shear  failures  could  be  characterized  as  “L-T”  failures, 
“L-N”  failures,  and  “T-N”  failures  corresponding  to  differences  between  principal  stresses  in  the 
L  and  T  directions,  the  L  and  N  directions,  and  the  T  and  N  directions;  these  stress  differences 
resolve  into  shear  stresses  in  the  three  planes  lying  at  45  degrees  to  the  three  pairs  of  coordinate 
directions. 

Application  to  a  laminate  containing  0°  and  90°  laminae  under  biaxial  loadings  in  the  0°  and 
90°  directions  is  illustrated  in  figure  23.  In  the  figure,  £«  is  the  tensile  failure  strain  of  the  fiber 
and  £o'  is  the  compressive  failure  strain  of  the  0°  lamina  which  may  be  less  than  8o  because  of 
microbuckling  which  precedes  the  fiber  rupture  under  compression.  (A  method  for  modifying 
the  failure  criterion  for  matrix  failure  under  transverse  tensile  stresses  has  also  been  described  in 
references  17-19  and  elsewhere.)  The  strain  coordinates  Ex  and  £y  are  laminate  strains  which  are 
assumed  to  be  the  same  as  the  strains  in  the  fiber.  In  general  the  Hart-Smith  criterion  is  intended 
for  fiber  dominated  lammates  which  contain  more  than  two  reinforcement  directions.  A 
procedure  for  applying  the  criterion  to  more  general  laminates  which  include  ±45°  reinforcement 
as  well  as  to  loadings  which  include  in-plane  shear  is  discussed  in  the  pertinent  references  (see 
reference  17,  for  example). 

The  failure  envelope  (called  truncated  maximum  strain  criterion  by  Hart-Smith)  shown  in  figure 
23  in  terms  of  strains  in  the  0/90  laminate  is  based  on  superposition  of  the  individual  lamina 
failure  envelopes  such  as  the  one  shown  in  figure  24  for  the  0  degree  lamina,  and  is  essentially 
the  same  as  the  Sudden  failure  method  (no  matrix  failure  is  assumed)  in  conjunction  with  the 
Maximum  Strain  criterion.  The  only  difference  between  figure  23  and  the  conventional 
Maximum  Strain  criterion  is  the  45°  cutoff  lines  in  the  2nd  and  4th  quadrants.  The  lamina  failure 
envelope  shown  in  figure  24  represents  the  projection  onto  lamina  strain  space  of  the  fiber  failure 
surface  based  on  the  three  shear  failure  modes.  The  discussion  in  references  17-19  indicates  that 
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nGURE23.  HART-SMITH’S  TRUNCATED  MAXIMUM  STRAIN  FAILURE 
ENVELOPE  FOR  LAMINATES 

the  45  degree  lines  in  the  second  and  fourth  quadrants  in  figure  24  correspond  to  L-T  failures, 
while  those  in  the  first  and  third  quadrants  corresponding  to  the  lines  at  angles  -a  to  the 
horizontal  and  -p  to  the  vertical  represent  T-N  and  L-N  failures,  respectively.  With  the  out-of¬ 
plane  stress  Gn  equal  to  zero  in  typical  laminate  applications,  the  latter  conditions  amount  to 
constant  Gt  and  constant  Gl  cutoffs,  respectively.  (As  indicated  in  the  notation  at  the  upper  left 
of  figure  24,  a  and  P  are  related  to  the  Poisson  ratios  of  the  lamina  as  required  to  produce  these 
constant  stress  conditions.) 

The  translation  from  the  fiber  failure  surface  to  the  lamina  strain  surface  depends  on  the 
assumption  made  by  Hart-Smith  that  the  longitudinal  and  transverse  strains  in  the  lamina  are  the 
same  as  those  in  the  fiber.  For  the  longitudinal  strains  this  is  obviously  a  valid  assumption. 
However,  for  the  transverse  strains  it  is  suspect,  since  except  for  aramid  fibers,  typical 
reinforcements  are  much  stiffer  than  the  usual  polymer  matrices  so  that  the  fiber  strains  will  tend 
to  be  considerably  less  than  the  effective  lamina  strain  except  for  high-volume  fractions  for 
which  the  fibers  are  nearly  in  contact.  This  will  have  the  effect  of  distorting  the  lamina  strain 
envelope  in  figure  24  and  the  corresponding  laminate  strain  envelope  in  figure  23  to  some  extent, 
although  a  precise  assessment  of  the  modification  would  require  a  micromechanics  analysis. 
There  are  certain  issues  that  need  to  be  resolved  in  connection  with  the  Hart-Smith  approach  to 
failure.  The  most  crucial  is  that  of  whether  or  not  there  are  45-degree  cutoffs  of  the  lamina 
failure  envelope  in  the  second  and  fourth  quadrants.  As  indicated  previously,  without  the  cutoffs 
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FIGURE  24.  FAILURE  STRAIN  ENVELOPE  PROPOSED  BY  HART-SMTTH  FOR  FIBERS 

the  Hart-Smith  criterion  reduces  to  the  maximum  Strain  or  Maximum  Stress  criterion.  Data  on 
failure  of  laminates  under  biaxial  loading  which  is  discussed  later  will  throw  light  on  this  issue. 

In  addition,  the  evidence  for  shear  failures  in  fibers  under  normal  stress  loading  as  well  as  the 
implications  of  low  shear  strength  in  ±45-degree  laminates  need  to  be  explored  more  fully.  The 
experimental  results  cited  by  Hart-Smith  for  these  phenomena  are  somewhat  limited  and  have 
not  been  confirmed  to  any  extent  by  other  workers.  Further  work  in  this  area  appears  to  be 
warranted.  In  this  regard  data  on  torsion  tests  of  ±45-degree  filament  wound  tubes  are  available 
[20]  which  may  throw  some  light  on  the  issue  of  low  shear  strength  in  ±45-degree  laminates. 

3.3  LAMINATE  FAILURE  ANALYSIS  UNDER  BIAXIAL  LOADING. 

Given  the  laminate  failure  analysis  methods  discussed,  the  six  lamina  failure  criteria  can  be 
characterized  by  their  ability  to  predict  failure  in  a  laminate  under  biaxial  loading.  Since  data  for 
quasi-isotropic  laminates  under  biaxial  g^x  -  Gyy  loading  are  readily  available,  it  is  appropriate  to 
start  the  evaluation  here. 
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A  quasi-isotropic  [0/±45/90]s  laminate  is  examined  using  the  ply-by-ply  discount  method  in 
conjunction  with  the  Parallel  Spring  Model  (PSM)  for  laminate  stiffness  reduction.  In  this 
analysis,  the  appropriate  lamina  moduli  reduce  directly  to  zero  at  the  individual  ply  failures.  The 
laminate  is  assumed  to  reach  ultimate  failure  when  any  ply  within  the  laminate  fails  by  fiber 
breakage.  Figure  25  shows  the  predictions  of  the  six  lamina  failure  criteria  assuming  the 
laminate  is  made  up  of  the  material  system  found  in  table  1.  This  laminate  is  subjected  to  pure 
<^xx  -  ctyy  biaxial  loading  with  x^y  =  0. 

. Hill-Tsai 

Maximum  Stress  &  Strain,  Hashin-Rotem,  Hashin 
-  -  -  -  Tsai-Wu 


-1500  -1000  -500  0  500  1000  1500 

Ultimate  (MPa) 

XZ 


HGURE  25.  FAILURE  ENVELOPE  OF  [0/±45/90]s  LAMINATE  USING  PLY-BY-PLY 
DISCOUNT  METHOD  AND  PSM  STIFFNESS  REDUCTION 

Before  the  ply-by-ply  discount  method  can  be  properly  assessed  with  different  lamina  failure 
criteria,  it  is  important  to  understand  the  sensitivity  of  this  method  to  the  applied  loads. 
Figure  26  examines  the  Hill-Tsai  envelope  from  figure  25.  Just  the  tensile  Gxx  quadrants  I  and  IV 
are  shown  here  since  they  suitably  illustrate  the  characteristics  to  be  discussed.  Indeed,  this 
failure  envelope  contains  both  linear  and  quadratic  regions  with  significant  jumps  in  predicted 
failure  loads  as  the  biaxial  loading  sweeps  through  the  two  quadrants.  As  the  biaxial  Oxx  -  ^yy 
load  changes,  the  stresses  and  strains  in  each  lamina  change.  This  fact  compounded  with  large 
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differences  in  tensile  and  compressive  strengths  for  the  material  used  produces  different  ply 
failure  combinations  in  each  region  of  the  failure  envelope. 


Hill-Tsai 


a  Ultimate  (MPa) 

XX 


FIGURE  26.  FAILURE  ENVELOPE  OF  [0/±45/90]s  LAMINATE  USING  PLY-BY-PLY 
DISCOUNT  METHOD  AND  HILL-TSAI  LAMINA  FAILURE  CRITERION 

The  first  quadrant  is  completely  symmetrical  about  45°  (a«  =  Cyy),  due  to  the  nature  of  this 
quasi-isotropic  laminate.  Regions  A  and  D  are  quadratic  because  when  the  fiber  dominant  ply 
fails,  there  is  still  coupling  among  all  components  in  that  ply.  For  example,  in  region  D,  the  0° 
ply  still  carries  Gu,  G22,  and  X72  when  it  finally  fails.  In  regions  B  and  C  on  the  other  hand,  every 
ply  in  the  laminate  experiences  matrix  failure  before  the  final  ply  failure.  Therefore,  G22  =  I12  = 
0  for  every  ply  in  the  laminate.  This  will  result  in  a  linear  relationship  since  there  is  no  coupling 
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of  the  stress  components,  only  <3n  exists  in  each  lamina.  Comparison  of  regions  D  and  E  shows 
how  differences  in  intermediate  ply  failures  do  not  necessarily  alter  the  ultimate  failure  of  the 
laminate.  The  second  ply  failure  in  regions  D  and  E  take  place  in  the  ±45  degree  plies  in  both 
mode  2  and  3.  The  Parallel  Spring  Model  used  in  this  analysis  method  does  not  discriminate 
between  types  of  matrix  failure;  thus,  the  laminate  stiffness  is  reduced  similarly  in  both  regions. 
In  region  F,  the  final  failure  in  the  90°  ply  is  linear  because  there  was  previous  failure  in  the  90° 
ply’s  matrix,  therefore  ^22  =  t:j2  =  0  for  that  ply.  Region  G  shows  how  a  laminate  in  certain 
loading  situations  can  be  very  fiber  dominant.  In  region  H,  the  addition  of  an  initial  ±45°  ply 
failure  slightly  changes  the  character  of  the  previous  curve  in  region  G. 

Again,  a  comparison  of  all  the  lamina  failure  criteria  used  with  this  particular  laminate  failure 
analysis  method  is  shown  in  figure  25.  The  Tsai-Wu  lamina  failure  criterion  yields  a 
combination  of  quadratic  and  linear  regions,  which  is  similar  to  the  Hill-Tsai  criterion  for  the 
same  reasons.  The  characteristic  of  the  Tsai-Wu  failure  envelope  for  this  laminate  is  different 
from  the  Hill-Tsai  prediction  because  of  the  linear  terms  contained  in  Tsai-Wu,  though  similar 
trends  are  seen  between  the  two  criteria.  The  limit  criteria  (Maximum  Stress  and  Maximum 
Strain)  produce  a  nearly  linear  laminate  failure  envelope  as  expected.  The  Maximum  Stress 
criterion  contains  a  small  jump  in  predicted  failure  load  near  the  tensile  uniaxial  loading  cases. 
This  occurs  because  Maximum  Stress  considers  the  ±45°  ply  as  the  critical  ply  of  the  laminate  in 
this  particular  loading  situation.  Ultimate  failure  is  predicted  when  the  ±45°  ply  has  matrix 
failure.  Unrealistic  laminate  strains  occur  if  the  analysis  attempts  to  continue  by  looking  for  a 
fiber  failure. 

Since  Maximum  Strain,  Hashin,  and  Hashin-Rotem  criteria  predict  virtually  the  same  failure 
envelope  as  Maximum  Stress,  they  are  plotted  as  one  curve.  Maximum  Strain  produces  a  slightly 
larger  failure  envelope  due  to  Poisson’s  effect.  The  Hashin  and  Hashin-Rotem  criteria  match  up 
almost  exactly  with  the  limit  criteria  in  the  analysis  of  this  particular  laminate.  Because  Hashin 
and  Hashin-Rotem  allow  for  stress  interaction  between  022  and  X;2,  the  jumps  in  failure  seen  in 
Maximum  Stress  and  Maximum  Strain  do  not  occur;  though  not  seen  in  the  figure,  the  curve 
would  maintain  the  same  slope  through  this  critical  region. 

It  is  clearly  shown  in  figure  25  that  except  in  the  third  stress  quadrant,  all  failure  criteria  yield 
roughly  the  same  laminate  strengths.  We  note  that  in  the  first  quadrant  all  plies  have  failed  in 
matrix  mode  before  the  laminate  fails  by  fiber  breakage.  Although  the  ply  failure  sequences 
predicted  by  the  failure  criteria  could  be  different,  the  end  results  are  the  same.  On  the  other 
hand,  in  the  third  stress  quadrant,  there  are  no  ply  matrix  failures  preceding  fiber  failure,  and  all 
the  stress  components  are  fully  “active”  in  the  failure  criteria,  resulting  in  very  different  laminate 
strength  predictions  between  the  fully  interactive  criteria  (Hill-Tsai  and  Tsai-Wu)  and  other 
criteria. 

Figure  27  shows  the  comparison  between  the  sudden  failure  method  and  the  ply-by-ply  discount 
method  for  the  quasi-isotropic  laminate  under  biaxial  loading.  The  Maximum  Stress  criterion  is 
used  in  the  analysis.  As  expected,  the  sudden  failure  method  predicts  higher  laminate  strength 
except  in  quadrant  HI  where  ply-by-ply  discount  laminate  failure  matches  the  sudden  failure 
assumptions. 
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Maximum  Strain  with  Parallel  Spring  Model 
Maximum  Strain  with  Sudden  Failure  Model 


a  Ultimate  (MPa) 

XX 


FIGURE  27.  FAILURE  STRESS  ENVELOPE  FOR  AN  AS4-3501-6  [0/±45/90]s 
LAMINATE  UNDER  BIAXIAL  LOADS  USING  PLY-BY-PLY 
DISCOUNT  (THE  PSM)  AND  SUDDEN  FAILURE  METHODS 

Figure  28  shows  the  failure  strain  envelope  corresponding  to  the  failure  stresses  of  figure  27.  In 
the  strain  plane,  the  predictions  of  the  ply-by-ply  discount  method  and  the  sudden  failure  method 
agree  everywhere  except  for  a  few  locations  where  the  ply-by-ply  discount  method  exhibits  some 
kinks.  A  careful  examination  of  the  computer  output  reveals  a  special  failure  process  as  follows. 
For  the  loading  condition  in  these  regions  (kinks),  the  0°  (or  90°)  ply  suffers  matrix  failure  first. 
When  the  load  reaches  some  critical  value,  the  +45°  plies  fail  in  matrix  mode,  followed 
immediately  by  the  fiber  failure  of  the  0°  (or  90°)  ply  without  any  increase  of  loading.  In  other 
words,  the  matrix  failure  of  ±45°  plies  leads  to  catastrophic  failure  of  the  laminate.  Thus,  the 
ultimate  failure  strain  of  the  laminate  for  the  loading  ratios  in  that  range  is  the  failure  strain 
corresponding  to  the  matrbc  failure  of  +45°  plies,  which  is  smaller  than  the  fiber  failure  strain. 
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MMrimum  Strain  with  Parallel  Spring  Model 


Maximum  Strain  with  Sudden  Failure  Model 
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HGURE  28.  FAILURE  STRAIN  ENVELOPES  CORRESPONDING  TO  THE  FAILURE 
STRESS  ENVELOPE  OF  FIGURE  27 


3.3.1  Comparison  of  Data  for  Biaxial  Loading. 

Recently  Swanson  et  al.  [19,  21]  performed  biaxial  testing  on  an  AS4/3501-6  graphite-epoxy 
[0/±45/90]s  laminate  using  tubular  specimens.  The  ply  properties  given  by  Swanson  are  listed  in 
table  6.  Note  that  they  are  slightly  different  from  those  given  by  table  1.  Also,  the  longitudinal 
modulus  Ej  is  taken  as  the  average  of  the  initial  modulus  and  the  secant  modulus  reported 
in  [21]. 

TABLE  6.  MODULI  AND  STRENGTH  VALUES  OF  AS4/3501-6  GRAPHITE-EPOXY 
SYSTEM  IN  [19, 21]. 


El 

134.60  GPa 

■ 

X 

1986.0  MPa 

E2 

11.03  GPa 

■ 

X' 

-1193.0  MPa 

Gj2 

5.52  GPa 

■ 

Y 

47.9  MPa 

^12 

0.28 

■ 

Y' 

-168.0  MPa 

Ply  Thickness: 

0.13  mm 

■ 

S 

95.7  MPa 

Figure  29  shows  Swanson’s  data  in  the  first  and  fourth  quadrants  of  the  Cxx  -  Oyy  plane  plotted 
against  various  failure  envelopes.  The  predicted  envelopes  are  generated  with  the  six  lamina 
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failure  criteria  in  conjunction  with  the  ply-by-ply  discount  method.  Comparing  with  the  data,  it 
is  evident  that  the  predictions  of  all  six  lamina  failure  criteria  agree  with  the  data  very  well. 
However,  as  discussed  in  section  3.3,  the  fully  interactive  criteria  such  as  Hill-Tsai  and  Tsai-Wu 
are  very  sensitive  to  the  sequence  of  lamina  matrix  failures  and  can  be  significantly  affected  by 
the  sudden  reduction  of  matrix  stiffness,  resulting  in  jumps  in  strength,  as  shown  in  figure  29. 


- Hill-Tiai 

Mixiinnin  Stren  &  Strain,  TT«jihiw.RntflTn,  Hashin 
—  —  Taai-Wu 


-1000  -750  -500  -250  0  250  500  750  1000 
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FIGURE  29.  COMPARISON  OF  ULTIMATE  STRESS  ENVELOPES  WITH 

EXPERIMENTAL  DATA  FOR  A  [0/±45/90]s  LAMINATE  UNDER 
BIAXIAL  LOADS 

Swanson  and  Qian  [14]  also  tested  laminate  tubes  in  the  c^x  -  Cyy  plane.  The  [0/±45/90]s, 
[03/±45/90]s,  and  [0/(±45)2/90]s  data  presented  in  reference  14  (I  and  IV  quadrants  only)  matched 
well  with  the  Maximum  Strain  analysis.  Considering  that  all  lamina  failure  criteria  are  in  close 
agreement  in  these  quadrants  for  the  ti/4  laminate,  it  is  clear  that  they  should  all  be  adequate  for 
these  laminates. 

3.3.2  Biaxial  Failure  in  the  Strain  Plane. 

The  failure  strains  corresponding  to  the  failure  stress  envelopes  of  figure  29  are  plotted  in 
figure  30.  It  is  interesting  to  note  that  the  failure  strain  envelopes  predicted  by  Maximum  Stress, 
Maximum  Strain,  Hashin,  and  Hashin-Rotem  criteria  essentially  coincide  with  the  limits  set  by 
the  ultimate  tensile  strain  (Xg)  and  compressive  strain  (XgO  of  the  unidirectional  composite. 
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However,  these  limit  strains  should  not  be  automatically  taken  as  the  ultimate  strains  of  the 
laminate  at  failure.  For  the  7t/4  quasi-isotropic  laminate  under  biaxial  loading  in  the  third  stress 
quadrant  (compressive  Gxx  and  Gyy),  there  are  no  ply  matrix  failures  before  ultimate  laminate 
failure.  Thus,  there  is  no  stiffness  reduction  before  laminate  failure,  and  the  calculated  failure 
strain  is  the  laminate  failure  strain. 


- mn-Tud 

Maadmum  Strew  &  Strain,  Hwhin-Botem,  TTaAin 
“  —  Twi-Wu 

□  Experimental  Data  from  Swanam  and  Traak  [21] 
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FIGURE  30.  COMPARISON  OF  ULTIMATE  STRAIN  ENVELOPES  WITH 

EXPERIMENTAL  DATA  FOR  A  [0/+45/90]s  LAMINATE  UNDER 
BIAXIAL  LOADS 


Under  tensile  biaxial  loading  (tensile  Gxx  and  Gyy),  all  plies  in  the  laminate  have  suffered  matrix 
failure  before  laminate  final  failure.  The  drastic  reduction  (to  zero)  of  ply  transverse  and  shear 
stiffnesses  tends  to  overestimate  the  ultimate  strains.  Therefore,  the  strain  failure  envelope  as 
shown  in  figure  30  should  not  be  used  directly  in  laminate  strength  design  without  accounting  for 
the  stiffness  reductions. 

In  figure  31,  the  Hart-Smith  criterion  is  compared  with  the  Maximum  Strain  criterion  (with  the 
ply-by-ply  discount  method)  and  the  experimental  data  of  Swanson  and  Trask  [21]. 
Hart-Smith  [17]  used  different  strength  data;  otherwise  his  prediction  would  have  coincided  with 
the  Maximum  Strain  criterion  except  for  the  45°  shear  cutoffs.  From  the  data,  the  45°  cutoff 
does  not  seem  to  be  necessary. 
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■ .  Maximcm  Strain 

— —  Hart-Smith  [17] 

□  Experimental  Data  from  Swanaon  and  Traak  [21] 
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FIGURE  3 1 .  COMPARISON  OF  MAXIMUM  STRAIN  AND  FLART-SMITH’S  ULTIMATE 
STRAIN  ENVELOPES  WITH  EXPERIMENTAL  DATA  FOR  A  [0/±45/90]s 
LAMINATE  UNDER  BIAXIAL  LOADS 

3.3.3  Biaxial  Testing  Data  For  Glass  Woven  Fabric  Composite. 

Recently,  Wang  and  Socie  [23-24]  performed  biaxial  testing  on  NEMA  G-10  E-glass  plain 
woven  fabric  composite  laminates.  They  used  both  tubular  specimens  and  flat  square  specimens. 
Strength  data  were  obtained  for  all  four  biaxial  loading  quadrants.  Their  test  results  clearly 
indicate  that  the  laminate  failure  strain  envelope  is  bounded  by  the  uniaxial  ultimate  strains  of  the 
laminate.  Specifically,  they  found  that  failure  in  one  direction  was  not  affected  by  loading  in  the 
transverse  direction.  Hence,  the  failure  strain  envelope  appears  to  be  rectangular  as  predicted  by 
the  Maximum  Strain  and  Maximum  Stress  criteria.  For  this  composite  laminate,  the  45°  shear 
cutoff  as  proposed  by  Hart-Smith  was  not  observed. 

3.4  LAMINATE  STRENGTH  ANALYSIS  FOR  UNIDIRECTIONAL  OFF-AXTS  LOADING. 

Laminate  strength  data  are  available  for  coupon  specimens  under  uniaxial  loading.  In  most  cases 
though,  free  edge  stresses  control  the  initiation  and  final  failure  of  these  laminates.  Data  as  such 
are  not  suitable  for  evaluating  the  laminate  failure  analysis  methods  as  attempted  here.  However, 
by  placing  film  adhesive  at  the  interfaces  of  laminate  coupon  specimens,  it  is  possible  to  suppress 
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these  3-D  effects.  This  type  of  coupon  specimen  is  used  to  generate  additional  laminate  strength 
data  for  the  purpose  of  evaluating  the  failure  criteria. 

3.4. 1  Generation  of  Laminate  Failure  Curve  for  Off-Axis  Loading. 

Theoretically  predicting  laminate  strength  under  off-axis  loading  is  similar  to  the  lamina  case 
(see  figure  5).  All  plies  of  the  laminate  are  simply  rotated  the  same  amount  as  the  angle  of 
loading.  The  newly  created  laminate  is  then  subjected  to  a  unidirectional  load  in  the  global  x- 
direction.  For  example,  a  [0/±45]s  laminate  under  10°  off-axis  loading  would  be  analyzed  as  a 
[10/-i-55/-35]s  laminate  under  a  unidirectional  load  in  the  x-direction.  The  failure  curve  for  a 
laminate  under  such  loading  is  generated  by  calculating  the  ultimate  strength  at  each  off-axis 
loading  angle  in  a  desired  range.  Typical  laminate  symmetry  requires  only  a  portion  of  the  entire 
360°  range  to  be  analyzed. 

3.4.2  Selection  of  Laminates  and  Off-Axis  Lxjading  Angles. 

A  variety  of  laminates  was  theoretically  examined  in  advance  of  actual  strength  testing.  Only  the 
ply-by-ply  discount  method  was  used  to  determine  ultimate  laminate  strength.  Laminate  layups 
and  loading  angles  were  chosen  so  as  to  provide  a  comparison  of  the  six  lamina  failure  criteria. 
The  focus  of  the  experiment  was  to  examine  whether  these  failure  criteria  could  predict  the 
correct  trend  of  the  laminate  strength  versus  the  loading  angle.  It  was  not  necessary  to  evaluate 
the  ability  of  the  lamina  failure  criteria  to  predict  laminate  strength  over  the  entire  possible  range 
of  loading  angles. 

Table  7  shows  the  laminates  and  off-axis  loading  angles  selected.  Included  in  these  tests  were 
unidirectional  laminate  specimens  for  principal  material.  Five  or  more  tests  were  performed  at 
each  off-axis  loading  angle  to  provide  accurate  mean  results.  In  the  table,  “A”  indicates  the 
location  of  an  adhesive  film  in  the  laminate. 

TABLE  7.  LAMINATES  AND  OFF-AXIS  LOADING  ANGLES  TESTED 


Laminate 

Off-Axis  Angles  Tested 

[0]8 

0°  &  90° 

[0/A/+45/A/-45/A/90]s 

0°  -  22.5°  every  7.5° 

[90/A/0/A/90/A/0]s 

0°  -  7.5°  every  1.5°,  15°,  22.5° 

[0/A/-i45/A/-45]s 

0°  -  30°  every  7.5°,  26°,  45° 

[90/A/+30/A/-30]s 

0°  -  22.5°  every  7.5° 

3.4.3  Consideration  of  Curing  Stresses  and  in  Situ  Lamina  Strength. 

The  mismatch  of  thermal  expansion  coefficients  of  the  fiber  and  matrix  can  cause  residual 
thermal  stresses  (curing  stresses)  to  form  within  the  laminate  during  the  manufacturing  process. 
Typically  these  stresses  place  the  matrix  in  tension.  These  stresses  can  be  significant;  thus,  they 
should  not  be  neglected  if  prediction  of  matrix  failure  is  of  importance. 
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Another  issue  in  laminate  analysis  is  the  in  situ  strength  of  constituent  laminae.  Flaggs  and  Kural 
[25]  found  that  the  in  situ  transverse  strength  Y  could  be  as  high  as  2.5  times  the  unidirectional 
transverse  strength.  This  can  clearly  be  a  cause  of  inaccurate  strength  predictions. 

As  it  turns  out,  curing  stresses  in  graphite/epoxy  composite  laminates  are  typically  of  a 
magnitude  similar  to  the  unidirectional  transverse  strength  Y.  Since  the  in  situ  lamina  strength  is 
not  available,  in  this  study,  the  effect  of  curing  stresses  is  neglected.  The  rationale  is  that  the 
underestimation  of  in  situ  strength  Y  will  be  offset  by  neglecting  the  effect  of  curing  stresses. 

3.4.4  Laminate  Coupon  Specimens. 

The  material  used  for  testing  was  AS4/3501-6  graphite-epoxy  from  Hercules.  A  0.13-mm-thick 
ply  of  film  adhesive  was  added  at  each  lamma  interface  except  the  middle  interface  because  of 
symmetry.  The  adhesive  was  FM  1000  marketed  by  AmericanCyanamid.  The  elastic  and 
strength  properties  of  this  material  were  determined  by  Sun  and  Zhou  [7]  and  are  listed  in  table  8. 

TABLE  8.  MATERIAL  PROPERTIES  FOR  FM  1000  FILM  ADHESIVE 


^ult 

38  MPa 

E 

1.724  GPa 

G 

0.648  GPa 

V 

0.33 

0.127  mm 

Guidelines  given  by  ASTM  [26]  were  followed  in  the  preparation  of  the  coupon  specimens. 
Specimens  were  cut,  using  a  water  jet,  from  12-  by  12-inch  square  laminate  panels  manufactured 
in  the  manner  recommended  by  Hercules.  All  specimens  were  0.75  inch  wide  with  a  6.75  inch 
gage  section.  This  yielded  a  9  to  1  aspect  ratio.  Water-jet  cutting  gave  smooth  edges  which  were 
free  from  initial  matrix  cracks  and  delamination.  X-ray  photography  was  used  to  verify  the  free 
edge  condition  for  a  few  individual  specimens. 

End  tabs  were  used  to  protect  them  from  being  damaged  in  the  test  machine  grips.  A  standard 
epoxy  adhesive  was  used  for  this  purpose.  For  those  laminate  specimens  with  anisotropic 
stiffness,  it  was  important  not  to  overlook  the  effect  of  shear-extension  coupling.  Highly 
anisotropic  laminate  coupon  specimens  can  fail  prematurely  due  to  stress  concentrations  in  the 
tab  region  if  rectangular  tabs  are  used.  In  order  to  accommodate  the  deformation  induced  by  the 
extension-shear  coupling,  oblique  end  tabs  suggested  by  Sun  and  Chung  [27]  were  used.  See 
figure  32.  The  oblique  angle  ^  was  derived  using  extensional  laminate  stiffnesses  Ai;  and  as 

cot  ^  =  -  Aic/Aii  (27) 
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end  tab 


FIGURE  32.  EXAMPLE  OF  OBLIQUE  END  TABS 

It  was  demonstrated  by  Sun  and  Chung  [27]  that,  using  oblique  end  tabs,  an  almost  uniform  state 
of  stress  corresponding  to  uniaxial  loading  could  be  generated  in  off-axis  laminate  coupon 
specimens  rigidly  gripped  during  loading.  In  quasi-isotropic  laminates,  extension-shear  coupling 
is  absent  and  standard  rectangular  tabs  were  used. 

Glass-epoxy  and  graphite-epoxy  quasi-isotropic  laminates  were  used  for  oblique  end  tabs.  For 
the  rectangular  tabs  ,  material  stiffness  is  not  as  critical,  and  chopped-fiberglass  circuit  board  was 
used.  The  tab  length  was  determined  by  assuming  an  ultimate  failure  load  and  then  sizing  the 
necessary  tab  area  based  on  the  strength  of  the  epoxy  system  used  for  tab  bonding.  In  all  cases, 
the  tabs  measured  0.75  inch  wide  and  a  minimum  of  1.25  inches  long.  Some  oblique  tabs  were 
extended  to  2.25  inches  in  order  to  assure  good  adhesion  to  the  coupon  specimen. 
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3.4.5  Testing  Procedure. 

Tests  were  performed  in  the  Composite  Materials  Laboratory  (CML)  at  Purdue  University.  All 
specimens  were  mechanically  loaded  in  tension  at  a  stroke  rate  of  0.1  inch/min  on  an  MTS  810 
servohydraulic  test  machine.  This  was  equal  to  a  strain  rate  of  approximately  1.5%  per  minute, 
which  adheres  to  ASTM  guidelines  for  fiber  composites.  Tests  were  performed  at  room 
temperature.  The  majority  of  tests  were  designed  to  determine  ultimate  load  and  in  some  cases 
ultimate  strain.  Hence,  only  load,  displacement,  and  strain  (if  applicable)  were  measured  during 
tensile  loading.  These  values  were  displayed  in  real  time  through  a  data  acquisition  package  on  a 
personal  computer. 

The  first  set  of  coupons  tested  were  0°  and  90°  unidirectional  laminates.  All  specimens  were 
strain  gauged  to  determine  the  appropriate  elastic  constants.  Table  9  lists  the  results.  Only  the 
moduli  £■/,  E2,  Vn  and  the  strength  values  X  and  7  were  determined  in  the  present  test.  The 
strengths  S,  X',  and  Y'  were  taken  from  a  similar  AS4/3501-6  material  system  given  by  Sun  and 
Zhou  [7].  The  shear  modulus  Gn  in  [7]  was  obtained  from  testing  [±45]s  laminate  and  is  lower 
than  6.9  GPa  as  obtained  by  Daniel  [28]  using  a  more  reliable  method.  This  value  is  listed  in 
table  9. 

TABLE  9.  MODULI  AND  STRENGTH  VALUES  FOR  THE  TESTED  AS4/3501  -6 
GRAPHITE-EPOXY  SYSTEM 


El 

153.7  GPa 

X 

2171.0  MPa 

E2 

11.0  GPa 

X 

-2013.0  MPa 

Gi2 

6.9  GPa 

Y 

67.0  MPa 

V12 

0.32 

Y 

-206.8  MPa 

Ply  Thickness: 

0.13  mm 

S 

110.3  MPa 

Ultimate  laminate  stress  for  the  unidirectional  laminate  coupon  specimens  (0°  and  90°)  was 
determined  by  dividing  the  measured  ultimate  load  by  the  cross-sectional  area  (width  by 
thickness).  The  8-ply  unidirectional  laminates  had  an  average  thickness  of  1.04  mm;  thus  the 
nominal  ply  thickness  was  assumed  as  0.13  mm  (0.0051  in.). 

For  laminates  containing  layers  of  film  adhesive,  the  calculation  of  ultimate  stress  must  account 
for  the  thickness  and  stiffness  of  the  film  adhesive  in  determining  the  true  graphite-epoxy 
laminate  strength.  First,  it  is  assumed  that  the  total  measured  ultimate  load  is  the  summation  of 
the  load  carried  by  the  composite  and  the  load  carried  by  the  adhesive,  i.e., 

Pexp=Pc+PA  (28) 


It  is  also  assumed  that  the  composite  and  adhesive  loads  can  be  separated  by  the  rule  of  mixtures: 


Pa  = 


Pa^a 


exp 


(29) 
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f 


(30) 


where  Ea  is  the  adhesive  modulus,  Ha  is  total  thickness  of  all  adhesive  layers  in  the  lanainate,  and 
An  is  the  total  extensional  stiffness  calculated  just  for  the  graphite-epoxy  plies.  Equation  30  is 
used  to  obtain  the  true  composite  ultimate  load  from  the  experimental  load.  The  laminate 
strength  is  then  determined  by  dividing  Pc  by  the  total  cross-sectional  area  of  the  composite 
plies.  Table  10  shows  all  the  averaged  ultimate  stress  data  for  the  tested  laminates. 

TABLE  10.  ULTIMATE  LAMINATE  STRESSES  (MPa) 


Off-Axis  Loading  Angle 

0° 

7.5° 

15° 

22.5° 

26° 

o 

O 

CO 

45° 

- 

r0/A/+45/A/-45/A/90]s 

765 

752 

774 

832 

- 

- 

- 

- 

[0/A/+45/A/-45]s 

843 

929 

mm 

1074 

818 

- 

[90/A/+30/A/-30]s 

908 

837 

807 

- 

- 

- 

Off-Axis  Loading  Angle 

0° 

1.5° 

3° 

4.5° 

6° 

7.5° 

15° 

22.5° 

[90/A/0/A/90/A/0]s 

1126 

1140 

1074 

1018 

861 

394 

288 

Though  the  majority  of  ultimate  failures  were  caused  by  fiber  rupture,  a  few  specimens  were 
matrix  dominated  in  which  fiber  failure  did  not  occur.  In  the  matrix  dominated  coupon  tests,  the 
specimen  never  fractured  into  two  (or  more)  pieces.  Instead,  the  deformation  continued  until 
reaching  a  strain  of  approximately  3%  at  which  the  test  was  stopped.  As  seen  in  figure  33,  a 
[90/A/0/A/90/A/0]s  laminate  loaded  at  22.5°  exemplifies  this  behavior.  For  comparison,  the 
stress-strain  curve  for  a  [0/A/-i-45/A/-45]s  laminate  loaded  at  45°  is  also  plotted,  which 
demonstrates  the  behavior  involving  typical  fiber  failure. 
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HGURE  33.  STRESS-STRAIN  CURVES  CHARACTERIZING  FIBER  AND  MATRIX 
FAILURES 


3.4.6  Comparison  With  Test  Data. 

Theoretical  predictions  of  the  different  lamina  failure  criteria  using  the  ply-by-ply  discount 
method  for  laminates  under  off-axis  loading  are  compared  with  the  experimental  data.  Both  the 
PSM  and  PSMs  are  considered  to  assist  in  analyzing  experimental  trends. 

Figures  34  and  35  show  theoretical  predictions  for  the  [0/A/+45/A/-45/A/90Js  laminate  obtained 
using  the  PSM  and  PSMs  stiffness  reduction  procedures,  respectively.  The  theoretical 
predictions  with  the  PSMs  are  slightly  higher  in  regions  where  transverse  matrix  cracking  has 
occurred.  This  is  because  the  transversely  failed  lamina  retain  shear  stiffness  which  results  in 
higher  lamina  stiffness  and,  consequently,  laminate  strength. 

The  experimental  data  (CML  data)  for  the  [0/A/+45/A/-45/A/90]s  laminate  shows  an  increase  in 
strength  as  the  off-axis  loading  angle  rotates  from  0°  to  22.5°.  Maximum  Stress,  Maximum 
Strain,  and  Hashin-Rotem  all  theoretically  predict  this  increase,  supporting  the  idea  of  separating 
fiber  failure  (governed  by  fiber  stress  Oif)  from  matrix  failure  (governed  by  matrix  stresses  G22 
and  X12). 
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HGURE  34.  COMPARISON  OF  ULTIMATE  STRENGTHS  FOR  A  [0/A/445/A/-45/A/90]s 
LAMINATE  UNDER  UNIDIRECTIONAL  LOADING  USING  DIFFERENT 
LAMINA  FAILURE  CRITERIA  (WITH  THE  PSM) 
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HGURE  35.  COMPARISON  OF  ULTIMATE  STRENGTHS  FOR  A  [0/A/+45/A/-45/A/90]s 
LAMINATE  UNDER  UNIDIRECTIONAL  LOADING  USING  DIFFERENT 
LAMINA  FAILURE  CRITERIA  (WITH  THE  PSMs) 
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Figures  36  and  37  compare  CML  data  with  the  theoretical  predictions  for  the  [90/A/0/A/90/A/0]s 
laminate,  using  the  PSM  and  PSMs  procedures,  respectively.  The  data  favor  the  use  of  the  PSMs, 
suggesting  that  this  laminate  at  small  off-axis  angles  retains  its  shear  stiffness  after  transverse 
matrix  cracking  in  the  90°  plies.  The  fully  interactive  criteria  of  Hill-Tsai  and  Tsai-Wu 
underestimate  the  strength  of  the  laminate.  The  Hashin  criterion  is  closer  to  the  data  than  the 
fully  interactive  criteria.  Again,  the  criteria  (Maximum  Stress,  Maximum  Strain,  and  Hashin- 
Rotem)  which  separate  fiber  failure  from  matrix  failure  best  match  the  data. 


Loading  Angle 


HGURE  36.  COMPARISON  OF  ULTIMATE  STRENGTHS  FOR  A  [90/A/0/A/90/A/0]s 
LAMINATE  UNDER  DIRECTIONAL  LOADING  USING  DIFFERENT 
LAMINA  FAILURE  CRITERIA  (WITH  THE  PSM) 
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HGURE  37.  COMPARISON  OF  ULTIMATE  STRENGTHS  FOR  A  [90/A/0/A/90/A/0]s 
LAMINATE  UNDER  UNIDIRECTIONAL  LOADING  USING  DIFFERENT 
LAMINA  FAILURE  CRITERIA  (WITH  THE  PSMs) 

Data  from  the  [0/A/+45/A/-45]s  laminate  are  displayed  in  figures  38  and  39.  The  difference 
between  the  PSM  and  the  PSMs  is  slight  due  to  the  fact  that  the  laminate  is  dominated  by  shear 
matrix  failure  (for  which  both  reduction  models  treat  the  same).  Overall,  Maximum  Stress, 
Maximum  Strain,  and  Hashin-Rotem  are  clearly  the  best  fit  for  the  data.  In  the  off-axis  region 
from  0°  to  20°,  the  Tsai-Wu  criterion  closely  fits  the  data.  However,  the  overall  trend  of  Tsai- 
Wu’s  predication  is  quite  different  from  the  experimental  data. 
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HGURE  38.  COMPARISON  OF  ULTIMATE  STRENGTHS  FOR  A  [0/A/+45/A/-45]s 

LAMINATE  UNDER  UNIDIRECTIONAL  LOADING  USING  DIFFERENT 
LAMINA  FAILURE  CRITERIA  (WITH  THE  PSM) 


Loading  Angle 


FIGURE  39.  COMPARISON  OF  ULTIMATE  STRENGTHS  FOR  A  [0/A/+45/A/45-45]s 
LAMINATE  UNDER  UNIDIRECTIONAL  LOADING  USING  DIFFERENT 
LAMINA  FAILURE  CRUERIA  (WITH  THE  PSMs) 
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Data  from  the  [90/A/+30/A/-30]s  laminate  are  displayed  in  figures  40  and  41.  As  with  the 
[0/A/+45/A/-45]s  laminate,  differences  between  the  PSM  and  the  PSMg  are  slight.  Once  again. 
Maximum  Stress,  Maximum  Strain,  and  Hashin-Rotem  match  both  the  magnitude  and  the  trend 
of  the  data. 


0  5  10  15  20  25  30  35  40  45  50  55  60 

Loading  Angle 


nOURE  40.  COMPARISON  OF  ULTIMATE  STRENGTHS  FOR  A  [90/A/+30/A/-30]s 
LAMINATE  UNDER  UNIDIRECTIONAL  LOADING  USING  DIFFERENT 
LAMINA  FAILURE  CRITERIA.  (WITH  THE  PSM) 
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HGURE  41.  COMPARISON  OF  ULTIMATE  STRENGTHS  FOR  A  [90/A/+30/A/-30]s 
LAMINATE  UNDER  UNIDIRECTIONAL  LOADING  USING 
DIFFERENT  LAMINA  FAILURE  CRITERIA  (WITH  THE  PSMs) 

3.5  OBSERVATIONS  ON  LAMINATE  FAILURE  CRITERIA. 

From  the  four  different  sets  of  laminate  strength  data  presented  in  the  previous  section,  it  is  clear 
that  the  interactive  criteria  (Hill-Tsai  and  Tsai-Wu)  and  Hashin  criterion  (which  couples  a;;  and 
X/2)  not  only  underestimate  ultimate  laminate  failure  but  cannot  correctly  predict  the  trend  of  the 
data.  On  the  other  hand.  Maximum  Stress,  Maximum  Strain,  and  Hashin-Rotem  criteria  all 
perform  quite  well. 

Simply  put,  those  criteria  which  separate  fiber  failure  completely  from  matrix  failure  are 
relatively  insensitive  to  inaccurate  lamina  strengths  Y  and  S.  There  is  some  matrix  strength 
sensitivity  in  these  criteria  from  the  effects  of  intermediate  ply  failures  (failure  preceding  ultimate 
fiber  failure).  However,  the  fully  interactive  criteria  and  the  Hashin  criterion  are  considerably 
more  sensitive  to  inaccurate  matrix  strengths.  This  can  be  illustrated  by  increasing  the  lamina 
matrix  strengths  used  in  the  theoretical  analysis  of  the  tested  laminates. 

The  sensitivity  of  the  interactive  criteria  to  matrix  strength  is  demonstrated  by  the  ultimate 
strength  curves  of  the  [0/A/+45/Ay-45/A/90]s  laminate  in  figure  42.  The  Hill-Tsai  failure 
criterion  is  plotted  using  three  different  lamina  shear  strength  values.  The  PSM  is  used  for  this 
analysis.  Thus,  the  curve  labeled  “Hill-Tsai  with  S  =  1.0”  is  the  same  as  the  Hill-Tsai  curve  from 
figure  34.  The  other  two  curves  in  figure  42  are  for  the  strengths  2  x  5  and  2.5  x  S,  respectively. 
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FIGURE  42.  COMPARISON  OF  ULTIMATE  STRENGTHS  FOR  A  [0/A/+45/A/-45/A/90]s 
WITH  DIFFERENT  LAMINA  SHEAR  STRENGTHS  USING  HILL-TSAI 
CRITERION 

All  other  material  constants  remain  the  same.  It  is  clear  that  by  simply  increasing  the  shear 
strength  of  the  lamina,  the  Hill-Tsai  criterion  takes  on  a  completely  different  trend,  predicting  an 
increase  in  strength  instead  of  a  decrease  as  the  off-axis  loading  angle  increases  from  0°  to  22.5°. 
With  these  high  shear  strengths,  this  criterion  now  correctly  predicts  the  trend  and  magnitude  of 
the  experimental  data. 

Another  issue  often  raised  is  whether  matrix  cracking  actually  occurs  in  composite  laminates 
with  well  dispersed  laminae.  To  answer  this  question,  we  tested  [0/90/0]  and  [0/90/0]s  laminates 
of  AS4/3501-6  graphite/epoxy  composite.  Coupon  specimens  were  tested  under  tension. 
Figure  43  clearly  shows  the  presence  of  matrix  cracks  in  the  90°  plies  at  the  load  about  95%  of 
the  laminate  strength.  Thus,  we  conclude  that  matrix  cracking  does  occur  and  the  ply-by-ply 
discount  process  in  laminate  failure  analysis  is  justified. 
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FIGURE  43.  EVIDENCE  OF  90°  PLY  MATRIX  CRACKS  IN  [0/90/0]  AND  [0/90/0]s 
LAMINATES  OF  AS4/3501-6  GRAPHITE-EPOXY  COMPOSITE 

4.  CONCLUSIONS. 

In  this  study,  the  following  conclusions  were  obtained. 

•  -  At  the  lamina  level,  those  criteria  (such  as  the  Hashin-Rotem  criterion)  which  separate  the 

fiber  failure  mode  from  the  matrix  failure  mode  are  the  most  reasonable  and  accurate. 
This  is  supported  by  test  data  and  a  micromechanical  consideration  which  indicates  that 
fiber  failure  and  matrix  failure  should  be  governed  by  different  failure  criteria. 

•  For  fiber-dominated  laminates,  Maximum  Stress,  Maximum  Strain,  and  Hashin-Rotem 
failure  criteria  outperform  other  criteria.  These  criteria  are  insensitive  to  variations  in 
matrix  strengths  (7  and  S)  which  are  very  difficult  to  obtain  in  situ. 
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•  The  interactive  failure  criteria  (Hill  Tsai,  Tsai-Wu,  and  Hashin)  are  sensitive  to  variation 
of  the  matrix-dominated  lamina  strengths  (i.e.,  Y  and  5).  Accurate  in  situ  composite 
strengths  are  critical  to  the  use  of  these  criteria.  Because  of  the  interaction  among  all 
stress  components,  sudden  switching  of  failure  modes  makes  the  failure  envelope  (in 
stress  or  strain)  very  jumpy. 

•  Experimental  results  indicate  that  matrix  cracking  does  take  place  even  in  laminates  with 
well  dispersed  laminae.  Thus,  the  ply-by-ply  discount  of  stiffnesses  in  failed  laminae  is 
justified. 

•  The  Parallel  Spring  model  for  stiffness  reductions  is  adequate  for  analysis  of  laminate 
strengths.  The  drastic  ply  stiffness  reduction  (the  concerned  stiffness  is  set  equal  to  zero 
after  ply  failure)  does  not  cause  appreciable  errors  in  the  predicted  laminate  strength  for 
fiber-dominated  laminates. 

•  At  this  point,  available  failure  data  does  not  appear  to  support  the  45-degree  fiber  shear 
failure  cutoffs  of  the  Hart-Smith  lamina  failure  criterion,  although  further  effort  on  this 
issue  appears  warranted.  In  addition,  data  supporting  Hart-Smith’s  contention  of  shear 
failures  in  tension  loaded  fibers  together  with  low  shear  strength  values  of  +45-degree 
laminates  is  somewhat  limited  and  needs  to  be  further  developed.  Failure  data  on  ±45- 
degree  filament  wound  tubes  tested  in  torsion  such  as  that  in  reference  20  is  relevant  to 
the  issue  of  low  shear  strength  in  ±45-degree  laminates  and  should  be  examined  further. 

5.  RECOMMENDATIONS. 

•  For  fiber-dominated  laminates,  maximum  stress,  maximum  strain,  and  Hashin-Rotem 
failure  criteria  should  be  used  to  obtain  most  reliable  laminate  strength  prediction. 

•  To  predict  lamina  matrix  failure  in  a  laminate,  the  in  situ  transverse  strength  (Y)  and 
shear  strength  (S)  should  be  used. 

•  For  matrix-dominated  laminates  and  loadings,  all  six  failure  criteria  can  be  used. 
Laminate  failure  should  be  declared  when  excessive  strains  occur. 

•  None  of  the  laminate  failure  criteria  investigated  here  based  on  the  assumption  of  a  2-D 
stress  state  are  valid  for  strength  prediction  at  free  edges  in  a  composite  structure. 
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APPENDIX  A— A  LIST  OF  FAILURE  CRITERIA  [1,2] 
The  Limit  Criteria: 
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,4  1  11  t,/ 

+  CfllCJ22(^  ^  o)"'"  o 


X^5  X  Y  S 


+ 


12 


G11G22 


2  _1 

VX  Y 


G  j  ]  ^22 

IT  ~ 


(G„G22-i:i2^)[cfn<722(^+ll)  +  ^C?ii^+ll«722^  -p(G„+G22)]- 
(G„^  +G22^  +<7iiG22  +1^12^)  =  1 


Malmeister: 

Py^ij  Pyid^ij^kl  ^ijklinn^ij^kl^iiin  ^ 

Tsai-Wu: 

F.a.+F.jG^aj  =1 

Fj^ij  +  F2G22  '^^11^11  "^^22^22  ^^12^11^22  ^66^^12  ”  ^ 

^1  X  ^  Y  Y^  ^  XX^  ^  ^22  "YY^  ^  ^ 

Huang-Kirmser: 

(FiGi)“  +  (F^g^gJ  +  (FykCfiGjG,  J  =  1 
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experimentally 

determined 


Hashin-Rotem: 
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APPENDIX  B— A  COMPUTER  CODE  FOR  STRENGTH  ANALYSIS  OF 

LAMINATED  COMPOSITES 


S  A  L  C  vl.l 


Strength  Analysis  of  Laminated  Composites 


User’s  Guide 


Written  in  FORTRAN 


Overview 


The  purpose  of  this  program  is  to  provide  a  thorough  analysis  of  the  failure  progression  leading 
to  ultimate  failure  in  laminated  composites.  The  program  utilizes  2-D  classical  laminated  plate 
theory  with  a  Ply-by-Ply  Discount  laminate  analysis  method.  Ultimate  laminate  stress  or  strain 
can  be  calculated  for  any  fiber  reinforced  material  system  under  any  combination  of  applied 
stresses  a^,  and 
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THEORETICAL  DISCUSSION 


1.  LAMINA  FAILURE  ANALYSIS. 

The  purpose  of  the  lamina  failure  criterion  is  to  determine  the  strength  and  mode  of  failure  of  a 
unidirectional  composite  or  lamina  in  a  state  of  combined  stress.  All  the  existing  lamina  failure 
criteria  are  basically  phenomenological  in  which  detailed  failure  processes  are  not  described 
(macromechanical).  Further,  they  are  all  based  on  linear  elastic  analysis.  The  failure  load 
(lamina  strength)  is  determined  by  evaluating  the  set  of  equations  provided  in  each  criterion. 
These  equations  compare  lamina  stresses  (or  strains)  to  material  strengths  in  the  principal 
directions. 

The  majority  of  the  lamina  failure  criteria  assume  a  2D  orthotropic  material  which  is  transversely 
isotropic.  Lamina  strengths  are  defined  in  a  composite  system  as  follows  (values  in  stress): 

X&X' :  tensile  and  compressive  strengths  in  fiber  direction. 

Y&Y':  tensile  and  compressive  strengths  in  transverse  direction  (perpendicular  to 

fibers). 

S :  shear  strength 

For  a  strain  based  analysis,  similar  strain  strengths  may  be  used:  X^ ,  Z/,  Y^ ,  Y^',  and  .  The 
ability  of  the  lamina  failure  criterion  to  determine  mode  of  failure  is  essential  in  bringing  this 
analysis  tool  to  the  laminate  level  (an  individual  lamina  failure  within  a  laminate  doesn’t 
necessarily  constitute  ultimate  failure).  Modes  of  failure  are  defined  as: 

Fiber  Breakage  (mode  1):  longitudinal  stress  (cii,)  or  longitudinal  strain  (e,,) 

dominates  lamina  failure. 

Transverse  Matrix  Cracking  (mode  2):  transverse  stress  (a,,)  or  transverse  strain  (£22) 

dominates  lamina  failure. 

Shear  Matrix  Cracking  (mode  3):  shear  stress  (X12)  or  shear  strain  (7,2)  dominates 

lamina  failure. 

It  is  important  to  point  out  that  both  mode  2  and  mode  3  are  matrix  failures.  The  two  modes  are 
separated  because  different  stress  components  cause  the  failure,  though  the  result  is  the  same. 
Hence,  mode  2  should  be  interpreted  as  matrix  cracking  due  to  G72,  and  mode  3  should  be 
interpreted  as  matrix  cracking  due  to  Xj,.  The  notation  is  for  convenience,  and  it  should  be 
assumed  that  the  resulting  matrix  crack  is  the  same  regardless  of  “failure  mode.” 
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1  ■  1  A  SURVEY  OF  LAMINA  FAILURE  CRITERIA. 


Lamina  failure  criteria  can  be  categorized  into  three  groups; 

•  T.imit  Criteria:  These  criteria  predict  failure  load  and  mode  by  comparing  lamina 
stresses  a,, ,  022 ,  and  Xj,  (or  strains  £„  ,  £22 ,  and  y,2 )  separately.  Interaction  between  the  stresses 
(or  strains)  is  not  considered. 

•  Interactive  Criteria:  These  criteria  predict  the  failure  load  by  using  a  single  quadratic  or 
higher  order  polynomial  equation  involving  all  stress  (or  strain)  components.  Failure  is  assumed 
when  the  equation  is  satisfied.  The  mode  of  failure  is  determined  indirectly  by  comparing  the 
stress/strength  ratios. 

•  Separate  Mode  Criteria:  These  criteria  separate  the  matrix  failure  criterion  from  the 
fiber  failure  criterion.  The  equations  can  be  dependent  on  either  one  or  more  stress  components; 
therefore,  stress  interaction  varies  from  criterion  to  criterion  within  this  group.  If  the  satisfied 
equation  contains  only  one  stress  component,  then  the  failure  mode  corresponds  to  that  particular 
direction;  otherwise,  the  failure  mode  can  be  determined  as  is  done  with  the  interactive  criteria  by 
comparing  stress/strength  ratios  of  the  satisfied  equation. 

•  Six  commonly  used  lamina  failure  criteria  are  available  in  the  program  SALC. 

I.imit  Criteria 


•  Maximum  Stress: 


•  Maximum  Strain; 


^.1 

fiber  failure 

X 

II 

transverse  matrix  cracking 

(1) 

^=1 

shear  matrix  cracking 

S 


—  =  1  fiber  failure 

X. 

—  =  1  transverse  matrix  cracking  (2) 

'W  T  ^ 


=  1  shear  matrix  cracking 
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Interactive  Criteria 


•  Hill-Tsai: 

•  Tsai-Wu: 


j  [  +  22  11  ^22^  22”  1 1^  22  ^6^^  12*  “  ^ 

2  2  2 

^1^11  ^2^22  ^11*^11  '*’  ^22*^  22  2Fj2CTjj022  "I"  F(i6’r|2  =  1 


(3) 

(4) 


where 

p  =  —  +  —  p  =  —  +  —  p  =  - L  F  =r  - F  =  —  F  _  experimentally 

'  ^  ’  2  Y  Y' ’  ’*  XX^  '  YY^  '  *’*  ^  ^  determined 


Separate  Mode  Criteria 

•  Hashin-Rotem: 

^=1 

X 

fiber  failure 

f— 

Y  S 

matrix  failure 

(5) 

•  Hashin: 

X  s 

fiber  failure  (tension) 

^.1 

X' 

fiber  failure  (compression) 

(— )'  +  (— )"  =  1 

Y  S 

matrix  failure 

(6) 

For  Maximum  Stress,  Maximum  Strain,  Hill-Tsai,  and  Hashin-Rotem,  the  criterion  is 
generalized  for  either  tensile  or  compressive  stresses:  the  corresponding  (tensile  or  compressive) 
strength  constant  must  be  chosen  based  on  the  sign  of  the  applied  stress.  The  Tsai-Wu  criterion 
is  designed  for  use  in  all  stress  quadrants  of  any  stress  plane,  thus  may  be  used  directly  without 
modification  for  different  stress  signs.  Tsai-Wu’s  criterion  requires  a  biaxial  test  to 
experimentally  determine  the  interaction  term  Fj,.  It  has  been  suggested  to  use  =  1/(2ZX'), 
which  reduces  Tsai-Wu  down  to  the  Hoffman  criterion.  Some  researchers  have  found  the  term  to 
be  insignificant,  and  suggest  setting  it  equal  to  zero. 

There  is  strong  evidence  that  when  a  unidirectional  composite  is  subjected  to  a  combined  -  x,, 
loading,  it  becomes  stronger  when  is  compressive.  More  specifically,  for  given  022  =  ±07,,,  the 
shear  stress  x,2  at  failure  corresponding  to  022  =  -G^  is  appreciably  greater  than  the  shear  stress  x,2 
corresponding  to  a22  =  +^o-  This  behavior  indicates  that  a  compressive  fiber/matrix  interfacial 
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normal  stress  (which  is  proportional  to  O22)  would  create  a  greater  fiber/matrix  interfacial  shear 
strength.  To  reflect  this  behavior,  the  matrix  failure  criterion  of  equation  5  may  be  modified  to 


■f=l 


(7) 


where 


O22<0 

I  0  0  22  >  0 


The  term  p,  plays  a  role  similar  to  friction  coefficients. 

2.  LAMINATE  FAILURE  ANALYSIS  METHODS. 

As  with  lamina  failure  analysis,  a  variety  of  laminate  failure  analysis  methods  have  been 
proposed.  Following  is  a  description  of  each  methodology: 

Plv-bv-Plv  Discount  Method 

This  is  a  very  common  method  for  laminate  failure  analysis.  The  laminate  is  treated  as  a 
homogeneous  material  and  is  analyzed  with  a  lamina  failure  criterion  at  a  mechanistic  level. 
Laminated  plate  theory  is  used  to  initially  calculate  the  state  of  stress  and  strain  in  each  ply  given 
the  global  loading  situation  and  the  material’s  elastic  and  strength  properties.  A  lamina  failure 
criterion  is  then  used  to  determine  the  particular  ply  which  will  fail  first  and  the  mode  of  that 
failure.  A  stiffness  reduction  model  is  used  to  reduce  the  stiffness  of  the  laminate,  due  to  that 
individual  ply  failure.  The  laminate  with  reduced  stiffnesses  is  again  analyzed  for  stresses  and 
strains.  Lamina  failure  criterion  predicts  the  next  ply  failure  and  laminate  stiffness  is  accordingly 
reduced  again.  This  cycle  continues  until  ultimate  laminate  failure  is  reached. 

There  has  been  a  number  of  definitions  proposed  on  how  to  determine  ultimate  laminate  failure. 
One  common  way  is  to  assume  ultimate  laminate  failure  when  fiber  breakage  occurs  in  any  ply. 
Another  way  is  to  check  if  excessive  strains  occur  (i.e.,  a  singular  laminate  stiffness  matrix). 
Matrix-dominated  laminates  such  as  [±45]s  may  fail  without  fiber  breakage.  Others  have 
suggested  a  “last  ply”  definition  in  which  the  laminate  is  considered  failed  if  every  ply  has  been 
damaged.  For  this  program,  the  laminate  is  loaded  until  fiber  breakage  occurs  in  any  ply,  unless 
the  reduced  stiffness  matrix  is  singular  which  denotes  a  matrix  dominated  ultimate  failure. 

Sudden  Failure  Method 

In  highly  fiber-dominated  composite  laminates,  the  laminate  stiffness  reduction  due  to 
progressive  matrix  failures  insignificantly  affects  the  laminate  ultimate  strength.  This  suggests 
that  in  such  laminates  the  progressive  stiffness  reduction  seen  in  the  previous  method  may  be 
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unnecessary,  and  laminate  failure  may  be  taken  to  coincide  with  the  fiber  failure  of  the  load 
carrying  ply  (the  ply  with  fibers  oriented  closest  to  the  loading  direction).  To  perform  this 
analysis,  a  lamina  failure  criterion  is  chosen  and  the  failure  load  is  determined  by  calculating  the 
load  required  for  fiber  failure  in  the  dominant  lamina.  No  stiffness  reductions  are  included  in  the 
process.  The  laminate  strength  predicted  by  the  Sudden  Failure  method  is  usually  higher  than  by 
the  Ply-by-Ply  Discount  method. 

Direct  Laminate  Method 

This  method  examines  the  laminate  as  a  whole,  using  effective  strength  values  in  the  laminate 
principal  directions.  An  equation  or  set  of  equations  using  these  values  and  the  applied  stresses 
predict  the  failure  of  the  laminate  in  a  similar  fashion  to  the  lamina  failure  criteria,  but  at  a 
laminate  level.  Of  course,  the  laminate  strength  values  are  applicable  to  just  the  particular 
laminate  being  analyzed,  hence  a  change  in  layup  could  require  a  different  set  of  laminate 
strengths.  There  seems  little  reason  to  justify  such  an  inflexible  phenomenological  approach  to 
laminate  analysis.  This  method  is  not  offered  in  SALC. 

2.1  STIFFNESS  REDUCTION 


After  an  individual  ply  within  the  laminate  has  failed,  the  following  two  methods  offer  a  way  to 
“discount”  the  failed  ply  and  reduce  the  laminate  stiffness  accordingly: 

Parallel  Spring  Model 

Each  lamina  is  modeled  with  a  pair  of  springs  representing  the  fiber  (longitudinal)  and  matrix 
(shear  and  transverse)  deformation  modes.  The  entire  laminate  is  modeled  by  grouping  together 
a  number  of  parallel  lamina  spring  sets  as  shown  in  Figure  B-1.  When  fiber  breakage  occurs,  the 
longitudinal  modulus  is  reduced.  When  matrix  cracking  occurs,  the  shear  and  transverse  moduli 
are  reduced.  The  value  to  which  the  moduli  are  reduced  is  arbitrary. 

This  model  is  also  capable  of  differentiating  between  types  of  matrix  failure  if  desired;  i.e.,  the 
transverse  and  shear  moduli  can  be  reduced  separately  depending  on  the  specific  type  of  matrix 
failure  mode.  The  model  which  reduces  for  fiber  failure  and  and  G,,  for  either  transverse 
or  shear  matrix  failure  is  denoted  the  PSM.  The  model  which  reduces  £,  for  fiber  failure,  E,  for 
transverse  matrix  failure,  and  E,  and  G,,  for  shear  matrix  failure  is  denoted  the  PSMs-  The  idea 
behind  the  PSMs  is  that  a  transverse  matrix  failure  doesn’t  necessarily  inhibit  the  ability  of  the 
lamina  to  carry  loading  in  the  shear  direction.  Creating  these  two  different  reduction  models  has 
little  micromechanical  basis,  and  is  done  mainly  for  curve  fitting  purposes. 
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Incremental  Stiffness  Reduction  Model 

To  avoid  the  sudden  jump  in  strain  at  ply  failure  seen  in  the  Parallel  Spring  Model,  a  model 
resembling  the  bilinear  hardening  rule  in  classical  plasticity  can  be  formulated.  Laminate 
stiffness  reduction  is  achieved  similar  to  the  Parallel  Spring  Model.  However,  it  is  assumed  that 
the  reduced  laminate  stiffness  governs  only  the  incremental  load-deformation  relations  beyond 
immediate  ply  failure.  This  model  in  not  available  in  SALC. 


For  most  fiber-dominated  composites,  setting  the  stiffness  constants  directly  to  zero  after  the 
corresponding  mode  of  failure  is  simple  and  unambiguous.  The  use  of  such  reduction  can  be 
justified  by  regarding  the  laminate  analysis  to  be  at  the  in-plane  (x,  y)  location  where  all  ply 
failures  would  occur.  Consider  a  90°  lamina  (within  a  laminate)  containing  a  number  of 
transverse  matrix  cracks,  as  shown  in  Figure  B-2.  The  90°  ply  still  retains  some  stiffness  in  the 
loading  direction  (£,  direction  locally).  However,  the  assumption  is  made  that  ensuing  0°  fiber 
failure  will  occur  at  the  weakest  point.  This  point  is  where  matrix  cracking  has  occurred  in  the 
90°  plies,  or  where  locally  =  0.  Thus,  it  is  acceptable  to  reduce  directly  to  zero  after 
transverse  matrix  cracking  for  an  ultimate  strength  analysis.  This  is  the  approach  used  in  SALC. 
Since  matrix  cracks  are  discrete,  between  two  cracks  a  failed  lamina  would  still  contribute  fully 
to  the  laminate  stiffnesses.  It  is  obvious  that  such  drastic  lamina  stiffness  reduction,  if  assumed 
to  be  true  over  the  whole  laminate,  would  greatly  overestimate  the  ultimate  strains  of  the 
laminate.  In  fiber-dominated  laminates,  the  effect  of  matrix  cracks  on  the  overall  laminate 
stiffnesses  is  usually  very  small.  It  is  reasonable  to  estimate  the  laminate  ultimate  strains  by 
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FIGURE  B-2.  SCHEMATIC  OF  LAMINATE  WITH  MATRIX  CRACKS 

using  the  virgin  laminate  stress-strain  relations  and  the  laminate  failure  stresses  obtained  from 
the  laminate  failure  analysis. 

Using  the  Program 

The  program  is  designed  to  calculate  ultimate  laminate  stress  and  strain  for  any  type  of  laminate. 
under  inplane  applied  stresses  and  Thermal  stresses  arising  from  the  manufacturing 

process  may  also  be  calculated  (but  not  included  in  the  failure  analysis).  Following  below  is  a 
step  by  step  explanation  of  how  to  run  the  program.  Input,  output,  or  file  names  associated  with 
the  program  are  identified  with  the  font  courier. 

1.  Input 

These  instructions  are  provided  for  a  unix  operating  system,  though  the  program  should  run 
similarly  on  any  platform.  The  executable  (program)  requires  these  additional  text  files  in  order 
to  run  (examples  at  end  of  guide): 

mat .  info  Text  file  containing  the  number  and  names  of  other  text  files  listing  material 

properties  to  be  used  by  the  code.  First  line  is  the  number  of  text  files  with 
the  names  of  the  material  files  following. 

lam.  info  Text  file  containing  the  number  and  type  of  laminates  to  be  used  by 

program.  First  line  is  the  number  of  laminates  followed  by  a  blank  line. 
Then  each  laminate  is  entered,  first  with  the  number  of  plies,  followed  by 
the  name  of  the  laminate,  followed  by  each  ply  angle  in  order  as  found  in  the 
laminate. 


?????  Any  file  listed  in  mat .  info  which  is  a  material  property  file.  See  example 

at  end  of  guide  for  order  of  material  constants. 

The  user  can  customize  these  files.  The  idea  for  this  setup  is  to  allow  each  user  to  generate  a  list 
of  material  properties  and  laminate  families  frequently  used. 

Once  the  program  is  executed,  it  reads  in  mat .  info  and  provides  a  list  of  material  files  to 
choose  from: 

Enter  material  data  file  for  analysis  (#)  ... 

Enter  in  the  file  #  desired.  Next  the  program  asks  what  type  of  laminate  to  analyze: 


Enter  #  of  laminate/ply  choice... 


Likewise,  enter  in  the  desired  laminate.  Next  the  program  asks  if  thermal  stresses  (strains) 
should  be  calculated.  NOTE,  these  are  not  included  in  the  failure  analysis.  Only  the  mechanical 
strains  are  considered  in  the  failure  analysis: 

Do  you  want  to  calculate  thermal  stresses? 
l:yes  other: no 

If  yes,  the  program  then  prompts, 

Enter  delta  T: 


Enter  the  effective  thermal  drop  during  the  curing  cycle.  Next  the  program  asks  what  type  of 
loading  to  place  the  laminate  under: 

Enter  Mechanical  Loading  Envelope  #: 

1 :  Nx/Ny  2 :  Nx/Nxy  3 :  Ny/Nxy  4 :  Pure  Thermal 

Enter  in  the  desired  biaxial  loading  envelope  in  options  #1-3.  Or  the  user  may  opt  to  just 
examine  the  thermal  stresses  and  strains  by  choosing  option  #4  (no  mechanical  loading  nor 
failure  analysis  is  performed).  Option  #1  assumes  Nxy  =  0.  Option  #2  assumes  Ny  =  0.  Option 
#3  assumes  Nx  =  0.  Once  a  biaxial  loading  plane  is  chosen,  the  program  prompts 

Enter  load-angle  range  and  step,  (i.e.,  0,90,10) 

Figure  B-3.  shows  a  schematic  of  how  to  apply  a  biaxial  load.  Consider  option  #1  (Nx  and  Ny 
applied  with  Nxy  =  0).  The  load  is  applied  as  a  vector  pointing  in  any  circular  direction.  By 
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changing  the  angle  (p)  of  the  vector,  the  components  Nx  and  Ny  are  altered.  These  components 
may  be  either  in  tension  or  compression.  Examples: 


•  pure  unidirectional  tensile  Nx:  enter 

•  pure  unidirectional  compressive  Ny:  enter 

•  complete  Nx-Ny  biaxial  failure  envelope:  enter 


0,0,1  (zero  to  zero  degrees,  one  step) 
270,270,1 

0,360,5  (failure  point  every  5°) 


In  the  last  example,  the  program  calculates  an  ultimate  laminate  failure  at  72  points  (360°/5°) 
around  the  Nx-Ny  biaxial  envelope.  The  user  may  choose  a  finer  or  coarser  number  of  steps,  or 
may  choose  just  to  examine  any  subsection  of  the  entire  envelope.  The  other  two  biaxial  planes 
would  operate  in  a  similar  fashion.  For  all  three  biaxial  cases,  the  first  component  listed  is  on  the 
x-axis  (0°  and  180°)  and  the  second  component  listed  is  on  the  y  axis  (90°  and  270°). 


Ny  (tension) 


Nx  (compression) 


Ny  (compression) 

HGURE  B-3.  SCHEMATIC  OF  APPLIED  LOADING  VECTOR 

Once  the  applied  loading  is  set,  the  program  asks  for  a  lamina  failure  criterion  and  stiffness 
reduction  model  to  be  used  in  the  failure  analysis: 


Enter  failure  criterion  and  reduction 

1:  Hill-Tsai 

2  :  Maximum  Stress 
3 :  Maximum  Strain 

4:  Tsai-Wu 

5 :  Hashin-Rotem 


model  { i . e . ,  1,1) 

1:  Parallel  Spring  Model 

2:  Parallel  Spring  Model (s) 

3 :  Sudden  Failure  Model 

4: 

5: 


B-10 


6: 
7  : 


Hashin 

all 


6: 

7: 


If  option  #7  is  entered  for  the  failure  criterion,  the  program  analyzes  the  chosen  laminate  with  all 
six  available  lamina  failure  criterion  (i.e.,  entering  7 , 1  does  an  analysis  using  all  six  lamina 
failure  criteria  along  with  the  Parallel  Spring  Model  for  stiffness  reduction).  See  the  theory 
section  for  explanation  of  reduction  models  #1  PSM  and  #2  PSMg.  If  option  #3  is  chosen  as  the 

reduction  model,  the  program  simply  calculates  the  “dominant”  ply’s  fiber  failure  strength  and 
ends  the  analysis  in  one  step. 

If  Tsai-Wu  is  chosen,  the  program  asks  for  the  interactive  coefficient  F,2  found  in  equation  4: 


Enter  F12  interaction  factor  for  Tsai-Wu: 
Return  =  default  of  0.0 


The  value  that  is  entered  is  then  multiplied  by  \I{XK').  This  is  done  since  those  terms  are 
frequently  found  in  F,2.  For  example,  if  0.5  is  entered,  =  \I{1XK')  which  would  reduce  Tsai- 
Wu  to  the  Hoffman  criterion  (not  offered  in  SALC).  As  explained  in  the  theory  section 
previously,  many  researchers  suggest  setting  F^2 zero. 

If  Hashin-Rotem  is  chosen,  the  program  asks  whether  or  not  to  include  the  interaction  factor  ii 
(as  described  in  equation  7): 

Include  mu  factor  with  Hashin-Rotem  criterion? 
l:yes  other: no 

If  yes,  the  program  prompts  for  the  value 

Enter  mu: 


Enter  a  value  for  p,  (typically  between  0.0  and  1.0). 

At  this  point,  the  program  is  capable  of  varying  the  angle  of  the  chosen  laminate’s  plies.  It  the 
loading  angle  j3  has  been  varied,  it  is  not  suggested  to  vary  any  ply  angles:  this  will  produce  a 
3-D  data  file  which  can  be  difficult  to  work  with.  The  user  is  prompted  with 

Do  you  want  to  vary  any  ply  angle  theta? 
l:yes  other: no 

If  yes,  the  program  asks 

Off-Axis  Laminate  Loading?  l:yes  other: no 
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Off-axis  laminate  loading  simply  rotates  the  entire  laminate  at  an  angle  0.  If  yes 

Enter  off-axis  range  and  step: 

Just  as  in  the  loading  case,  the  program  requires  an  angle  range  and  step.  For  example,  if  you 
entered  0,45,1  the  program  would  rotate  the  entire  laminate  in  1°  increments  from  0°  to  45° 
and  calculate  ultimate  laminate  failure  at  each  point. 

If  you  choose  “no”  to  off-axis  laminate  loading,  the  program  reminds  you  of  the  chosen  laminate: 

Here  is  the  original  laminate:  ([0/±45/90]s  forexample) 

Enter  angle  of  ply  to  vary. . . 


The  program  is  able  to  vary  only  one  ply’s  theta.  Enter  the  angle  of  the  ply  to  change.  Both  ± 
values  of  a  ply  may  be  changed  (i.e.,  ±60  or  ±45). 

Vary  both  signs  of  this  ply  (if  applicable)? 
l:yes  other: no 

This  allows,  for  example,  the  ±45°  plies  to  become  ±0°  plies  if  desired. 

Enter  range  of  theta  with  step: 


As  before,  enter  the  angle  range  and  step  of  the  ply  to  vary.  A  typical  example  of  this  feature 
would  be  to  change  a  [0/±45]s  laminate  to  [O/±0]s,  and  vary  0. 

Lastly,  the  program  asks 


Do  you  want  to  create  a  detailed  file  'analysis'? 
l:yes  other: no 

This  is  a  very  large  file  detailing  all  steps  in  the  failure  analysis.  Beware  that  this  file  can  grow 
very  large  (many  megabytes  if  not  careful)  and  slow  down  the  program.  If  either  the  loading 
angle  P  or  ply  angle  0  is  varied  over  10  steps  or  more,  it  is  suggested  that  analysis  not  be 
generated. 
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2.  OUTPUT 


The  program  creates  a  variety  of  output  files.  All  are  automatically  generated  with  the  exception 
of  analysis. 


out . file 


bistrs 


bistrn 


unstrs 


unstrn 


analysis 


Simple  output  file  which  gives  laminate  stresses  and  strains  at  each  ply 
failure  (including  ultimate  failure). 

Data  file  for  graphing.  Intended  for  plotting  biaxial  failure  envelopes. 
Contains  20  columns  in  the  order:  ply  angle  0,  loading  angle  p,  Hill-Tsai’s 
ultimate  <5^  and  Max.  Stress’  ultimate  <7^,  and  Max.  Strains’ 
ultimate  <3^  and  x^^;  Tsai-Wu’s  ultimate  <5^  and  x^^;  Hashin-Rotem’s 
ultimate  a^,  c^y  and  x^^;  Kashin’s  ultimate  a^,  and  x^y. 

Data  file  for  graphing.  Identical  to  bistrs  but  replaces  ultimate  stresses 
Oyy,  and  x^y  with  ultimate  strains  £^,  and  y^y. 

Data  file  for  graphing.  Intended  for  plotting  off-axis  laminate  failure  loads 
under  unidirectional  Nx.  Contains  eight  columns  in  the  order:  ply  angle  0, 
loading  angle  P,  Hill-Tsai’s  ultimate  a„,  Max.  Stress’  ultimate  <7^,  Max. 
Strains’  ultimate  (?„,  Tsai-Wu’s  ultimate  Hashin-Rotem’s  ultimate 
Hashin’s  ultimate  a^. 

Data  file  for  graphing.  Identical  to  unstrs  but  replaces  ultimate  stress 
with  ultimate  strain  e^. 

Very  detailed  file  containing  step  by  step  analysis  of  laminate  failure. 
Includes  all  stiffness  matrices,  stresses  and  strains  in  each  ply,  and  failure 
analysis  summary. 
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APPENDIX 


Following  are  examples  of  SALC’s  input  and  outputs.  In  order  they  are: 


mat . info 

available  material  input 

example . data 

material  system  file  explaining  order  of  constants 

aae555 . data 

available  material  system  (in  psi) 

exp . data 

available  material  system  (in  MPa) 

lam. info 

available  laminate  data  file 

exl 

example  of  input  for  [0/±45/90]s  laminate  under  tensile  Nx 

out . file 

output  from  exl 

analysis . 1 

detailed  output  from  exl 

ex2 

example  of  input  for  [90/0] s  laminate’s  thermal  stress  calculation 

analysis . 2 

detailed  output  from  ex2 

exp. data 
aae555 .data 


20e6 

EX 

1  -  Se6 

E2 

le6 

C312 

.29 

vl2 

-  005 

ply  tbic 

310o3 

X 

-3lO«3 

X' 

9e3 

Y 

-30e3 

Y' 

14e3 

S 

X 

xe 

1 

Xe' 

1 

Ye 

1 

Ye' 

X 

Se 

2e-6 

alpbalX 

lSe-€ 

alpba22 

0-0 

alphLal2 

AAES55 

Class  Not- 

Any  additional  c 

placed 

after  tbe 

label  to  each  value;  tHe  actua 
shoulci  not  contain  any  labels  . 


example-data 


■^^lues  in  p»i  , 

mater-ial  system  may  be 
s  -  This  printout  g^ives 
file  used  witb  SAic 


NOTE;  if  ultimate  strain 
calculate  tbe  value  if  a 
strain  is  calculate  from 
values  given . 


values  are  unavailable, 

of  1  is  used.  Tbe 
tbe  ultimate  stress  and 


SALC  will 
ultimate 
modulus 
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20e6 

1.5e6 

le6 

.29 

.005 

310e3 

-310e3 

9e3 

-30e3 

14e3 

1 

1 

1 

1 

1 

2e-6 

15e-6 

0.0 

AAE555  Class  Notes  Data:  stress  values  in  psi. 


Order  of  material  properties: 

El 

E2 

G12 

vl2 

ply  thickness 

X 

X' 

Y 

Y' 

S 

Xe 

Xe' 

Ye 

Ye' 

Se 

alphall 

alpha22 

aiphal2 
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153750 

11031 

6900 

.32 

.0129 

2171 

-2013 

€7 

-206.8 

110.3 

1 

1 

1 

1 

1 

2e-6 

15e-6 

0. 


exp.oata 


Tested  at  the  Composite  Materials  Lab  at  Purdue 
AS4/3501-6  Graphite  Epoxy  from  Hercules 
Fail  1993 

Stress  values  in  MPa 
Ply  thickness  in  cm 


University. 


strain  /  vl2 : 
.001  /  .333 
,002  /  ,320 
.003  /  .324 
.004  /  .317 
.005  /  .333 
.006  /  .327 
.007  /  .322 
.008  /  .318 
.009  /  .315 


Order  of  material  properties- 

El 

E2 

G12 

vl2 

ply  thickness 

X 

X' 

Y 

Y' 

S 

Xe 

Xe' 

Ye 

Ye' 

Se 

alphall 

alpha22 

alphal2 
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1 

O  degree  ply 

0.0 

8 

[0/+-45/90]s 

0.0 

45.0 

-45.0 

90.0 

90.0 

-45.0 

45.0 

0.0 

4 

[90/0]s 

90. 

0. 

0. 

90. 

6 

t90/0]2s 

90. 

0. 

90. 

0. 

0. 

90. 

0. 

90, 

6 

[0/+-60]s 

0 

60 

-60 

-60 

60 

0 

6 

[90/+-30]s 

90 

30 

-30 

-30 

30 

90 

6 

[0/+-451S 

0 

45.0 

-45,0 

-45.0 

45.0 

0 

6 

[+-45/901S 
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45 

--45 

90 

90 

-45 

45 

4 

[+-45]s 

45. 

-45, 

-45. 

45. 


10 

[0/90/90/90/90]s 

0 

90 

90 

90 

90 

90 

90 

90 

90 

0 
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Do  you  want  to  create  a  detailed  file  'analysis'? 
l;yes  other: no 


starting  Analysis  of:  [0/+-45/90]s 
Using  material  system:  aaeSSS.data 


out-file 


Hill-Tsai  w/  Parallel  Spring  Model 

Ply  Mode  SigmaXX  SigmaYY 


Beta:  0.00 

TauXY  eXX 


Theta : 
eYY 


0.00 

gXY 


90.00 

2 

51586.7 

0.0 

-45.00 

3 

62732.7 

0.0 

45.00 

3 

62732.7 

0.0 

0.00 

1 

103223.0 

o 

o 

Fiber  Failure 


0.0 

0.0065 

0.0 

0.0082 

0.0 

0.0087 

0.0 

0.0155 

-0.0019  0.0000 
-0.0023  0.0000 
-0.0027  0.0002 
-0.0051  0.0000 


Done ! 
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Enter  material  dlats  file  for  analysis 
1:  exp. data 
2:  aae555.data 

2 


ex2 
(#)  ... 


'  .  .A. 


Enter  #  of  laininate/ply  choice.  . . 
1 :  0  degree  ply 

2:  [0/+-45/90]s 

3:  [90/0]s 

4:  [90/0]2s 

5:  (0/+-60]s 

6:  [90/+-30]s 

7:  [0/+-45]s 

8:  [+-45/90]s 

9:  [+-45]S 

10:  [0/90/90/90/90]s 


Do  you  want  to  calculate  thermal  stresses? 
l:yes  other: no 

1 

Enter  delta  T: 

-250 

Enter  Mechanical  Loading  Envelope  #: 

1:  Nx/Ny  2:  Nx/Nxy  3:  Ny/Nxy  4:  Pure  Thermal 

4 

Enter  failure  criterion  and  reduction  model  (i.e.  1,1) 

1: 

2: 

3: 

4  : 

5: 

6: 

7  : 

Ir  1 

Do  you  want  to  vary  any  ply  angle  theta? 
l:yes  other:no 

2 

Do  you  want  to  create  a  detailed  file  'analysis"> 
l:yes  other :no 

1 

Starting  Analysis  of:  [90/0]s 
Using  material  system:  aae555.data 


Hill-Tsai 

Maximum  Stress 

Maximum  Strain 

Tsai“Wu 

Hashin-Rotem 

Hashin 

all 


rararrei  bpring  Model 
Parallel  Spring  Model (s) 
Sudden  Failure  Model 


Percent 

Complete : 

% 

0.0 

Percent 

Complete ; 

% 

0.0 

Percent 

Complete : 

% 

0.0 

Percent 

Complete : 

g. 

0.0 

Percent 

Complete : 

?: 

0.0 

Percent 

Complete : 

% 

0.0 

Percent 

Complete : 

?; 

0.0 

Percent 

Complete : 

% 

0.0 

Percent 
Done ! 

Complete : 

% 

0.0 
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.0000001873  0.0000046294  0.0000000000 

.0000000000  0.0000000000  0.0000500000 


ooo  oooooooooooo 


★  *•*■*•**★★****■***'*•*■*'★********★*******■*■■*■*★*■*■*■★'**★■*★*■*■★★•*■*****★******** 


* 

strength  Analysis  of  Laminated  Composties  (SALC) 


Version  1.1 

* 

June  1994 

* 

Brian  J.  Quinn  /  C.  T.  Sun  /  J.  Tao 

* 


■k 


********************************************************************** 


doubie  precision  thick,el(40),e2(40),gl2(40),vl2,theat(40), 
&fmr(3),ststrs(5),ststrn(5)Jangle,pi,temp,v21 
doubie  precision  angLang2,step,pangle,theatLtheat2,tstep, 
&  intfac,angvar,otheat(40),mu,guess,endstn(3),xfrnr(3) 
double  precision  dbstrs(6,3,3000),dbang(2,3000) 
double  precision  dbstrn(6,3,3000),alphl2(3),delt 
character*35  rmodel(l  0).fmodel(l  O),indatajamcho 
character*!  dothet panel 
integer  fmod,rmod,chfile,f  1  ,f2,lodsit,mufac,thanal 
common  panal 


***********  input  Data 


call  systemC'clear") 

open(unit=10,file='mat.info',access='sequentiar) 
write(6,*)"Enter  material  data  file  for  analysis  (#) ..." 
read(10,*)numfile 
do  15  k=Lnumfile 
read(10,10)indata 
10  format(a20) 

write(6,*)k,":  "Jndata 
15  continue 
read(5/)chfiie 
rewindO  0) 
do  25j=0,chfile 
read(l  0,20)indata 
20  format(a20) 

25  continue 
write(6,*) 

open(unit=15,fiIe=indataaccess='sequentiar) 
open(unit=25,file=’out.file') 
open(unit=40,file='analysis') 
open(unit=50Jle='lam.info',access='sequential') 
open(unit=70,file='bistrs') 
open(unit=7Lfile='bistrn') 
open(unit=80,file='unstrs’) 
open(unit=81  ,file='unstrn') 

read(15,*)  el(l),e2(l),gl2(l),vl2 
read(15,*)thick 
read(15,*)  (ststrs(i)J=l,5) 
read(l  5,*)  (ststrn(i)J=l  ,5) 

read(15,*)(alphl2(i),i=L3) 

call  CHOLAM(n,theat,lamcho) 

do27j=l,n 

el0)=e1(l) 

e20>e2(l) 

gl20)=gl2(l) 
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otheat(i)=theat(j) 
27  continue 
c 


delt=0.0 


write(6,*) 

write(6,*)"Do  you  wont  to  calculate  thermal  stresses?" 
write(6,*)"  1  ;yes  other:no" 
read(5,39)thanal 
if(thanal.eq.l)then 
write(6/)"Enter  delta  T:" 
read(5,*)delt 
endif 


if(ststrn(l  ).eq.  1 )  then 
ststmO  )=ststrsO  )/el  (1 ) 
endif 

if(ststm(2).eq.l)  then 
ststrn(2)=ststrs(2)/el  (1 ) 
endif 

if(ststrn(3).eq.l)  then 
ststrn(3)=ststrs(3)/e2(l ) 
endif 

if(ststrn(4).eq.l)  then 
ststm(4)=ststrs(4)/e2(l ) 
endif 

if(ststrn(5).eq.l)  then 
ststm(5)=ststrs(5)/gl  2(1 ) 
endif 


pi=3.14159265359 
intfac=0.0 
theotl  =0.0 
theat2=0.0 
tstep=l  .0 
angl=0.0 
ang2=0.0 
step=1 .0 

do  29  i= 1,6 
do29j=l,3 
do29k=l,1000 
dbstrs(j,j,l<)=0.0 
dbstm(i,j,l<)=0.0 
dbang(l,k)=0.0 
dbang(2,k)=0.0 
29  continue 

fmodel(l)='Hill-Tsai' 
fmodel(2)='Maximum  Stress' 
fmodel(3)='Maximum  Strain' 
fmodei(4)='Tsai-Wu' 
fmodel(5)='Hashin-Rotem' 
fmodel(6)='Hashin' 
fmodel(7)='ail' 
c 

rmodel(l  )='Parailei  Spring  Model' 
rmodel(2)='Parallel  Spring  Model(s)' 
rmodel(3)='Sudden  Failure  Model' 
rmodel(4)=' ' 
rmodel(5)=' ' 
rmodel(6)=' ' 
rmodel(7)=' ' 
c 

write(6,*) 
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write(6,*)"Enter  Mechanical  Loading  Envelope  #;" 
write(6,T  1 ;  Nx/Ny  2:  Nx/Nxy  3:  Ny/Nxy 
&"  4:  Pure  Thermal" 
read(5,39)lodsit 

if(lodsit.eq.4)  goto  100 

write(6,*)  ,  „  •, 

write(6/)"Enter  load-angle  range  and  step,  (i.e.  0,90,10) 
write(6,*) 

read(5,*)angl  ,ang2,step 
c 

write(6,*)  ,  , 

write(6,7'Enter  failure  criterion  and  reduction  model  (i.e.  1,1) 
write(6,*) 
do35j=l,7 

write(6, 30)j,f  modelO),j,  rmodelQ) 

30  format(i2,':  ',a35,2xi2,':  ',a35) 

35  continue 

read(5,  *)fmod,  rmod 

if(fmod.eq.7)  then 

fl=l 

f2=6 

else 

f  1  =fmod 
f2=fmod 
endif 
c 

if(fmod.eq.7.or.fmod.eq.4)  then 
write(6,*)  ,  ^  „ 

write(6,TEnter  FI  2  interaction  factor  for  Tsai-Wu: 

write(6,*)''  Return  =  default  of  0.0" 
read(5,37)intfac 

37  format(f6.3) 
write(6,*) 
write(6,38)intfac 

38  format('lnteraction  factor  =  ',f6.3) 
endif 


c  write(6,*) 

c  write(6,*)"Classify  final  failure:"  ,  „  ,  ,  ^ 

c  write(6,*)"  Return:  fiber  breakage  1 :  all  ply  s  damaged 

c  read(5,39)idamge 

39  format(il) 
c 

if(fmod.eq.7.or.fmod.eq.5)  then 
write(6,*)  ^  ^ 

write(6,Tlnclude  mu  factor  with  Hashin-Rotem  criterion? 
write(6,*)"  1  :yes  other:no" 
read(5,39)mufac 
if(mufac.eq.l )  then 
write(6,*)"Enter  mu:" 
read(5,*)mu 
guess=30000 

endif 

endif 

c 

41  write(6,*)  ,  ,  ^ 

write(6,*)"Do  you  want  to  vary  any  ply  angle  theta? 
write(6,*)"  1  :yes  other:no" 
read(5,40)  dothet 

40  format(al) 
if(dothet.ne.'T)  goto  100 

^  write(6,*)"Off-Axis  Laminate  Loading?  l:yes  otherino" 
read(6,39)ioffax 
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if(ioffax.eq.l)then 

write(6/)''Enter  off-axis  range  and  step:" 
read(5,*)theatl  ,theat2,tstep 
goto  1 00 
endif 


write(6,42)lamcho 

42  format('Here  is  the  originai  laminate:  ',a30) 
write(6,*) 


c 


c 


write(6,7'Enter  angle  of  ply  to  vary..." 
read(5/)angvar 

write(6, TVary  both  signs  of  this  ply  (if  applicable)?" 

write(6/)"  1  :yes  other:no" 

read(5,39)nsign 


write(6,*)  ^ 

write(6,*)"Enter  range  of  theta  with  step: 
read(5,*)theatl  ,theat2,tstep 
write(6,*) 


write(6,TDo  you  want  to  create  a  detailed  file  'analysis'?" 
write(6,*)"  1  :yes  other:no" 
read(5,n0)panal 
110  formatCal) 


******  Start  Analysis  ******************* 

write(25,n2)  lamcho 
write(6,112)  lamcho 

1 1 2  formatestarting  Analysis  of:  ',a30) 

write(25,*)"Using  material  system:  "Jndata 
write(6,*)"Using  material  system:  "Jndata 
write(6,*) 

if(panal.eq.'l')  then 
write(40,112)lamcho 

write(40,*)"Using  material  system:  "jndata 
write(40/) 
endif 

do  200  nfmod=f  1  ,f2 
c 

ihl 

c 

do  200  pangle=theatl,theat2,tstep 
if(dothet.eq.'r.and.ioffax.ne.l)then 
do  120k=l,n 
if(nsign.eq.l)then 

if(abs(otheat(k)).eq  .abs(angvar))  then 
if(otheat(k).lt.0.0)  then 
theat(k)=-l  .0*pangle 
else 

theat(k)=pangle 

endif 

endif 

0lS0 

if(otheat(k).eq.angvar)  then 
theat(k)=pangle 
endif 
endif 

120  continue 
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endif 

if(clothet.eq.'l  ’.and.ioffax.eq.l)  then 
do  135k=l,n 

theat(k)=otheat(k)+pangle 
135  continue 
endif 

150  do  200  Iangle=angl,ang2,step 
temp=langle*pi/l  80.0 

if(lodsit.eq.l)then 
fmr(l  )=dcos(temp) 
fmr(2)=dsin(temp) 
fmr(3)=0.0 

if(abs(langie).eq.90.0)  then 

fmr(l)=0.0 

endif 


c 

else  if(lodsit.eq.2)  then 
fmr(1  )=dcos(temp) 
fmr(2)=0.0 
fmr(3)=dsin(tennp) 
if(abs(langle).eq.90.0)  then 
fmr(l)=0.0 
endif 
c 

else  if(lodsit.eq.3)  then 
fnnr(l)=0.0 
fmr(2)=dcos(temp) 
fmr(3)=dsin(temp) 
if(abs(langle).eq.90.0)  then 
fmr(2)=0.0 
endif 


else  if(lodsit.eq.4)  then 
fmr(l)=0.0 
frnr(2)=0.0 
fmr(3)=0.0 
endif 


write(25,*) 

write(25, 1 30)fnnodel(nfmod),rmodel(rmod),langle,pangle 
write(25,125) 
write(25/) 

125  formate  Ply  Mode  SigmoXX  Sigma YY  TaoXV, 

&  7x/eXX  eYY  gXY^ 

if(panal.eq.'r)  then 

write(40, 1 30)fmodel(nfmod)jmodel(rmod),langle,pangle 
write(40,*) 

endif 

130  format(al5/w/',a26,'  Beta:’,f7.2,'  Theta:’,f7.2) 

call  STATUS(nfmod,fmod,theat1  ,theat2,tstep,angl ,ang2,step, 

&  pangle.langle) 

call  AN ALYS(el  ,e2,g  1 2,vl  2,v21  ,athiclctheat,fmr,ststrs,ststm, 

&  nfmod,rmocl,langle,pangleJntfac,idamgeJsmred,mu,guess,mufac, 
&  xfmr,endstn,alphl  Zdeltlodsitthanal) 


if(panal.eq.'l')  then 
write(40,*) 
write  (40,*) 
write(40,*) 


n**************  Q'f  oriolysis  ******************1* 

ii*************************************************ii 
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write(40,*) 

endif 

c 

do  190ii=l,3 

dbstrs(nfmod,iUj)=xfmr(ii) 

dbstrn(nfmod,ii,jj)=ends1n0i) 

190  continue 

dbang(l  Jj)=pangle 
dbang(2Jj)=langle 

ii=ii+i 

c 

200  continue 

^  call  DBASE(dbstrs,dbstradbang,  n,thiclcjj-1  ,ststrs) 

write(25,*) 

wiite(40,*) 

write(25,*)"Done!" 

write(6,TDone!" 

write(40,*)"Done!" 

c 

close(lO) 

close(15) 

close(25) 

close(40) 

close(50) 

close(70) 

close(71) 

close(80) 

close(81) 

c 

end 


subroutine  DBASE(dbstrs.dbstrn,dbang,n,thick,numel,ststrs) 


double  precision  dbstrs(6,3,3000),dbang(2,3000),thick 
double  precision  dbstrn(6,3,3000),ststrs(5),strnth(3) 

stmth(l)=ststrs(1) 
stmth(2)=ststrs(3) 
stmth(3)=ststrs(5) 

dol0nn=0J  .  ,, 

write(80+m,*)"blank  Theta  Beta  Hill-Tsai  Maximum-Stress  , 
8c"  Maximum-Strain  Tsai-Wu  Hashin-Rotem  Hashin" 
write(70+m,7'blank  Theta  Beta  Hill-Tsai  Hill-Tsai ", 

8c"  Hill-Tsai  Maximum-Stress  Maximum-Stress  Maximum-Stress  , 
8c"  Maximum-Strain  Maximum-Strain  Maximum-Strain", 

8c"  Tsai-Wu  Tsai-Wu  Tsai-Wu  Hashin-Rotem  Hashin-Rotem  , 

8c"  Hashin-Rotem  Hashin  Hashin  Hashin" 

10  continue 


c 


do  200  k=l,numel 
write(80, 1 00)  dbang(l  ,k),dbang(2, 
8c  (dbstrs(i,l,k)/n/thick,i=l,6) 

write(81 . 1 00)  dbang(l  ,k),dbang(2, 
8c  (dbstmC,  1  ,k),i=  1 ,6) 

write(70, 1 50)  dbang(l  ,k),dbang(2, 
8c  ((dbstrs(i,j,k)/n/ttiick,j=l  ,3),i=l  ,6) 

write(7 1 , 1 50)  dbang(l  ,k),dbang(2, 
8c  ((dbstm(i,j,k),j=l,3),i=l,6) 

100  format(2f8.2,2>c6fl5.5) 


k), 

k), 

k), 

k), 
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150  format(2f8.2,2x,18fl5.5) 

200  continue 
c 

return 

end 


subroutine  STATUS(nfmod,fmod,theatl  ,theat2,tstep,angl  ,ang2,step, 
&  panglejangle) 


double  precision  theatLtheat2,tstep,angl,ang2,step,pangle, 
&  langle,percntfacl,currnt,dem,num,temp 
integer  nfmod,fmod 

fool =1.0 

if(theatl  .ne.theat2)  then 
fad =1  +abs(theatl  -theat2)/abs(tstep) 
currnt=abs(pangle-theatl)/abs(tstep) 
endif 

if(angl  .ne.ang2)  then 
facl=l  +abs(angl-ang2)/abs(step) 
currnt=abs(langle-angl)/abs(step) 
endif 

if(fmod.eq.7)  then 
num=currnt+(nfmod- 1  )*fac  1 
dem=facl  *6.0 
else 

num=ourmt 

dem=facl 

endif 

do  200  x=l  .0,1 0.0 
temp=denn/10.0*x 
if(dint(nunn).eq.dint(temp))  then 
percnt=num/dem*l  00.0 
write(6, 1 00)percnt 

3  formatCPercent  Complete:  %  'f4.1) 
endif 

3  continue 

return 

end 


subroutine  CLPT(el  ,e2,gl  2,vl  2,v21  ,n,thick,theat,fmr,strsl  2, 

&  stml  2,lamstn,a,blow,alph  1 2,deltlodsit,thanal,thstl  2,thsnl  2) 

double  precision  a(3,3),b(3,3),d(3,3),thick,theat(40), 

&  fmr(3),el  (40),e2(40),gl  2(40),vl  2,v21  ,pi,temp,q(3,3),qbar(3,3) 
double  precision  tr,cl,c2,c3,c4,ql  I,ql2,q22,q66,sl,s2, 

&  s3,s4,qql  I,qql2,qq22,qq66,qql6,qq26,qt(3,3,40),tml(3,3,40) 
double  precision  tm2(3,3,40),t(40),z(40),zb(40),h,abd(6,6), 

&  lamstn(3),strsxy(3,40),strsl  2(3,40),stml  2(3,40),ainv(3,3) 
double  precision  alphl2(3),alphxy(3,40),delt 
&  thnTist(3),ntheim(3),thstxy(3,40),thstl  2(3,40),thsnl  2(3,40) 
integer  lodsitthanal 
character*!  panal 
common  panal 

pi=3. 141 59265359 
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blow=0.0 

do25i=L3 
nthermCO=0-0 
25  continue 

*********  Calculate  Thicknesses  and  Centroids  ******* 


h=n*thick 
z(l)=h/2. 
nl=n+1 
do  50  l=2,nl 


tCi-l)=thick 
z(l)=z(l)-thick*(i-l) 
zb(i-1)=(z(l)+z(i-l))/2. 
)  continue 


**********  Calculate  Q  Matrices 


***************** 
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qql6=(ql  l-ql2-2.*q66)*sl  *c3+(ql2-q22+2.*q66rs3*cl 
qq26=(ql  1-ql  2-2.*q66)*s3*cl  +(ql  2-q22+2.*q66)*sl  *c3 
qq66=(ql  1  +q22-2.*ql  2-2.*q66)*s2*c2+q66*(s4+c4) 

qt(l,l,k)=qqn 

qt(l,2,k)=qql2 

qt(2,Lk)=qql2 

qt(Lak)=qql6 

qt(3,l,k)=qql6 

qt(2,2,k)=qq22 

qt(2,3,k)=qq26 

qt(3,2,k)=qq26 

qt(3,3,k)=qq66 

c 

do  60  i= 1,3 
do60j=l,3 
qbar(i,j)=qt(i,j,k) 

60  continue 
c 

if(panal.eq.'l')then 

write(40,65)theat(k) 

65  formotCQbar  nnatrix:',f7.2) 
call  PMATRX(qbar) 
write(40,*) 
endif 

****  Form  Transformation  Matrices  ****** 


tml(l,l,k)=c2 

tml(Z2,k)=c2 

tm2(1,l,k)=c2 

tm2(2,2,k)=c2 

tml(l,2,k)=s2 

tml(2,l,k)=s2 

tm2(l,2,k)=s2 

tm2(2,l,k)=s2 

tml(l,3,k)=2,*sl*cl 

tm2(3,2,k)=2.*srcl 

tm2(3,l,k)=-tml(l,3,k) 

tml(2,3,k)=-tml(l,3,k) 

tm2(l,3,k)=sl*cl 

tml(3,2,k)=sl*c1 

tm2(2,3,k)=-tml(3,2,k) 

tm1(3,l,k)=-tml(3,2,k) 

tm2(3,3,k)=c2-s2 

tml(3,3,k)=c2-s2 

if(abs(theat(k)).eq.45,0)  then 
tml(3,3,k)=0.0 
tm2(3,3,k)=0.0 
endif 


*********  Caiculate 


Nthermal 


c 


75 


alphxy(l  ,k)=alph  1 2(1  )*c2+alph  1 2(2)*s2 

alphxy(2,k)=aiphl2(l)*s2+alphl2(2)*c2 
alphxy(3,k)=2.0'(alphl  2(l)-aiphl  2(2))*sl  *cl 


do75i=l,3 

do75j=l,3 

ntherm(i)=ntherm(i)+qt(i,j, 

continue 


k)*alphxy(i,  k)*t(k)*delt 


c 

100  continue 
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if(thanal.eq.1)then 

c 

if(panal.eq.'l')then 

write(40/)"Alpha's  in  x-y  coordinates;" 
write(40,7'  Ply  Alpha  x  Alpha  y  Alpha  xy" 
write(40/) 

write(40ri  05)  theat(k),alphxy(l  ,k),alphxy(Zk),alphxy(3,k) 
105  format(f7.1,2x,3el3.4) 

110  continue 

endif 
c 

if(panal.eq.'r)  then 
write(40/) 

write(40,*)'Thermal  loads;" 
write(4a*)"  Nx  Ny", 

&"  Nxy" 
write(40,*) 

write(40,80)ntherm(l),ntherm(2),ntherm(3) 

80  format(3fl6.2) 
write(40,*) 
endif 
c 

endif 


*********  Calculate  A,B,D  Matrices 


************ 


do200i=l,3 
do200J=l,3 
a(ij>0.0 
b(ij)=0.0 
d(i,j)=0.0 
200  continue 

do  250  i=l,3 


do250j=l,3 
do250k=l,n 
a(i  J)=a(i,j)+qt(i,j,  k)*t(k) 
b(i,j)=b(iJ)+qt(i,j,k)*t(k)*zb(k) 
d(U)=d(U)4it(ui)*(t(k)*zb(k)**2+t(k)‘*3/l^ 
250  continue 


call  matinv(a,ainv,blow) 
if(blow.eq.l  .0)  then 
if(panal.eq.'l')  then 
write(40/)"A  matrix;" 
call  PMATRX(a) 

write(40,*)  .  ,  „ 

write(40/)"The  A  matrix  has  become  singular;" 
write(40,*)'The  laminate  is  failed  by  matrix  cracking. 
write(40/) 
endif 


goto  500 
endif 


do300i=l,3 

do300j=l,3 

ii=i+3 

jj=j+3 

abd(ij)=a(i,j) 

abd(i,jj)=b(lj) 

abd(ii,j)=b(ij) 

abd(ii,jD=d(i,j) 
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300  continue 
c 

if(panal.eq.'l')  then 

write(40,*)"Effective  laminate  moduli:" 
write(40,325)"Ex  =  ",(a(l,l)*a(2,2)-a(L2)**2)/h/a(2,2) 
write(40,325)"Ey  =  ",(a(l,l)*a(2,2)-a(l,2)**2)/h/a(Ll) 
write(40,325)'Vxy  =  ",a(h2)/a(2,2) 
write(40,325)'Vyx  =  ",a(L2)/a(l,l) 
write(40,325)"Gxy  =  ”,a(a3)/h 
325  formotCa  1 0,f  1 5.4) 

write(40,*) 
endif 
c 

if(panal.eq.'l')  then 
write(40,*)"A  matrix:" 
call  PMATRX(a) 
write(40/) 

write(40/)"A  Inverse  matrix:" 
call  PIMTRX(ainv) 
write(40,*) 

write(40,*)"B  matrix:" 
call  PMATRX(b) 
write(40,*) 

write(40,*)"D  matrix:" 
call  PMATRX(d) 
write(40,*) 
endif 
c 

c********  Calculate  Laminate  Strains  (xy)  ******** 
c 

do350i=l,3 

lamstn(i)=0.0 

thrmstO)=0.0 

do350j=1,3 

lamstn(i)=lamstn(i)+ainv(i,j)*fmrO) 

thrmst(i)=thrmst(i)+ainv(iJ)*nthermO) 

350  continue 
c 

if(thanal.eq.  1 )  then 
if(panal.eq.'T)  then 
write(40,*) 

write(40/)"Effective  Laminate  Thermal  Coefficients:" 
write(40,*)"  Alpha  x  (bar)  Alpha  y  (bar)", 

&"  Alpha  xy  (bar)" 
write(40,*) 

write(40,390)thrmst(l)/delt,thrmst(2)/delt,thrmst(3)/delt 
390  format(3el7.4) 
write(40,*) 
endif 
endif 
c 

c*******  Calculate  Laminar  Stresses  (x-y)  ********** 

c 

do400l<=l,n 

do400i=l,3 

strsxyCi,k)=0.0 

ths5^(lk)=0,0 

do400j=l,3 

strsxy(i,k)=strsxy(i,k)  +  Iamstn0)*qt(i,j,k) 
thst)^(i,k)=thst^(i,k)  +  qt(i,j,k)*(thrmst(i)-delt*alphxyO,k)) 

400  continue 
c 

c**********  Calculate  Laminar  Stresses  and  Strains  (1-2)  ****** 
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do450k=l,n 

do450i=l,3 


strsl2(i,k)=0.0 

stml2(i,k)=0.0 

thstl20,k)=O.O 

thsnl2(i,k)=0.0 

do450j=l,3 

strsi  2(i,k)=strs1 2(i,k)+strsxyak)*tml  (i,j,k) 
stm12(lk)=stm12(i,k)+lannstna)*tm2(lj,k) 

thstl  2(i,k)=thstl  2(i,k)+ths1xyO-k)*tml  (i,j,k) 
thsnl  2(i,k)=thsn1 2(lk)+thrmstO)*tm2(i,j,k) 
450  continue 


if(panal.eq,'T)  then 
if(thanal.eq.l)then 
ifClodsit.eq.4)  then 

call  PTHERM(thrmstthstxy,thstl  2,thsnl  Zatheat) 

0lS@ 

call  Pl\/lECH(lamstn,strsxy,strsl  2,stm1 2,atheat) 
call  PTHERM(thmnst,thstxy,thstl  2,thsnl  Zatheot) 
endif 

0lS@ 

call  PMECH(lannstastrsxy,strsl  Zstrni  Zatheat) 
endif 
endif 


c 


500  return 
end 


subroutine  ANALYS(iel  ,ie2Jgl  Zvl  Ziv21  ,athick,^eatjfmr,ststrs, 

&  ststrn,fmodel,rmodel,iangle,pangle,intfac,idamge,ismred,mu,guess, 
&  mufac,xfmr,endstn,alphl  Zdeltlodsitthanal) 


double  precision  e1(40),e2(40),gl2(40),vlZthick,theat(40),fmr(3), 
&  ststrs(5),sts1m(5),strsl  2(3,40),strm  2(3,40),\^Liv21, scalar, 

$  scalar0,scalar00,fmr0(3),endstn0(3),endstn00(3),ratio,thec^ 

double  precision  crrat(3,40),iel  (40),scalar,mu,guess,thsnl  2(3,40), 
&  ie2(40),igl  2(40),ifmr(3),lamstn(3),a(3,3) 
double  precision  endstn(3),langle, 

&  pangle,intfac,it,xfmr(3),alph  1 2(3)^elt,thst  1 2(3,^) 
integer  frnodel,rmodel,fply,pmode(40),pply(40),rtemp,lodsit, 

&  ttianalfblyO 
character*!  panal 
common  panal 


Set  up  Data  for  Subroutine 


********** 


fload=l  .0 
nloop=l 
mode=0 
1^1y=0 
fply0=0 
lt=n*thick 
scalar00=l  .0 
scalarO=l  .0 

do  !0j=!,3 
fmr(i)=ifmrO) 
fmrO0)=fiTir(j) 
endstn00)=0. 
endstn00Q)=0. 
10  continue 
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v21=iv21 


do20j=l,40 

eia)=ie10) 

e2a)=ie2(i) 

gl2®=igl2© 

20  continue 

*******  Determine  Stresses  and  Strains  ************ 

0  call  CLPT(el ,e2,gl 2,vl  2,v21  ,n,thick,ttieat,fmr,strs1 2,strnl 2, 

&  lamstn,a,blow,alphl  2,deltlodsit,thanal,thstl  2,thsnl  2) 

if(lodsit.eq.4)  then 
mode=5 
goto  666 
endif 

******  Check  for  Unrealistic  Strains  (singular  A)  ***** 
if(blow.eq.1 .0)  goto  1000 

********  Call  Failure  Criterion  ******************* 

pmode(nloop)=mode 

pply(nloop)=fply 

if(fmodel.eq.l)  then 

call  HLTSAI(strsl  2,strnl  2,ststrs,ststran,fply,mode, 

&  crrat,scalar,scalarO,theatthstl  Zthanal) 
else  if(fmodel.eq.2)  then 
call  MXSTRS(strsl  2,stml2,ststrs,ststran,fply,mode, 

&  crratscalar,scalarO,theat,thstl  2,thanal) 
else  if(fmodel.eq.3)  then 
if(rmodel.eq.4)  then 
rtemp=ismred 
else 

rtemp=rmodel 

endif 

call  MXSTRN(strsl  Zstrnl  2,ststrs,ststm,n,fply,mode, 

&  crratscalar,scalarO,theatnloop, 

&  pmode,pply,rtemp,thsnl2,thanal) 
else  if(fmodel.eq.4)  then 
call  TSAIWU(strsl  2,strnl  Zststrs,ststran,fpIy,mode, 

&  crrat,scalar,scalarO,theat,lntfac,thstl  2,thanal) 
else  if(fmodel.eq.5)  then 
call  HASROT(stTsl  2,strnl  2,ststrs,ststran,fpiy,mode, 

&  crrat,scalar,scalarO,theat,mufac,mu,guessdhstl  2,thanal) 
else  if(fmodel.eq.6)  then 
call  HASHIN(strs1 2,stml  2,ststrs,ststm,n,fply,mode, 

&  crratscalar,scalarO,theat,thstl  2,thanal) 
endif 

pmode(nloop)=mode 

pply(nloop)=fply 

*******  Call  Stiffness  Reduction  Model  ************** 

if  (rmodel  .eq .  1 )  then 
call  PSMA(el  ,e2,gl  2,fply,theatn,mode) 
else  if(rmodel.eq.2)  then 
call  PSMB(el,e2,gl2,fply,theat,amode) 
else  if(rmodel.eq.3)  then 


B-40 


mode=l 

call  SFM(fply,theat,crrat,ascalar) 
endif 

250  continue 
c 

c********  Calculate  Strains  after  Ply  Failure  ********** 
c 

do  100j=1,3 

endstn(j)=lannstnO)*scalar 
100  continue 

Q 

********  Record  the  Maximum  Loads  and  Strains  ******** 


if(scalar00.1e.l  .O.and.scalarO.gt.l  .0)  then 

write(40,*)  x  j...  o  .  w 

write(40/)"****  Stress  Redistribution  at  the  Same  Loading 

wiite(40,*) 

theat0=theat(fply0) 

do333j=l,3 

fmrO0)=fmr0) 

endstnO0)=endstnOO(i) 

333  continue 
endif 


scalar00=scalar0 

fply0=fply 

cio444j=l,3 

endstnOO0)=endstn0) 

444  continue 

c******  Increase  FMR  and  Print  Info  ************** 
c 

do300i=l,3 

fmr(i)=fmr(i)*scalar 

xfmr(i)=fmr(i) 

300  continue 


call  PDATA(theatfply,mode,fmr,athiclcendstn) 
c 

if(panal.eq.'r)  then 
write(40,*) 

write(40,*)"Laminate  Strains:" 
write(40,*)"Before  ply  failure  After  ply  failure" 
write(40/) 

write(40,360)lamstn(l  ),endstn(l ) 
write(40,360)lamstn(2),endstn(2) 
write(40,360)lamstn(3),endstn(3) 

360  formcit(2x,fl  0.5, 15x,f  10.5) 
write(40,*) 
endif 

Q 

c****.  Check  for  Fiber  Failure/Unidirectional  Failure  **** 
c 

if(rmodel.eq.3)  then 
mode=4 
goto  1000 
endif 


if(mode.eq.l  .or.n.eq.l)  then 
goto  1000 
endif 


c 

if(idamge.eq.l)  then 
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itemp=0 
do  370  j=l,n 
do  365  k=l,nloop 

if(abs(theat(j)) .  eq  .a  bs(theat(pply(k))))  then 
itemp=itemp+l 
endif 

365  continue 

if(itemp.eq.O)  goto  380 
itemp=0 
370  continue 
c 

goto  1 000 
endif 
c 

380  continue 
nloop=nloop+1 
c 

if(panal.eq.'l')  then 
write(40,*) 

write(40,*)" - Reduced  laminate  stiffness... 

write(40,*) 

endif 


goto  50 


***********************************^ 


r************ 


C 
C 

1 000  continue 
c 

ratio=(abs(fmr(l))+abs(fmr(2))+abs(fmr(3)))/ 

$  (abs(fmr0(l  ))+abs(fmrO(2))+abs(fmrO(3))) 
if  (ratio.le.l.0)then 
dolll  j=l,3 

endstn0)=endstnO0) 

111  continue 
endif 

^  wrtte(40,T$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$3$$$$" 
write(40/)"$  $" 

write(40,*)"$  FINAL  FAILURE  STRAINS  $" 
write(40,*)"$  $" 

write(40,*)"$  $" 

write(40,222)endstn(l ) 
write(40,222)endstn(2) 
write(40,222)endstn(3) 
write(40,*)"$  $" 

write(k)/)"$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$S$$$$$$$$$$'' 
222  format(10x,fl0.5) 
c 

call  PSBEN  D(mode,  ratio.theatO,  endstn) 
c 

666  return 
end 


c - 

subroutine  HLTSAi(strsl  Zstrnl  2,ststrs,ststraafply,mode, 
&  crrat,scalar,scaiarO,theat,thstl  Zthanal) 


doubie  precision  strsl2(3,40),strnl2(3,40).xx(40),yy(40), 
&  ststrs(5Xststrn(5),c1  ,c2,c3,c4crrat(3,40),temp,scalar 
double  precision  theat(40),thstl  2(3,40),scalar0 
integer  fply,mode,thanal 
character*!  panal 
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common  ponol 

Determine  Sign  of  Strength  Constants  ******* 
call  STRNTH(strsl  2,ststrs,axxyy,theat) 

********  Find  Critical  Ratio  and  Load  ************ 


if(panal.eq.’l ')  then 
write(40,*) 

write(40,*)"Hill-Tsai  critical  ratios:" 
write(40,T  Ply  N-critical" 
write(40,*) 
endif 


do  100i=l,n 
cl=strsl2(l,0/xx0) 
c2=strsl2(2J)/yyO) 
c3=strsl  2(2J)/xx(i) 

c4=strsl2(3,i)/ststrs(5)  o 

crrat(lj)=dsqrt(l  .0/(cl  ;‘2^2*;2-cl  *c3+c4‘ 
iffpanal.eq.'l ')  write(40,50)theat(i),crrat(1j) 
50  format(f7.2,2x,f20.6) 


‘2)) 


100  continue 
c 


fply=l 

temp=crrat(l,l) 
do  150  i=l,n 

if(crrat(l/i).lt.temp)  then 
temp=crrat(lj) 
fply=i 
endif 

150  continue 

^  if(crrat(l  ,fply).lt.  1 .0)  then 

scalar0=l  ./crrat(l  ,fply) 
crrat(l,fply)=l  .0 
endif 

scalar=crrat(l  ,fply) 


c 

c*********  Find  Mode  of  Failure 
^  do200k=l,3 

strsl  2(l<,fply)=strsl  2(k,fply)*scalar 
200  continue 


cl  =strsl  2(1  ,fply)/xx(fply) 
c2=strsl  2(2,fply)/yy(fply) 
c3=abs(strsl  2(3,fply)/ststrs(5)) 
mode=l 

if(c2.gt.cl  .and.  c2.gt.c3)  then 
mode=2 
endif 

if(c3.gt.cl  .and.  c3.gt.c2)  then 
mode=3 
endif 


if(panal.eq.'l ')  then 

write(40/)  ,,  „ 

write(40,*)"Critical  ply  strength  ratios:" 
write(40,*y'  Fiber  Transverse", 
&"  Shear" 
write(40/) 

\Amte(40,250)cl  ,c2,c3 
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250  format(3fl4.4) 
write(40,*) 
endif 
c 

return 

end 


subroutine  MXSTRS(strsl  2,stm1 2,ststrs,ststran,fply,nnode, 
&  crrat,scalar,scalarO,theat,thstl  2,thanal) 


double  precision  strsl2(3,40),strnl2(3,40),xx(40),yy(40), 
&ststrs(5),ststrn(5),crrat(3,40),temp,scalar,theat(40) 
double  precision  thstl  2(3,40),scalar0 
integer  fply,mode,thanal 
character*]  panal 
common  panal 

*********  Determine  Sign  of  Strength  Constants  ******* 

call  STRNTH(strsl  2,ststrs,n,xx,yy,theat) 

********  Find  Critical  Ratio, Load, Ply, and  Mode  ******** 

if(panal.eq.'r)  then 
write(40,*) 

write(40,*)"Maximum  Stress  critical  ratios:" 

write(40,*)"  Ply  Fiber  Transverse", 

&"  Shear" 

write(40,*) 

endif 

do  100i=l,n 

crrat(l  ,i)=strs]  2(1  ,i)/xx(i) 

crrat(2,i)=strs]  2(2,i)/yy(i) 

crrat(3,  i)=abs(strs  1 2(3,  i))/ststrs(5) 

if(panal.eq.’l ')  write(40,50)  theat(i),crrat(l,i),crrat(2,i), 

&  crrat(3,i) 

50  format(f7.2,2x,3fl5.8) 

100  continue 
c 

temp=crrat(l,l) 

1ply=l 

mode=l 

do200i=l,n 

do200j=l,3 

if  (crratO,i).gt.temp)  then 
temp=crrat(j,i) 
fply=i 
mode=j 
endif 

200  continue 
c 

if(crrat(mode,fply).gt.l  .0)  then 
scalarO=crrat(mode,fply) 
crrat(mode,fply)=l  .0 
endif 

crrat(mode,f  ply)=  1 . 0/crrat(mode,fply) 
scalar=crrat(mode,fply) 
c 

if(panal.eq.'l ')  write(40,*) 

return 

end 
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subroutine  MXSTRN(strsl  Zstrnl  2,ststrs,ststrn,n,fply,moda 

&  crrat,scalar,scalar0,theatnloop,pmode,pplyjeclmodthsnl2,thanal) 


double  precision  strsl2(3,40),strnl2(3,40),xx(40),yy(40), 
&ststrs(5).ststm(5),crrat(3,40),temp,scalar,theat(40), 
&thsnl2(3,40),scalar0  „ 

integer  fply,mode,pmode(40),pply(40),redmod,thanal 
character*!  panal 
common  panal 

Determine  Sign  of  Strength  Constants  ******* 
call  STRNTH(strnl  Zststrn,axx,yy,theat) 

********  Find  Critical  Ratio,Load,Ply,and  Mode 


if(panal.eq.'1 ')  then 
write(40,*) 

write(40,*)"Maximum  Strain  critical  ratios: 

write(40,*)"  Ply  Fiber  Transverse", 
&"  Shear" 
write(40,*) 
endif 


do  lCX}i=l,n 

crrat(l  ,i)=stTnl  2(1  ,D/xx(i) 
crrat(2,i)=strnl  2(2,i)/yy(i) 
crrat(3,D=abs(strnl2(3,i))/ststm(5) 
c  if(panal.eq.’l  ’)  write(40,50)  theat(i),crrat(l  ,i),crrat(2,i), 
c  &  crrat(3,i) 

50  format(f7.Z2x,3fl5.8) 

100  continue 


temp=0.0 

fpiy=l 

mode=l 

do200i=l,n 

do200j=l,3 

do  150  k=l,nloop 

if  (theat(i).eq.theat(pply(k)).and.pmode(k).ne.O)  then 
if0.ne.l)then 

if(redmod.eq.2.and.pmode(k).eq.2)  then 
crrat(2,0=O.O 

CTrat(ZD=0.0 

crrcit(3,D=0.0 

endif 

endif 

endif 

150  continue 
c 

if  (crratO,i).gt.temp)  then 
temp=crratQ,i) 
fply=i 
mode=j 
endif 

200  continue 


if(panal.eq.'l')  then 
do300i=l,n 

write(40,50)  theat(i),crrat(l  ,i),crrat(2,i), 
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&  crrat(3,i) 

300  continue 
write(40,*) 
endif 

if(ciTat(mode,fply).gt.l  .0)  then 
scalarO=crrat(mode,fply) 
crrat(mode,fply)=l  .0 
endif 

crrat(mode,fply)= 1 .0/crrat(mode,fply) 
scalar=crrat(nnodedply) 
c 

return 

end 


subroutine  TSAIWU(strsl  Zstrnl  2,ststrs,ststraafply,mode, 
&  crrat,scalar,scaiat0,theat,intfac,thst1 2,thanal) 


double  precision  strsl2(3,40),strnl2(3,40),xx(40),yy(40), 
&ststrs(5),ststm(5),crrat(3,40),temp,scalar,theat(40) 
double  precision  fl,f2,fl  I,f22,fl2,f66,cl,c2,c3,c4,c5,c6, 
&  a,b,c,prodJntfac,thstl  2(3,40),scalar0 
integer  fply,mode,thanal 
character*!  panel 
common  panal 

*******  Find  Critical  Ratio  and  Load  ************ 

if(panal.eq.'l')then 
write(40,*) 

write(40, 1 0)intfac 
10  format('interaction  factor  =  ',f6.3) 
write(40,*) 

write(40/)'Tsai-Wu  critical  ratios:" 
write(40/)"  Ply  N-critical" 
write(40,*) 
endif 

f  1 = 1 .0/ststrs(l )  +  1 .0/ststrs(2) 
f2=l  .0/ststrs(3)  +  1 .0/ststrs(4) 
fll=-1.0/ststrs(!)/st5trs(2) 
f22=-l  .0/ststrs(3)/ststrs(4) 
f  1 2=intfacMstrs(l  )/ststrs(2) 
f66=1.0/ststrs(5)**2 


do  100i=l,n 

cl=fl*strsl2(l,i) 
c2=f2*strsl2(2J) 
c3=fll*strs  12(1 ,0**2 
C4=f22*strsl2(2,i)**2 
c5=2.0*fl  2*strsl  2(1  ,i)*strsl  2(2,i) 
c6=f66*strsl2(3,i)**2 

a=1.0 

b=(cl  +c2)/(c3+c4+c5+c6) 
c=-l  .0/(c3+c4+c5+c6) 
prod=b*b-4.0*a*c 
crrat(l  ,i)=(-b+dsqrt(prod))/2/a 
if(panal.eq.'r)  write(40,50)theat(i),crrat(l,i) 
50  format(f7.2,2>cf20.6) 
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TOO  continue 
fply=1 

temp=crrat(l,l) 
do  150i=l,n 

if(crrat(l,i).lt.temp)  then 
temp=crrat(l,i) 
fi^y=i 
endif 

150  continue 

if(crrat(l  ,fply).lt.  1 .0)  then 
scalarO=l  ./crrat(l,fply) 
ciTat(l,fply)=1.0 
endif 

scalar=crrat(l  ,fply) 

Q 

^*********  Deternnine  Sign  of  Strength  Constants 

Q 

call  STRNTH(strsl  2,ststrs,n,xx,yy,theat) 
c 

c*********  Find  Mode  of  Failure 
^  do200k=l,3 

strsl  2(k,fply)=strsl  2(lcfply)*scalar 
200  continue 

cl  =strsl  2(1  ,fply)/xx(fply) 
c2=strsl  2(2,1ply)/yy(fply) 
c3=abs(strs  1 2(3,fply)/ststrs(5)) 
nnode=l 

if(c2.gt.cl  .and.  c2.gt.c3)  then 
nnode=2 
endif 

if(c3.gt.cl  .and.  c3.gt.c2)  then 
nnode=3 
endif 
c 

if(panal.eq.'l ')  then 

write(40,*)  „ 

write(40,*)"Critical  ply  strength  ratios:" 
write(40,*)"  Fiber  Transverse", 

&"  Shear" 
write(40,*) 

\vrite(40,250)cl  ,c2,c3 
250  fomnat(3fl4.4) 
write(40,*) 
endif 
c 

return 
end 
c 


subroutine  HASROT(strsl2,stml2,ststrs,ststrn,afply,rTiode, 

&  crrat,scalar,scalarO,theat,mufac,mu,guess,thstl  zthanal) 


double  precision  strsl  2(3,40),strnl  2(3,40),xx(40),yy(40y 
&  ststrs(5),ststm(5),cl  ,c2,crrat(3,40),temp,scalar  scalarO 
double  precision  theat(40),nnu,sol,guess,thstl  2(3,40) 
integer  fply,mode,mufac,thanal 
character*!  panal 
common  panal 
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*********  Determine  Sign  of  Strength  Constants  ******* 
call  STRNTH(strsl  2,ststrs,n,xx,yy,theat) 

********  Find  Critical  Ratio  and  Load  ************ 


if(panal.eq.'T)  then 
write(40,*) 

write(40,*)"Hashin-Rotem  critical  ratios:" 

write(40,*)"  Piy  Fiber  Matrix" 
write(40,*) 
endif 

do  100i=l,n 
crrat(l  ,i)=strsl  2(1  J)/xx(i) 
if(mufac.ne.l  .or.strsl2(2,i).ge,0.0)  then 
cl=strs12(2,i)/yy(i) 
c2=abs(strsl  2(3,i)/ststrs(5)) 
crrat(2,i)=dsqrt(cl  **2+c2**2) 
crrat(3,i)=crrat(2,i) 

0lS0 

call  NEWTON(strsl2(2J),strsl  2(3,i),mu,ststrs(5),yy(i). 
&  guess,sol) 
crrat(2,i)=1 .0/sol 
endif 


c 

if(panal.eq,'l ')  write(40,50)  theat(i),crrat(Li),crrat(2,i) 
50  format(f7.22x,2fl5.8) 

100  continue 
c 


temp=crrat(l,l) 
fply=l 
mode=l 
do  150i=l,n 
do  150j=l,2 

if  (crrat(j,i).gt.temp)  then 
temp=crrat(]J) 
fply=i 
mode=j 
endif 

150  continue 
c 

if(crrat(mode,fpiy).gt.l  .0)  then 
scalarO=crrat(mode,fply) 
crrat(mode,fply)=l  .0 
endif 

crrat(mode,fply)= 1 .0/crrat(mode,fpiy) 
scalar=crrat(mode.fply) 


c*********  Find  Mode  of  Failure 
^  do200k=l,3 

strsl  2(k,fply)=strsl  2(k,fply)*scalar 
200  continue 
c 

cl  =strsl  2(1  ,fply)/xx(fply) 
c2=strsl  2(2,'fply)/yy(fply) 
c3=abs(strsl  2(3,fply)/ststrs(5)) 
c 

if(mode.eq.2)  then 
if(c3.gt.c2)  mode=3 
endif 


c 
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if(panal.eq.'l ')  then 
write(40,*) 

write(40,*)"CriticaI  ply  strength  ratios:" 
write(40,*)"  Fiber  Transverse", 
&"  Shear" 
write(40,*) 
write(40,250)cl  ,c2,c3 
250  fomnat(3fl4.4) 
write(40,*) 
endif 
c 

return 

end 


subroutine  HASHIN(strs1  Zstrnl  2,ststrs,ststrn,n,fply,mode, 
&  crrat,scalar,scalarO,theat,thstl  2,thanal) 


double  precision  strsl  2(3,40),strnl  2(3,40),xx(40),yy(40), 

&  ststrs(5),ststrn(5),cl  ,c2,c3,crrat(3,40),temp,scalar,scalar0 

double  precision  theat(40),thstl  2(3,40) 

integer  fply,mode,thanal 

character*!  panal 

common  panal 

*********  Determine  Sign  of  Strength  Constants  ******* 
call  STRNTH(strs1 2,ststrs,n,xx,yy,theat) 

********  Find  Critical  Ratio  and  Load  ************ 


c 


c 


if(panal.eq,'r)  then 
write(40,*) 

write(40,*)"Hashin  critical  ratios:" 

write(40,*)"  Ply  Fiber  Matrix" 
write(40,*) 
endif 


do  100i=l,n 

cl=strsl2(l,i)/xx0) 

c2=strsl2(2,i)/yy(i) 
c3=abs(strsl  2(3,i)/ststrs(5)) 


if(strsl  2(1  ,D.gt.0.0)  then 
crrat(l  ,i)=dsqrt(cl  **2+c3**2) 


else 

crrat(l,i)=cl 

endif 


crrat(2,i)=dsqrt(c2**2+c3**2) 

c 

if(panal.eq.'l ')  write(40,50)  theat(i),crrat(l,i),crrat(2,i) 
50  format(f7.Z2x,2fl5.8) 

ICXD  continue 
c 

temp=crrat(l,l) 
fjoly=l 
mode=! 
do  150i=l,n 
dol50j=l,2 

if  (crratO-i).gt.temp)  then 
temp=crrat(j,i) 
fply=i 
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mode=j 

endif 

150  continue 
c 

if(crrat(mode,tply).gt.l  .0)  then 
scalarO=crrat(mode,fply) 
crrat(mode,fply)=l  .0 
endif 

crrat(mode,fply)=1.0/crrat(mode,fply) 

scalar=ciTat(mode,fply) 

c 

c*********  Find  Mode  of  Failure  ***** 
c 

do200l<=L3 

strsl  2(k,fply)=strsl  2(l<,fply)*scalar 
200  continue 
c 

cl  =strsl  2(1  ,fply)/xx(fply) 
c2=strsl  2(2,1ply)/yy(fply) 
c3=abs(strsl  2(3,fply)/ststrs(5)) 
c 

if(mode.eq.l)  then 
if(c3.gt.cl)  mode=3 
endif 

if(mode.eq.2)  then 
if(c3.gt.c2)  mode=3 
endif 


*★★★***★**** 


if(panal.eq.'r)  then 
write(40,*) 

write(40,*)"Critical  ply  strength  ratios:" 
write(40,*)''  Fiber  Transverse", 

&"  Shear" 
write(40,*) 
write(40,250)c  1  ,c2,c3 
250  format(3fl4.4) 
write(40,*) 
endif 
c 

return 

end 

c 

c - 

subroutine  STRNTH(value,stren,  n,xx,yy,theat) 

c - 

double  precision  value(3,40),stren(5),xx(40),yy(40),theat(40) 
character*!  panal 
common  panal 
c 

do  100i=l,n 
c 

if(value(l,i).gt.0.0)  then 
xx(i)=stren(l ) 

©Is© 

xx(i)=stren(2) 

endif 

c 

if(value(2,i).gt.0.0)  then 
yy(i)=stren(3) 
else 

yy(i)=stren(4) 

endif 

c 

100  continue 
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r 


return 

end 


subroutine  PSMA(el  ,eZgl  2,fply,theatn,mode) 


double  precision  el(40),e2(40),gl2(40),theat(40) 
integer  fply,mode 

do  100k=1,n 

if(abs(theat(fply)).eq.abs(theat(k)))  then 
if(theat(fply)  .eq.  theat(k))  then 
e2(k)=0.0re2(k) 
gl2(k)=0.0rgl2(k) 
endif 

0  continue 

return 

end 


subroutine  PSMB(el,e2,gl2,fply,theatn,mode) 


double  precision  el(40),e2(40Xgl2(40),theat(40) 
integer  fply,mode 

dol00k=l,n 

if(abs(theat(fply)).eq.abs(theat(k)))  then 
if(theat(fply)  .eq.  theat(k))  then 
if(mode.eq.3)  then 
e2(k)=0.0re2(k) 
gl2(k)=0.0rgl2(k) 
gIsg 

e2(k)=0.0re2(k) 

endif 

endif 

0  continue 

return 

end 


subroutine  SFM(fply,theat,  crrat,  ascalar) 


double  precision  theat(40),crrat(3,40),temp,scalar 
integer  tply 

temp=crrat(l,l) 
fply=l 

do  lCXDi=l,n 

write(40,*)crrat(l  J),temp,fply 
if(crrat(lj).gt.tennp)  then 
temp=crrat(l,i) 
fpiy=i 
endif 
too  continue 
c 

scalar=temp 
if(scalar.lt.l  .0)  then 
scalar=l  .0/scalar 
endif 
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return 

end 


subroutine  PSBEND(mode,ratio,theatO,endstn) 

c - 

c 

double  precision  endstn(3),ratiodheat0 
integer  mode 
if(ratio.le.  1.0)  then 
write(25,5)theat0 

5  formate***  The  failure  of  ',F6.1,'  deg  ply  is  catastrophic!') 
write(25, *)"***  FINAL  FAILURE  STRAINS:  " 
write(25,10)(endstnO),  i=l,3) 

10  format(51x3f9.4) 
else  if(mode.eq.l)  then 
write(25,*)"Fiber  Failure" 
else  if(mode.eq.5)  then 

write(25,*)''Only  Thermal  Stresses  Calculated  -  No  Failure", 
&"  Analysis  Performed" 
else  if(mode.eq.4)  then 
write(25,*)"Sudden  Failure" 
else  if(mode.eq.2.or.mode.eq.3)  then 
write(25,*)"Matrlx  Failure" 
else 

write(25,*)"lnfinite  Loop  or  Error  Occurred" 
endif 
c 

return 

end 

c 

c - 

subroutine  PDATA(theat,fply,mode,fmr,n,thick,endstn) 

c - 

c 

double  precision  theat(40),fmr(3),thicl<,endstn(3) 
integer  fply,mode 
character*!  panal 
common  panal 


write(25, 1 0)  theat(fply),mode,fmr(l  )/n/thick,fmr(2)/n/thick, 
&  fmr(3)/n/thick,(endstn(i),  1=1,3) 

10  format(f7.2,i4,2x3fl2.1,2x,3f9.4) 
c 

if(panal.eq.'r)  then 

write(40,*)"Ply  Failure:" 

write(40,*)"  Ply  Mode  N(xy)  Sigma(xy)" 
write(40,*) 

write(40,20)  theat(fply),  mode,fmr(l  ),fmr(l  )/n/thick 
write(40,30)  fmr(2),fmr(2)/n/thick 
write(40,30)  fmr(3),fmr(3)/n/thick 
20  format(f7.2,i4,2x,2fl2.1) 

30  format(13x2fl2.1) 
endif 
c 

return 

end 

c 

c - 

subroutine  PMATRX(m) 

c - 

c 

double  precision  m(3,3) 
c 
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subroutine  PIMTRX(m) 


write(40, 1 0)m(l ,  1  ),m(l  ,2),m(l  ,3) 
write(4a  1 0)nn(2, 1  ),m(2,2),m(Z3) 
write(40, 1 0)m(3, 1  ),m(3,2),m(a3) 
10  formcit(17>c3f20.4) 
c 

return 
end 
c 
c- 

c - 

c 

double  precision  m(3,3) 

^  write(40, 1 0)m(l ,  1  ),m(l  ,2),  m(l  ,3) 
write(40, 1 0)m(2, 1  Xm(Z2),m(2,3) 
write(40, 1 0)m(3, 1  ),m(3,2),m(3,3) 
10  formotO  5x3^0. 10) 
c 

return 

end 


c 

c- 


subroutine  PMECH(lamstn,strsxy,strsl  2,stml  2,n,theat) 


double  precision  lamstn(3),strsxy(3,40),strs  12(3,40), 
&stml2(3,40),ttieat(40) 
c 

write(40,*) 

write(40,*)"Mechanical  Stresses  and  Strains  on  each  ply:" 
write(40,*)"  Ply  Strain(xy)  Strain(12)  ", 

&"  Stress(xy)  Stress(l  2)" 
do  100i=l,n 
write(40,*) 
do  100j=l,3 

write(40, 1 0)theatO),lamstnO),strnl  20,i),strsxy0',D'Strsl  20, i) 

10  format(f7.Z2x,2fl2.7,3x2fl5.3) 

100  continue 
c 

200  return 
end 
c 

c - 

subroutine  PTHERM(thrmst,thstxy,thstl  2,thsnl  Zatheat) 


double  precision  thrmst(3),thstxy(3,40),thstl  2(3,40), 
&thsnl2(3,40),theat(40) 
c 

write(40,*) 

write(40,*)"Thermal  Stresses  and  Strains  on  each  ply:" 
write(40,*)"  Ply  Strain(xy)  Strain(12)  ", 

&"  Stress(xy)  Stress(12)" 
do  100i=l,n 
write(40,*) 
do  100j=l,3 

write(40, 1 0)theatO),thrmstO),thsnl  2Q,i),thstxyO,D-fr>stl  20,i) 
10  format(f7.Z2x2fl2.7,3x2fl5.3) 

100  continue 
c 

200  return 
end 
c 
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c - 

subroutine  MATINV(a,ai,blow) 
c - 


c 

c 


c 


double  precision  a(3,3),ai(3,3),deta 


deta  =  a(Ll)*(  a(2,2)*a(3,3)-a(2,3)*a(3,2) ) 

#  -  a(l,2)*(  a(2,l)*a(3,3)-a(2,3)*a(3,l) ) 

#  +  a(L3)*(  a(2,l)*a(3,2)-a(2,2)*a(3,l) ) 


write(40,*)''deta='',deta 
if(abs(deta).lt.lOOO)  then 
blow=l  .0 
goto  1 00 
endif 


c 


c 


c 

c 


c 

100 


deta=l  .0/deta 

ai(l,l)  =  (a(2,2)*a(3,3)  -  a(2,3)*a(3,2))*deta 
ai(l,2)  =-(0(2,1  )*a(3,3)  -  a(3,l)*a(2,3))*deta 
ai(l,3)  =  (a(2,l)*a(3,2)  -  a(2,2)*a(3,l))*deta 
ai(2,l)  =-(a(l,2)*a(3,3)  -  a(l,3)*a(3,2))*deta 
ai(2,2)  =  (a(l,l)*a(3,3)  -  a(l,3)*a(3,l))*deta 
ai(2,3)  =-(a(l,l)*a(3,2)  -  a(1.2)*a(3,l))*deta 
ai(3,l)  =  (a(l,2)*a(2,3)  -  a(2,2)*a(l,3))*deta 
ai(3,2)  =-(a(l,l)*a(2,3)  -  a(2,l)*a(l,3))*deta 
ai(3,3)  =  (a(l,ira(2,2)  -  a(l,2)*a(2,l))*deta 
ai(2,l)  =  ai(l,2) 
ai(3,l)  =  ai(l,3) 
ai(3,2)  =  ai(2,3) 


return 

end 


c 

c - 

subroutine  CHOLAM(n,theat.lamcho) 

c - 

c 

double  precision  theat(40),theata(20,40) 
character*30  lamtyp(20),lamcho 
integer  choice,  nplys(20) 
c 

read(50,*)nlams 

c 

do  200  j=l,nlams 
read(50,*)nplys(i) 
do  100  i=l,nplysO)+l 
if(i.eq,l)then 
read(50, 1 0)lamtyp(j) 

1 0  format(a30) 

0lS@ 

read(50,*)theataO,i-l ) 
endif 

100  continue 
200  continue 
c 

write(6,*)"Enter  #  of  iominote/ply  choice..." 

do  300j=l,nlams 

write(6,250)j,lamt/p(j) 

250  formotoa':  ',a30) 

300  continue 
c 

read(5,*)choice 

n=nplys(choice) 

lamcho=lamtyp(choice) 
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GO  GO 


do400k=l,n 

theat(k)=theata(choice,  k) 
400  continue 
c 

return 

end 


subroutine  NEWTON(sig22,tau,mu,S.Y,guess,soI) 


double  precision  sig22,taamaS,Y,solab,c,d,g,num,dem,ans, 
&  delguess 

a=mu**2*sig22*M/Y**2 
b=-2.0*S*nnu*sig22**3/Y**2 
c=S**2*sig22**2/Y**2+tau**2-mu**2*sig22**2 
d=2.0*S*mu*sig22 

g=guess 
ans=0.0 
del=g-ans 

do  100i=l,50 

num=a*g*M+b*g**3+c*g**2+d*g-S**2 
dem=4*a*g**3+3*b*g**2+2*c*g+d 
ans=g-num/dem 
del=g-ans 
g=ans 

if(abs(del).lt.0.01)  goto  150 
100  continue 
c 

150  sol=ans 
c 

return 

end 
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